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FINAL  REPORT 
on 

TRANSIENT  INTERACTION  OF  ELECTROMAGNETIC  PULSES  IN  DIELECTRICS 

AND  MICROWAVE  BIOPHYSICS 

by 

Theodore  C.  Guo  and  Wendy  W.  Guo 

The  Catholic  University  of  America 
Washington,  D.  C.  20064 


1.  BACKGROUND 

"  Due  to  recent  progress  in  developing  equipments  that  can  generate  short  microwave  and 
millimeter  wave  pulses,  there  has  been  an  increasing  proliferation  of  microwave  pulse 
transmitters,  some  with  short  pulse  width  (0.1  microsecond)  and  extremely  high  intensity 
(100-1000  megawatts).  Microwave  pulse  transmitters  are  used  extensively  by  the  military  for 
communication  and  remote  control;  using  microwave  pulses  as  directive  energy  weaponry  and  as 
means  of  transporting  energy  has  also  been  contemplated.  Electromagnetic  pulses  (EMP)  are 
also  emitted  in  nuclear  blasts  and  from  EMP  simulators.  All  this  production  of  microwave 
pulses  affects  the  operation  of  military  personnel  in  non-combat  environment  as  well  as  in 
battle  fields.  Therefore  minimizing  microwave  damage  is  central  to  successful  operations  of  all 
military  units.  Understanding  basic  interactions  between  microwave  pulses  and  dielectric 
materids  will  contribute  greatly  to  the  protection  of  human  subjects  from  microwave  damage 
and  to  the  development  of  preventive  measure.  . 

Until  recently  most  analyses  on  microwave  effects  have  been  based  on  continuous  wave 
approach.  Recognizing  the  importance  of  basic  understanding  of  the  interactions  between  short 
radiation  pulses  and  dielectric  materials,  and  its  potential  application  to  radiation  hazards  and 
radiation  treatment  to  biological  materials,  the  Walter  Reed  Army  Institute  of  Research 
(WRAIR)  started  a  program  at  the  Catholic  University  of  America  (CUA),  beginning  6/1/85,  to 
study  the  transient  interaction  of  electromagnetic  pulses  in  dielectrics  and  microwave 
biophysics.  The  program  was  administrated  by  the  Office  of  Naval  Research  (ONR).  Originally 
the  program  was  intended  for  three  years  to  end  on  5/30/88.  However,  due  to  budget 
constraint,  the  third  year  funds  \\'as  not  provided  and  the  program  prematurely  ended  on 
9/30/87.  This  document  repons  the  program  progress  as  of  the  end  of  funding  and  constitutes 
a  final  report  of  the  program.  The  main  body  of  the  report  provides  a  progress  summary  with 
technical  details  given  in  the  appendices  which  are  publication  reprints. 


2.  PROGRAM  OBJECTIVES 

The  long-term  objectives  of  this  program  are  as  follows; 

2.1.  To  understand  the  basic  physics  of  interaction  between  microwave  radiation  pulses  and 
dielectric  materials,  including  biological  systems. 
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2.2.  To  define  and  derive  characteristic  physical  quantities  relating  to  the  absorption  and 
energy  conversion  of  short  radiation  pulses  in  dielectric  materials,  especiily  water 
dominated  media,  such  as  biological  subjects,  and  to  define  biological  hazard  of  microwave 
pulses. 

2.3.  To  understand  the  basic  interactions  of  microwave  radiation  pulses  with  biological 
subjects,  and  to  understand  the  mechanisms  of  secondary  interactions  resulting  from 
microwave  pulses,  to  define  and  to  quantify  physical  and  biological  parameters  of  these 
interactions  with  respect  to  dielectric  properties  and  pulse  parameters,  such  as  carrier 
frequency,  pulse  width,  and  peak  power. 

2.4.  To  develop  dosimetric  techniques  and  hazard  specification  that  are  applicable  for  transient, 
near-field,  and  high-field  regimes,  and  to  develop  applications  of  microwaves,  including 
imagery  and  target  organ  andysis. 


3.  SCOPE  OF  WORK 

The  program  consists  of  the  following  three  subjects: 

3.1.  Transient  Dielectric  Relaxation  and  Absorption 

This  task  is  to  study  the  interaction  of  microwave  radiation  pulses  with  dielectrics  and 
biological  systems,  and  the  material  response  to  ultra-short  pulses  in  the  transient  regime. 

3.2.  Microwave  to  Acoustic  Energy  Conversion  Induced  by  Microwave  Pulses 

This  task  is  to  study  mechanisms  through  which  pressure  waves  may  be  generated  by 
microwave  pulses,  and  to  derive  a  relationship  between  microwave  pulses  and  the 
generated  pressure  waves. 

3.3.  Microwave  Scattering,  Inverse  Scattering,  Dosimetry,  and  Imagery 

This  task  is  to  develop  non-invasive  techniques  for  dosimetry  of  dielectric  bodies  under 
microwave  exposure,  and  to  develop  algorithm  for  three-dimensional  medical  imaging  using 
low-level  microwaves. 


4.  PROJECT  PROGRESS 

During  the  first  year  of  this  project,  research  was  performed  on  developing  a  non-invasive 
microwave  dosimetry  technique,  and  on  continuing  previous  theoretical  studies  on  transient 
dielectric  interaction  and  microwave  to  acoustic  energy  conversion.  As  progress  was  being 
made  in  WRAIR  on  the  installation  of  a  high-power  pulsed  microwave  generator,  our  effort  in 
the  second  year  was  focused  on  developing  theoretical  basis  to  support  forthcoming  experiments 
in  WRAIR  on  biological  effects  of  transient  and  high  power  microwave  damage.  On  the 
transient  effects,  we  have  developed  a  theoretical  basis  for  experimental  observation  of 
temporal  non-linearity  in  the  transient  regime,  which  may  be  measurable  and  applicable  to  the 
experimental  setup  being  planned  at  >VRAIR.  On  microwave  acoustics,  we  incorporate 
electrostrictive  effect  in  the  pervious  formulation,  and  established  criteria  to  predict  the 
dominant  mechanism,  thermoacoustic  or  electrostrictive,  for  generating  pressure  waves  by 
microwave  pulses.^-^  On  dosimetry,  we  have  a  major  breakthrough  on  developing  an  algorithm 
for  evaluating  the  vector  fields  both  inside  and  outside  a  three-dimensional  body  of  arbitrary 
geometry  and  dielectric  profile.-’  We  have  also  developed  a  formula  for  computing  the  vector 
fields  inside  the  dielectric  body  from  measurement  of  scattered  fields  in  a  limited  region 
outside  the  body.  Our  success  in  developing  the  non-invasive  dosimetry  algorithm  also  led  us 
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to  discover  a  non-diffractional  imagery  technique,  which  may  provide  images  with  resolution 
limited  only  by  signal-to-noise  ratio,  and  not  by  wavelength  of  the  probing  microwave  or  by 
geometrical  configuration  of  the  receiving  antenna.'^ 

The  effort  of  this  program  during  the  two  funded  years  have  resulted  in  publication  of  six 
papers.^-^'^’^'^'^  The  following  gives  a  technical  summary  of  the  progress  and  status  in  each 
task  of  this  program.  Details  of  research  results  are  given  in  the  appendices,  which  are 
reprints  of  our  publications. 


4.1.  Transient  Dielectric  Response 
Theory  of  Transient  Response 

Our  effort  on  this  subject  has  been  focused  on  developing  theoretical  basis  to  support 
forthcoming  experiments  in  WRAIR  on  nonlinear  biological  effects  of  transient  and  high  power 
microwave  damage.  The  objective  was  to  devise  a  theoretical  scheme  for  experiments  that  are 
able  to  isolate  the  effects  of  short  and  high-intensity  microwave  pulses  from  low-intensity 
continuous  waves  (CW)  on  biological  materials.  One  may  directly  analyze  the  materials  and 
compare  the  damage  between  pulse-wave  exposure  and  CW  exposure.  However,  from  the 
observation  of  the  damage  alone,  it  is  not  possible  to  conclude  if  the  observed  pulse  or 
high-intensity  effects  are  direct  primary  electromagnetic  effects  or  indirect  non-electromagnetic 
secondary  effects  that  are  induced  by  the  microwave  pulses.  Therefore,  instead  of  directly 
analyzing  damage  to  biological  materials,  we  studied  the  approach  of  analyzing  the  pulse 
propagation  and  the  field  strength  of  the  microwave  pulse. 

Our  previous  theoretical  studies  have  shown  that,  due  to  a  non-steady  state  and 
irreversible  process  in  the  transient  regime,  dielectric  relaxation  cannot  be  parameterized  by 
time-independent  parameters,  such  as  the  damping  coefficient  in  the  transient  regime. 
However,  for  experimental  and  practical  purpose,  macroscopic  parameterization  is  desirable. 
Therefore  we  propose  the  following  time-dependent  model  in  the  transient  regime  for  the 
dielectric  polarization  current,  J(t),  which  is  the  time-derivative  of  the  dielectric  polarization 

P(tywjo 


x(t)-^  + J(t) 

- 

-  S'dT 

(1) 

(1 

-  exp(-t/T;)] 

(2) 

The  model  was  conceived  from  the  reasoning  that  the  damping  coefficient  vanished  initially 
and,  as  molecules  approached  an  equilibrium  or  a  steady  state,  it  gradually  increased  to  a  final 
value,  viz.,  Debye’s  steady-state  value.  The  time  it  takes  to  reach  the  steady-state  value  is 
characterized  by  xr,  the  function  H(t)  is  the  Heavyside  function  which  is  zero  for  t<0,  1/2  for 
t=0,  and  1  for  t>0.  With  this  model,  it  was  found  that  the  dielectric  is  nonlinear  in  the 
transient  regime  in  a  way  that  any  frequency  component  of  the  the  dielectric  response  may 
depend  on  the  action  of  different  frequencies.  Thus,  given  an  electric  field  as  an  action,  then 
the  dielectric  polarization  in  the  frequency  regime,  P(a)),  depends  on  E(q))  in  the  following  way: 


P(w)  =  J  xC^.wOEfto'j-dto' 


(3) 


where  x(cd,(jo')  is  the  dielectric  susceptibility.  In  mathematical  terms,  dielectric  susceptibility  is 
not  diagonal  in  the  frequency  regime.  This  means  that,  in  the  time  regime,  response  at  a 
certain  time,  t',  caused  by  an  action  at  an  earlier  time  t,  does  not  depend  on  the  difference 
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t-t'  only.  Therefore,  linear  superposition  of  action-response  relationship  is  no  longer  valid. 
This  is  expected  because,  in  the  transient  regime,  the  dielectric  is  in  a  non-steady  state  and 
undergoes  irreversible  processes.  We  also  derived  the  function  x(w,0)')  from  eqs.  1  and  2,  and 
obtained:^-® 

^  (O'  r(-ioni+Xj/Xoyr(-i(si'x,+l)  y  1  _ 

^  2tU\o  to+iO' r(-zaycj+l)-r(-ia)'x7-K7/To)  iio  (w'-oyinjxdi-iai'xi+n+xi/xo) 


g  Xi  _  CO'  r(-/{l)T/-K//To)-r(-/CO'Tj-Hl) 

2^'  Xo'  ((j>i-jO)(o)'-ci>-iO)  T(-ja)t/+l)-r(-Wt;-K7/to+l) 

X  jF2(i'(o)'-fo)X;, -iw'+tt/to,  1; /(a)'-fO)X7+l,-/a)'X7+X7/Xo-i-l;  1) 


where  T  is  the  gamma  function  and  mF,  is  the  generalized  hypergeometric  function^^ 

There  are  several  experimental  indications  that  electromagnetic  pulses  with  pulse  width  of 
the  order  of  nanoseconds  or  shorter  may  produce  transient  effects  in  biological  materials.  One 
of  these  experimental  indications  concerns  the  life-times  of  excited  vibrational  modes  of  DNA 
chains.  Based  on  the  line-widths  of  vibrational  excitations  by  microwaves  observed  recently, 
the  life-times  are  of  the  order  of  120  nanoseconds.'^-'^  The  time  it  takes  to  reach  an 
equilibrium  state  or  steady  state  must  be  longer  than  these  life-times.  Another  experimental 
evidence  of  transient  phenomenon  is  the  recent  observation  of  time-dependency  of  dielectric 
su.sceptibility  of  an  inorganic  solid  under  high  field.'^  The  observation  was  derived  from 
applying  electric  fields  in  the  range  of  megavolts  per  centimeter  to  thin-film  of  amorphous 
aluminum  oxide,  AhOj.  It  was  found  that  tiie  dielectric  susceptibility  depends  on  time  for  as 
long  as  hundred  seconds.  While  the  value  of  Xo,  which  is  the  steady-state  Debye’s  relaxation 
time,  is  known  for  most  dielectrics,  the  exact  values  of  Xi  have  never  been  measured.  The 
above  experimental  observations  indicate  that  non-steady  state  transient  regime  exists  in 
dielectrics  and  it  may  last  longer  than  the  microwave  pulse  width.  Availability  of  the  values  of 
Xi  will  provide  much  knowledge  on  dielectric  response  to  high-intensity  shon  microwave  pulses. 
Therefore,  our  research  was  focused  on  devising  some  theoretical  base  for  an  experimental 
measurement  of  Xy  for  some  materials,  especially  for  solution  of  biological  materials  in  liquid 
water.  For  dielectrics  composed  of  large  molecules,  it  is  expected  that  X;  will  be  quite  long 
comparing  to  pulse  width  of  existing  microwave  pulse  systems.  From  the  observation  of  the 
excited  vibrational  modes  of  DNA  chains,'^-'^  it  indicates  that  Xi  must  be  of  the  order  of  120 
nanoseconds  or  longer.  This  suggests  that  transient  non-steady  state  effect  may  play  an 
important  role  in  response  of  DNA  solutions  to  nanosecond  microwave  pulses. 


Experimental  Proposal 

The  above  result  suggests  that  any  transient  effect  or  high  field  effect  in  dielectrics  is 
expected  to  produce  non-linear  harmonic  and  anharmonic  generation  in  the  electromagnetic 
spectrum.  Conventional  non-linear  analyses  employ  network  analyzers  to  study  the  output 
spectrum  from  dielectric  specimens  which  are  exposed  to  a  monochromatic  wave;  non-linearity 
is  inferred  if  harmonic  generations  are  observed  in  the  output  spectrum.  This  method  cannot 
be  applied  to  pulsed  fields  since  the  input  contains  a  wide  spectrum  of  frequencies.  Therefore 
different  method  must  be  undertaken  to  analyze  transient  effects  and  high  field  pulse  effects. 
To  this  end,  we  have  found  a  macroscopic  effect  of  the  transient  state  model  that  is  applicable 
to  pulse  inputs;  the  effect  may  also  be  measurable  and  applicable  to  an  experimental  setup  that 
has  been  planned  at  WRAIR.  To  briefly  summarize,  our  approach  is  to  analyze  the  scattering 
matrix  of  a  dielectric  specimen  under  an  incident  pulsewave.  The  scattering  matrix  is 
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constructed  from  the  measurements  of  input  and  output  waves  and  analyzed  in  the  frequency 
domain.  Transient  effects  or  high-field  pulse  effects  are  inferred  if  the  scattering  matrix  have 
off-diagonal  elements.  This  scattering  matrix  approach  appears  to  be,  in  a  very  broad  sense, 
similar  to  the  inverse  scattering  approach  that  we  have  developed  in  microwave  dosimetry  and 
imagery.  Indeed,  all  physical  techniques  which  infer  the  characteristics  of  unknown  objects  by 
analyzing  their  action-response  data  are,  in  a  broad  sense,  inverse  scattering  and  imagery.  The 
difference  is  in  the  physical  characteristics  that  are  being  "imaged".  In  the  scattering  matrix 
approach  described  here,  it  is  the  temporal  and  spectral  dielectric  response  that  is  "imaged", 
whereas  in  conventional  imaging,  the  spatial  distribution  of  some  physical  characteristics  is 
imaged. 


Microwave  -> 


D 

Distance  of  microwave  path 


Figure  1.  Sketch  of  experimental  setup  to  measure  the  field  strength  along  the  path 
of  a  microwave  pulse  in  a  solution  of  a  biological  material  in  water. 


To  describe  the  scattering  matrix  approach  in  more  detail,  let  us  consider  a 
one-dimensional  propagation  in  an  experimental  setup  as  illustrated  in  Figure  1.  Take  two 
points  along  the  propagation  path  of  a  microwave  pulse,  say  xi  and  xz.  The  material  between 
these  two  points  may  be  regarded  as  the  specimen;  the  electric  field  of  the  microwave  pulse  at 
X;  may  be  considered  as  an  input  field  to  this  specimen,  and  that  at  xz  may  be  considered  as 
the  output  field  from  the  specimen.  Both  the  input  field,  denoted  by  and  output  field, 
denoted  by  are  functions  of  time.  By  making  time-resolved  measurements  of  these  two 
functions  at  a  certain  sampling  rate,  one  obtains  two  column  vectors  in  the  time  domain;  these 
two  vectors  may  be  Fourier- transformed  to  the  frequency  space.  We  shall  construct  a 
scattering  matrix  from  the  measurements  of  a  set  of  input  fields  and  the  corresponding  set  of 
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output  fields.  If  there  are  total  of  N  samples  in  the  time-resolved  measurements,  then  these 
two  vectors  are  of  N-dimensional.  In  order  to  construct  the  scattering  matrix,  one  must  then 
prepare  N  input  fields  and  N  output  fields.  So,  measurements  of  the  fields  of  N  pulses  must  be 
made.  Let  and  where  n  =  1,  2,  ...,  N,  be  the  sets  of  output  vectors  and  input 

vectors,  respectively.  From  these  two  sets  of  fields,  we  may  then  construct  the  input  matrix 
and  output  matrix  as  follows:  Lining  up  the  N  input  vectors,  one  by  one  as  columns  of  a 
matrix,  one  obtains  a  square  matrix  which  is  the  input  matrix;  similarly,  lining  up  the  N 
output  vectors,  one  by  one  as  columns  of  a  matrix,  one  obtains  the  output  matrix, 

It  can  then  be  proven  that  the  scattering  matrix  is  the  product  of  these  two  matrices; 

S  =  FW.[FW]'.  (4) 

The  scattering  matrix  obtained  this  way  has  the  property  that 

=  (n  =  l,...,N).  (5) 

It  is  remarked  that  the  N  microwave  pulses  must  be  linearly  independent,  viz.,  the  input  matrix, 
F^‘^,  must  be  non-singular.  In  other  words,  the  determinant  of  F'^'^  must  not  vanish.  Instead 
of  making  measurements  on  N  input  pulses,  one  may  alternatively  use  one  single  pulse,  and 
make  time  resolved  measurements  at  N+l  successive  points  of  equal  spacing  along  the  path  of 
the  pulse.  Thus,  the  pulse  measured  at,  say  x„  is  the  output  pulse  of  the  previous  layer  of 
material  and  the  input  pulse  of  the  next  layer  of  material.  Since  the  dielectric  is  assumed  to 
be  homogeneous,  the  layers  of  material  between  all  successive  pairs  of  points  are  identical.  So, 
the  pulses  at  x;,  xj,  ...,  xa^  are  N  input  pulses  to  the  specimen,  and  the  pulses  at  X2,  xj,  ..., 
xw+/  are  corresponding  output  pulses.  Owing  to  dielectric  dispersion  and,  perhaps,  also 
transient  or  other  non-linear  effects,  these  N  input  pulses  will  be  linearly  independent,  and 
may  then  be  used  to  construct  the  input  matrix. 

Based  on  the  transient  state  model  that  was  discussed  earlier  (cf.  eqs.  1-4),  we  have 
studied  the  properties  of  the  scattering  matrix  S  constructed  from  the  field  measurement  data 
as  described  in  the  last  paragraph.  We  found  a  distinct  property  which  may  be  used  to  isolate 
transient  effect  from  continuous  wave  (CW)  effect.  If  there  is  any  transient  response  in  the 
biological  material,  or  any  dielectric  material,  between  the  two  points  x/  and  X2,  then  the 
matrix  S  will  not  be  diagonal  in  the  frequency  domain.  Conversely,  if  there  is  no  transient 
effect,  then  the  scattering  matrix  must  be  diagonal  in  the  frequency  domain.  Thus,  this 
property  may  be  used  to  isolate  transient  response  from  CW  effect.  This  approach  is  applicable 
to  the  measurement  setup  planned  at  WRAIR.  While  we  believe  that  the  basic  principle  of  this 
approach  is  sound,  there  is  still  a  technical  problem  that  needs  further  investigation.  The 
problem  concerns  the  finite  sampling  rate  and  finite  time  duration  of  measurement.  Ideally, 
theory  requires  that  infinitely  many  pulses  be  measured  at  both  x/  and  X2,  and  that 
measurements  be  made  continuously  in  time.  Finite  sampling  of  measurement  on  finite  number 
of  pulses  will  then  cause  some  "spreading"  of  the  diagonal  elements  of  the  matrix  S  into 
off-diagonal  part,  thereby  produce  some  mixing  of  the  CW  effect  with  transient  effect.  The 
effect  of  this  spreading  may  be  treated  as  a  noise  in  the  measurement  data.  The  magnitude  of 
this  noise  level  remains  to  be  investigated. 

The  effect  of  off-diagonal  matrix  elements  in  transient  response  is  equivalent  to  harmonic 
generation  in  non-linear  circuit  theory.  Therefore,  the  principle  of  the  above  approach  also 
applies  to  high-field  effect,  however,  there  is  a  problem  that  needs  separate  consideration. 
Any  high-Field  effect  must  give  rise  to  a  dependency  in  the  output  fields,  therefore,  some 
normalization  method  may  dlso  be  needed  on  the  output  fields,  as  well  as  the  input  fields.  The 
advantage  of  this  method  over  conventional  method  of  analyzing  harmonic  generation  is  that  it 
is  not  limited  to  monochromatic  input  source;  pulse  input  of  arbitrary  waveform  may  be 
employed. 
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4.2.  Stress  Waves  and  Pressure  Waves  Induced  by  Microwave  Pulses 

Acoustic  waves  generated  by  pulsed  microwaves  have  been  cited  as  a  mechanism  for 
microwave  hearing.^^-"®-^^  More  recently  it  has  been  demonstrated  that  acoustic  waves  are 
transduced  in  dielectric  objects  simulating  the  ocular  lens  when  exposed  to  pulsed 
microwaves.'®  The  effect  has  also  been  cited  as  the  operant  mechanism  for  cellular  damage  in 
studies  of  the  murine  ocular  lens  in  vitro.'^  Three  theories  on  the  microwave  auditory 
mechanism  have  been  reviewed  and  compared  by  Lin,^"  viz.,  the  radiation  pressure  theory,  the 
electrostrictive  theory,  and  the  thermoacoustic  theory.  Among  them  it  was  found  that  only  the 
thermoacoustic  mechanism  could  produce  elastic  waves  of  magnitude  large  enough  to  explain  the 
experimental  observations.  Comparing  to  the  magnitude  of  acoustic  pressure  generated  in 
typical  biological  tissues,  the  electrostrictive  effect  was  found  to  be  about  two  orders  of 
magnitude  smaller,  while  the  radiation  pressure  three  orders  of  magnitude  smaller.^' 

Previous  theories  on  microwave  electrostrictive  process  and  microwave  thermoacoustic 
process  are,  however,  based  on  a  linear  wave  equation  with  a  generating  function  derived  from 
the  microwave  pulses.  Pressure  waves  are  thereby  generated  in  a  way  similar  to  a  forced 
harmonic  oscillator.  In  the  case  of  the  thermoacoustic  process,  the  generating  function  is 
derived  from  an  inhomogeneous  heating  by  the  microwaves,  whereas  in  the  case  of  the 
electrostrictive  process,  it  is  derived  from  stress  due  to  dielectric  polarization.'',^®  This 
approach  fails  to  consider  the  proper  balance  of  the  distribution  of  the  absorbed  microwave 
energy  among  the  internal  thermd  energies  and  the  bulk  kinetic  energy,  and  the  effect  of 
thermal  or  dielectric  discontinuities  at  medium  interface.  Also  neglected  is  the  coupling 
between  the  thermoacoustic  process  and  electrostrictive  process,  which  may  result  from 
dependency  of  dielectric  permittivity  on  various  thermodynamical  coefficients,  such  as  mass 
density,  temperature,  and  pressure;  the  transient  effect  of  ultra-short  electromagnetic  pulses 
may  further  enhance  the  electrostrictive  effect.  We  approached  these  problems  by  maWng  a 
thorough  formulation  of  the  coupling  of  microwave  pulses  to  pressure  waves,  with  particular 
emphasis  on  material  discontinuity  and  electromagnetic  transient  effect.  Elastic  wave  equations 
were  then  derived  for  both  thermoacoustic  and  electrostrictive  effects,  and  some  models  of  air- 
water  system  were  used  to  make  numerical  computation  and  estimate  the  generated  pressure 
waves  in  water.  Contrary  to  previous  theory  that  thermoacoustic  pressuie  waves  were 
generated  mainly  by  inhomogeneous  distribution  of  microwaves  in  dielectrics,'''  it  was  found 
that  pressure  waves  are  generated  whenever  there  i'  a  discontinuity  in  thermal  or  dielectric 
parameters  in  the  medium.  As  to  the  ratio  of  electrostrictive  effect  to  thermoacoustic  effect, 
it  depends  on  the  rise  time  of  the  microwave  pulse,  being  greater  than  one  for  rise  time 
shorter  tha.;  a  nanosecond,  and  smaller  than  one  for  rise  time  longer  than  a  nanosecond.  This 
result  is  in  contrary  to  previous  estimate  based  on  a  simpler  theory,  which  gave  a  ratio  of  the 
order  of  10 ',  and  independent  of  the  rise  time  of  microwave  pulses. 

The  following  paragraphs  provide  a  summary  of  our  results  in  thermoacoustic  and 
electrostrictive  effects.  Details  of  the  formulation  are  given  in  Appendices  A,  B,  and  F  and 
other  references  cited  therein. 


Thermoacoustic  Effect 

Our  approach  to  derive  the  coupling  from  microwave  pulses  to  acoustic  waves  started  from 
the  basic  laws  of  thermodynamics:  the  conservation  of  mass,  the  conservation  of  momentum,  the 
conservation  of  energy,  and  the  thermodynamic  equation  of  state;the  thermodynamic  equation  of 
state  was  represented  by  the  following  expression  of  internal  energy  in  terms  of  some 
thermodynamical  coefficients: 

d,U/  =  [-Cp/(pPp)  +  (l/p')p.>0;VyV-v]d,p  +  (Cv/Pp)(Kr)vd,p,/  (7) 
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which  we  derived  from  the  first  and  second  laws  of  thermodynamics.^'*  In  the  above  equation, 
p  is  the  mass  density,  v  is  the  bulk  velocity,  py  is  the  stress  tensor,  Cv  and  Cp  are  the 
specific  heats  (per  unit  mass)  at  constant  volume  and  constant  pressure,  respectively,  pp  is  the 
isobaric  thermal  expansion  coefficient,  and  (Kr)i>  is  the  isothermal  compressibility  tensor.  By 
mutual  substitution  among  equations  representing  these  laws,  one  may  obtain  the  wave 
equations  for  any  of  the  thermodynamic  quantities.  The  boundary  conditions  may  be  obtained 
by  first  transforming  these  equations  to  the  Lagrangian  specification  (in  which  the  boundary 
surface  is  stationary)  and  then  integrating  each  term  across  a  thin  layer  of  the  boundary 
interface.'*-^'*  If  the  strain  tensor  is  isotropic,  one  obtains  the  following  wave  equation  for  the 
stress  tensor,  p,/; 


d,[(pJpr^iCJCp)iKr).^:p.j]  -  a,[(p./p)"^(l/p)^py] 

-  a.:(p./P)"^(P./Cp)P3  (8) 

where  Cv  and  Cp  are  the  specific  heats  (per  mass)  at  constant  volume  and  constant  pressure, 
respectively,  pp  the  isobaric  thermal  expansion  coefficient,  Kr  the  isothermal  compressibility,  p 
the  mass  density,  P  the  microwave  specific  absorption  rate  (SAR),  and  the  subscript  o  signifies 
that  the  quantity  is  evaluated  at  its  ambient  equilibrium  value. 

To  illustrate  the  implications  of  the  wave  equations  and  their  respective  boundary 
conditions,  we  consider  a  dielectric  sphere  of  radius  a  surrounded  by  air  at  1  atmospheric 
pressure.  A  square  pulse  of  microwave  of  duration  x  and  amplitude  Po  is  incident  upon  the 
sphere.  The  sphere  is  assumed  to  be  small  so  that  the  microwave  absorpdon  exhibits  no  spatial 
variation.  In  the  linear  approximation,  the  solutions  of  pressure  waves  in  the  frequency  domain 
are  then: 

p/(r,(u) 


and 

P2(r,w) 


where  j„(x)  and  h„(x)  are,  respectively,  the  spherical  Bessel  and  Hankel  functions  of  the  first 
kind  of  order  n,  the  subscripts  1  and  2  label  the  dielectric  and  the  air  media,  respectively, 
Cij  -  [Cp/pCvK7]*'^  tancp  =  (pc)2/(pc)/,  and  P(o))  is  the  Fourier  transform  of  the  square  pulse 
SAR,  which  is  given  by: 


=  -i(^) 


P((d) 

‘  o)-i0 


■j^(k;r) 
lio(k2a)' 


-  I 


jp(k/a)  -  tantp  •  •  j/(k;a) 


(9) 


c^ppp,  P((0)  h„(k2r) 


yk/Q) 


j<,(k2fl)  -  tantp  •  •  j/Ckifl) 


-  1 


(10) 


=  (>» 

Similar  results  may  also  be  obtained  for  the  bulk  velocity  at  either  side  of  the  interface.  The 
total  acoustic  energy  coupled  into  the  air  may  be  obtained  by  calculating  the  work  done  by  the 
dielectric  on  the  air  at  the  interface.  To  the  first  order,  the  result  is: 

Ea,r  =  Ana^po  ]  v/(a,t)dt  =  Ana^poiln)"^  ■  V2(a,co=0) 


=  (27r)*'2  ^  p„[-^],  ■  P(co=0) 


(12) 
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Since  Eati  =  (27c)^'^(4Ka73)prP(a>=0)  is  the  total  microwave  energy  absorbed  through  the  entire 
pulse,  coupling  efficiency  of  the  absorbed  microwave  energy  to  the  acoustic  energy  in  the  air  is 
then: 


r[air 


Eaijj 


(13) 


Take  the  dielectric  to  be  water  at  30  °C,  so  =  2.8x10"'  '  and  Cp  =  4.186x10^  erg/gm-°C. 
With  the  equilibrium  pressure  of  the  air,  po,  being  1  atm  =  l.OlxKP  dynes/cm^  the  coupling 
efficiency  is  then  Tio.>  =  6.7X10  '’  for  a  sphere  of  1  cm  diameter. 


The  minimum  pressure  amplitude  generated  in  the  water  may  be  estimated  by  the 
fundamental  harmonic  in  eq.  9.  Taking  an  optimal  pulse  width  of  x  =  Jt/cOo,  it  gives  a  pressure 
amplitude  of  (2a/7T)(cp^/Cj;Pa.  Using  the  values  of  p;,  and  Cp  for  30  °C  as  cited  above,  a 
peak  SAR  of  K  =  15  kW/gm  results  in  a  pressure  wave  of  amplitude  equal  to  0.1  bar  in  the 
water, which  is  10%  of  the  initial  ec'iilibrium  pressure.  This  fraction  of  variation  in  pressure 
will  also  result  in  at  least  equal  fraction  of  error  in  the  linear  approximation.  Similar 
calculation  was  also  made  for  a  one-dimensional  air-water  system.^^  It  was  found  that, 
accounting  for  only  the  n  =  ±1  term,  a  peak  SAR  of  6.5  kW/gm  results  in  a  pressure  wave  of 
amplitude  equal  to  0.1  bar  in  the  water;’*’  if  up  to  the  first  11  terms  are  included,  a  peak  SAR 
of  4.5  kW/gm  suffices  to  generate  pressure  waves  of  such  amplitude. 


Electrostrictive  Effect 

We  approached  the  problem  of  electrostrictive  effect  by  including,  in  the  equations  of 
conservation  of  momentum  and  conservation  of  energy,  the  following  electrostrictive  tensor:^’ 

The  total  stress  on  the  dielectric  is  the  sum.  of  material  stress  and  electrostrictive  stress: 


SiJ  —  p,y  +  0,J  (15) 

However,  it  is  the  material  stress  p,/  which  is  responsible  for  the  material  strain,  and 
therefore  accountable  in  the  thermd  internal  energy.  Indeed,  in  the  presence  of  an 
electromagnetic  field,  the  material  stress  automatically  adjusts  itself  in  reaction  to  the 
electrostrictive  stress.  Therefore,  in  deriving  the  equation  of  state  which  is  represented  by  an 
expression  of  internal  ene'gy  in  terms  of  some  thermodynamical  coefficients  (cf  eq.  7),  only 
the  material  stress  is  included.  Upon  some  isotropic  assumption,  we  then  derive  the  following 
wave  equation  for  the  total  stress  tensor:'-^ 

+  ^(Kr).ya,s.y) 


p.9<P  -t- 


(16) 


where  the  subscript  o  signifies  that  the  quantity  is  evaluated  at  iis  ambient  equilibrium  value, 
and  o  is  the  isotropic  electrostrictive  pressure: 
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For  an  isotropic  dielectric,  s.y  =  -s5.y  and  (Kr),y  =  Kt/3,  so  ]L(Kr)u  (sum  over  i)  is  Kr,  the 
above  wave  equation  is  then  further  reduced  to 


Equation  18  is  a  typical  linear  wave  equation  with  the  velocity  c: 


where  u  denotes  the  acoustic  velocity  in  the  dielectrics  in  the  absence  of  the  electromagnetic 
waves.  The  electromagnetic  waves  change  the  acoustic  velocity  through  the  electrostrictive 
pressure  a.  To  estimate  the  magnitude  of  this  correction,  we  consider  water  at  room 
temperature  as  the  dielectric,  for  which  =  2.8x10''  "C  ‘,  Cp  =  4.186x10^  erg/gm-*C,  p  = 
1  gm/cm^,  and  e  =  80.  We  also  assume  that  the  dielectric  permittivity  is  linearly  proportional 
to  the  mass  density,  so  (0e/3p)r  =  e/p,  then  eq.  17  gives  a  =  (elEP)/(247r).  Assume  a  plane 
electromagnetic  wave  of  intensity  I  =  1  kW/cm^  (which  is  equivalent  to  an  electric  field  of 
amplitude  2.9x10'  volts/m),  then  its  field  magnitude  is  lEP  =  8nl/(cyje)  =  0.94  erg/cm^,  where  c 
is  the  speed  of  light  in  vacuum.  Then  a  =  1  dyne/cm^  and  the  fractional  change  to  the 
acoustic  speed  is  (oPp)/(pCp)  =  6.7x10  which  is  insignificant.  While  the  fractional  change  of 
acoustic  speed  is  negligibly  small  for  most  dielectrics  under  currently  available  microwave 
sources,  it  may  be  appreciable  for  some  dielectrics  under  lower  frequency  electromagnetic 
fields. 

As  illustrated  at  the  end  of  the  last  paragraph,  for  most  dielectrics  and  for  microwave 
intensity  up  to  gigawatts/cm^  the  term  o/p  on  both  sides  of  eq.  18  may  be  neglected.  Eq.  18 
may  then  be  rewritten  as 


where  u  =  [pC.Kr/Cp]'''^  is  the  speed  of  sound  in  the  dielectric  in  the  absence  of  the 
electrostrictive  force.  The  three  terms  inside  the  brackets  on  the  right  hand  side  are  the 
generating  forces  of  pressure  waves.  The  first  term  represents  the  thermoelastic  effect,  the 
second  term  is  the  electrostrictive  force,  and  the  third  term  comes  from  the  electromagnetic 
energy  of  dielectric  polarization.  The  above  equation  will  be  our  basis  for  comparing  the 
relative  strength  of  the  three  forces  on  microwave  to  pressure  wave  coupling.  From  the 
expression  ot  these  three  forces  in  the  above  equation,  one  sees  that  time  variation,  namely, 
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transient  effect,  of  the  microwave  pulses  has  an  important  effect  on  the  electrostrictive  force 
and  on  the  polarization  energy.  Here  we  may  make  a  preliminary  estimation  of  the  relative 
sizes  of  these  three  coupling  forces.  First  we  express  them  in  terms  of  lEP.  Let  a  be  the 
attenuation  coefficient  of  microwave  intensity  as  defined  by  the  equation  I(x)  =  Le  ".  Then 
the  first  term  is  poP  =  al  =  (acVe/87c)IEP.  If  the  dielectric  permittivity  is  linear  in  p,  then 
(3e/9p)r  =  e/p,  so  eq.  17  gives  a  =  (elEP)/(24jt)  and  the  electrostrictive  term  gives 
[(peCyKr)/(247tPp)]9,lEP.  Denote  by  u  the  rise  time  of  a  microwave  pulse,  then  3,IEP  ~ 
lEP/b.  Thus,  the  relative  magnitudes  of  the  electrostrictive  and  polarization  terms  with 
respect  to  the  thermoelastic  term  are,  respectively. 


Electrostrictive  effect 
Thermoelastic  effect 


pCvK^e  1. 
Scajip  u 


Polarization  effect  _  -y/e  1 

Thermoelastic  effect  ~  ca  t. 


(21) 


One  sees  that  the  second  and  the  third  driving  forces  may  be  larger  than  the  thermoelastic 
term  if  the  microwave  is  rapidly  modulated,  contrary  to  the  estimate  by  previous  theory  that 
the  electrostrictive  effect  is  two  orders  of  magnitude  smaller^'’.  For  water  at  30  ‘C,  p  = 
1  gm/cm^,  Q  =  1.013  Cv  =  4.186x10^  erg/gm“C,  Kr  =  4.46x10"'^  cmVdyne,  =  2.8x10"*  °C  L 
and,  at  3  GHz,  a  =  0.88  cm  *  and  e  =  80.  The  above  ratios  are  then,  respectively: 


Electrostrictive  effect 
Thermoelastic  effect 


0.75  nanosec 

tr 


Polarization  effect 
Thermoelastic  effect 


0.34  nanosec 

tr 


(22) 


One  may  conclude  that,  for  water,  a  rise  time  of  I  nanosecond  is  the  marginal  modulation  rate, 
beyond  which  the  electrostrictive  term  and  the  polarization  term  contribute  more  to  the 
pressure  waves  than  the  thermoelastic  term. 


Comparison  with  Experimental  Measurements 

Pressure  measurements  were  made  at  WRAIR  in  air-water  system  similar  to  the  models  in 
our  theoretical  computation.  Pressure  waves  were  measured  in  both  air  and  saline  liquid  for 
two  types  of  vessels  placed  in  an  WR  975  exposure  system.  The  first  vessel  was  ca.  one  inch 
on  each  side  with  a  1/8  inch  wall  thickness.  The  second  vessel  was  circular  in  cross  section, 
one  inch  in  outer  diameter,  1/8  inch  wall  thickness,  and  made  of  acrylic  plastic.  Both  were 
filled  with  phosphate  buffed  saline.  These  vessels  were  placed  into  the  WR  975  via  the  same 
waveguide  below  cut-off  window  as  the  actual  lenticular  exposure  chamber.  The  two  vessels 
extended  beyond  the  waveguide  to  provide  a  liquid  column  for  coupling  of  the  acoustic  wave 
into  regions  where  the  hydrophone  would  not  be  subject  to  direct  influence  of  the 
electromagnetic  field.  Each  vessel  was  separately  matched  to  the  microwave  source  with  the 
triple  stub  tuner.  Typical  return  losses  were  ca.  -15  to  -20  dB  for  either  vessel.  All  liquid 
based  pressure  measurements  are  expressed  in  dB  relative  to  0  dB  =  1  micro  Pascal,  and  all  air 
based  pressure  measurements  are  made  in  dB  relative  to  0  dB  =  20  micro  Pascal.  Both  air  and 
water  pressure  measurements  were  tested  for  direct  electromagnetic  field  effects.  The  results 
of  the  measurements  appeared  to  be  in  agreement  with  our  theoretical  calculation.** 
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4.3  Microwave  Scattering,  Inverse  Scattering,  and  Dosimetry 

An  important  practical  aspect  of  microwave  biological  damage  concerns  the  dosimetry  of 
microwaves  in  biological  bodies.  Biological  bodies  are  complicate  dielectric  structures  and 
multiple  scattering  and  reflection  from  dielectric  interfaces  often  make  it  impossible  to  predict 
microwave  dose  distribution  in  the  bodies.  Experimental  measurements  using  insertion  technique 
are  used,  but  the  instruments  often  produce  unpredictable  disturbances  on  the  electromagnetic 
field  and  thus  render  the  measurements  inaccurate.  Therefore  non-invasive  measurement 
techniques  are  much  desired.  Our  effort  in  this  subject  has  been  to  develop  an  inverse 
scattering  approach  to  accomplish  microwave  dosimetry  by  measuring  scattered  fields. 

There  are  two  aspects  in  this  subject.  One  is  the  determination  of  scattered  fields  from 
the  dielectric  profile  of  a  dielectric  object,  and  another  concerns  the  determination  of  the 
electric  fields  inside  the  target  from  measurement  of  the  scattered  fields  in  a  limited  region 
outside  the  target.  Both  aspects  have  been  under  intensive  investigation  by  many  researchers. 
On  the  determination  of  scattered  fields  from  dielectric  profiles  of  scattering  objects,  previous 
approaches  by  other  investigators  included;  1)  Bom  approximation,  of  which  the  accuracy  was 
limited,^*  2)  the  moment  method  by  Richmond,  which  applied  to  only  two-dimensional  objects, 
viz.,  objects  with  cylindrical  symmetry,^®-’®  and  3)  linear  parameter  technique  with  moment 
method,  which  approximated  the  field  inside  the  scatterer  by  a  linear  combination  of  some  basis 
functions  and  then  numerically  solved  the  linear  coefficients.  Therefore  our  main  effort  on 
this  aspect  has  been  to  extend  Richmond’s  moment  method  to  three-dimensional  object  with 
arbitrary  geometry  and  dielectric  profile.  Our  effort  has  been  quite  successful.  We  have 
developed  the  following  algorithm  to  compute  the  scattered  field  of  arbitraiy  three-dimensional 
dielectric  ubjects  from  their  dielectric  profile.  As  to  the  second  problem,  viz.,  determination  of 
electric  fields  inside  a  target  from  measurement  of  scattered  fields  in  a  limited  region  outside 
the  target,  it  is  the  problem  of  inverse  scattering,  and  is  also  related  to  microwave  imaging; 
once  the  field  inside  the  target  is  known,  one  may  then  obtain  a  dielectric  image  in  terms  of 
the  profile  of  dielectric  permittivity  of  the  target.  On  this  subject,  we  also  developed  two 
inverse  scattering  approaches.  The  first  approach  utilizes  a  method  of  "soft  focusing"  of 
scattering  data  to  reconstruct  qualitative  dielectric  images  with  negligible  computing  time.  The 
second  approach  utilizes  a  method  of  converting  the  integral  wave  equation  to  matrix  equation 
and  reconstructs  images  by  matrix  inversion.  This  metitod  provides  images  with  resolution 
limited  only  by  signal-to-noise  ratio  and  computer  power,  and  not  by  wavelength.  Both 
algorithms  achieve  three-dimensional  imaging  by  measuring  scattered  waves  in  the  near  zone 
where  the  wave  structure  is  correlated  to  the  "depth"  of  the  target. 


Microwave  Scattering 

The  objective  of  this  problem  was  to  develop  a  method  to  compute  scattered  microwave 
fields  from  an  arbitrary  three-dimensional  inhomogeneous  object.  Our  approach  is  to  convert 
the  integral  equation  into  a  matrix  equation  by  digitizing  the  target  space  and  the  measurement 
space.  First,  denote  by  Vo  the  region  occupied  by  the  dielectric  object,  and  by  Va  a  region 
outside  of  Vo  in  which  scattered  fields  are  measured.  We  divide  both  Yo  and  Va  into  a  number 
of  small  cells,  say,  N  cells,  and  denote  by  x.  and  y,-,  respectively,  the  centers  of  the  Jth  cell 
in  Va  and  Vo.  Let  be  a  N-dimensional  vector  representing  the  incident  field  in  the  N 
cells  inside  the  dielectric  body  Vo,  and  let  fa^*'  be  a  N-dimensional  vector  representing  the 
scattered  fields  in  the  N  cells  inside  the  measurement  region  Va.  Then  we  have  developed  the 
following  formula  to  compute  fa*"^^  from  the  incident  field 

lf„«>  =  G„(S  '  -  (23) 
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where  S,  Go,  and  Ga  are,  respectively,  dielectric  profile  of  the  scattering  body.  Green’s  matrix 
inside  Vo,  and  correlated  Green’s  matrix  between  Va  and  V<,;  is  the  incident  field  inside 
the  target,  which  is  known  and  readily  obtainable  in  analytic  form.  Note  that,  since  we  are 
dealing  with  vector  fields,  each  element  of  the  N-dimensional  vectors  is  itself  a  3-dimensional 
vector,  and  each  element  of  the  N  x  N  matrices  is  a  3  x  3  matrix.  Explicitly,  the  matrix  S  is 
given  by: 

Sij  -  ^ [e(y,)  -  e.]  6,  (y.cVo)  (24) 


As  to  Go  and  Ga,  we  have  been  successful  in  developing  analytical  formulas  for  these  matrices, 
and  thereby  avoided  making  any  Bom  approximation  or  linear  parameterization.  Denoting  by 
km  the  wave  number  in  the  background  medium,  a  the  radius  of  a  sphere  of  volume  equ^  to 
that  of  the  cells,  y,  and  yj,  respectively,  the  centers  of  the  ith  and  yth  cells  inside  Vo,  and 
X,  the  center  of  the  ith  cell  in  Va,  the  results  are,  for  the  matrix  Go: 


i=y- 


i*y- 

\ 


(Go),  =  l  J 


(y„  y;  c  Vo) 


(Go),  =  a  ^^P^|*^yy|‘'y^'^-[cos(kma)  -  |^•sin(kma)]  (y.,  y/cVo) 

fn  .  1  I  1  ...1  (-V.-y/)(y.-y/).,  3  ,  3  , 

^  ikmlyry/l  (/kmly.-y;l)^  lyi-yj-lyryj  ikmly.-yyl  (dcmlyry/l)^  ^ 


(25) 


and,  for  the  correlated  matrix  Ga  between  Va  and  Vo: 


(Ga),  =  a  ^^P^|*^|^|‘'y^'^-[cos(k.nQ)  -  p^-sin(kma)]  (x,  cVa,  y;  e  Vo) 


(26) 


x(  fn  -  ^  .  I  .  ^  1  -  (x--y>)(x<-y>)fi  3  3  , 

^  /knlxryj  (/k„lx.-y;l)^  ^  lx.-3ol-lx.-i)l  ^  /kmlx.-y;l  (ikjx.-yjiy  ^ 


Details  of  description  of  these  quantities  and  derivation  of  the  formulas  are  given  in 
Appendix  C. 


Microwave  Dosimetry  and  Quantitative  Microwave  Imaging  by  Matrix  Inversion 

The  Green’s  matrix  approach  described  above  also  allowed  us  to  solve  the  second  aspect 
of  dosimetry  problem,  viz.,  the  determination  of  the  electric  field  inside  the  target  from 
measurement  of  the  scattered  field  in  a  limited  region  outside  the  target.  We  have  derived  the 
following  equations  relating  the  scattered  fields  in  Vo  and  in  Va  to  the  total  fields  inside  the 
dielectric  body,  Vo: 


lfo^'^>  =  GoSlfo> 

(27) 

|f„w>  =  G„Slfo> 

(28) 
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where  f<,  is  a  N -dimensional  vector  representing  the  total  fields  in  the  N  cells  of  the  dielectric 
body.  Applying  the  inverse  of  Ga  to  both  sides  of  the  second  equation  and  substituting  the 
result  to  the  first  equation,  and  then  adding  the  incident  fields  to  both  sides,  one  gets: 


lfi>  =  lf/>>  +  GoGaW^^>  (29) 

which  gives  the  total  field  inside  the  dielectric  body  in  terms  of  the  measured  scattered  field 
in  any  region,  say  Va,  outside  the  body.  With  this  formula,  we  have  accomplished  the  principal 
goal  of  non-invasive  dosimetry. 

The  above  approach  may  also  be  extended  to  microwave  imaging.  Once  the  field  is  known 
inside  the  target  Vo,  one  may  then  derive  the  dielectric  profile  inside  the  target.  To  this  end, 
we  first  derived  the  following  relationship: 


S[lfo«>  +  GoGa-'lfa'^J>]  =  Ga'lfaf^^> 


(30) 


from  which  we  obtained  the  dielectric  profile  inside  the  target  in  terms  of  the  scattered  field 
outside  the  target: 


Si,  =  ~[e(.v,)-e«]  = 

C-m 


2  (GaW>(X,) 

-i-  2  (Go),/(Ga  O;*foi^'Xxt) 


(31) 


The  above  equation  gives  the  dielectric  permittivity  of  the  /“*  cell,  Su,  in  the  target  in  terms 
of  the  scattering  fields  in  any  scattering  space,  Va,  outside  the  scatterer,  and  the 

incident  fields  inside  the  scatterer.  The  incident  fields  inside  the  scatterer,  fo^'^X;),  are 
known  analytically. 

There  is  still  one  problem  that  needs  further  study.  It  concerns  the  data  stability  of  the 
inverse  of  the  matrix  Ga-  Our  dosimetry  algorithm,  as  described  by  the  above  equation, 
requires  computation  of  the  inverse  of  Go.  If  this  matrix  is  almost  singular,  then  computation 
of  its  inverse  may  be  unstable  with  respect  to  noises  and  errors  in  the  data  of  its  matrix 
elements.  Therefore,  further  effort  on  non-invasive  dosimetry  must  be  focused  on  developing 
techniques  to  stabilize  the  inverse  of  Go. 


Microwave  Dosimetry  and  Qualitative  Microwave  Imaging  by  "Soft  Focusing" 

The  objective  of  inverse  scattering  is  to  reconstruct  the  target  from  the  scattered  field. 
From  Maxwell’s  electromagnetic  theory,  if  one  knows  the  scattered  field  everywhere  in  space, 
the  source  of  the  field,  which  is  the  induced  charge-current  distribution  in  the  target,  can  be 
derived  completely.  However,  in  practice,  one  can  only  measure  the  scattered  field  exterior  to 
the  target  and  only  at  a  limited  number  of  points  in  space  which  are  often  confined  in  a  small 
region.  The  question  is  then  how  much  information  on  the  scattering  target  one  can  infer 
based  on  a  limited  knowledge  of  the  scattered  field.  Here  we  remark  that  the  term  "scattering 
field"  differs  from  the  conventional  definition  in  the  sense  that  it  refers  to  a  field  anywhere 
outside  of  the  target;  it  does  not  need  to  be  far  away  from  the  target.  Indeed,  the  fields  for 
image  reconstruction  must  be  close  to  the  target  in  order  to  contain  information  on  the  depth 
of  the  target. 


Page  14 

Transient  Interaction  of  Electromagnetic  Pulses  in  Dielectrics  and  Microwave  Biophysics 


5/88  Report  to  WRAK  via  ONR  Contract  No.  N00014-85-K -0475  The  Catholic  University  of  America  -  Theodore  C.  Guo  and  Wendy  W.  Guo 


Our  approach  to  inverse  scattering  with  a  limited  number  of  data  is  to  use  a  soft  focusing 
technique  to  focus  the  scattering  data  back  to  the  target.  Mathematically,  soft  focusing  is 
similar  to  regular  focusing  of  radiation,  which  we  call  "hard  focusing".  In  hard  focusing,  the 
substance  is  some  sort  of  radiation,  the  tool  is  a  lens,  and  the  focal  point  is  where  the 
radiation  energy  converges.  In  parallel,  the  substance  of  soft  focusing  is  a  set  of  information 
data,  the  focusing  tool  is  some  algorithm,  and  the  focal  point  is  where  the  information 
consolidates.  The  soft  focusing  technique  is  based  on  an  inverse  scattering  theorem  that  we 
have  formulated;^-^^ 

!“///[  ^  ^  ^  ^  ?  J.-E.«„(x,)  (32) 

Vo 

where  the  quantities  are  defined  as  below; 
c  Speed  of  light  in  vacuum. 

Vo  Space  occupied  by  the  target. 

Em,  Xm  Dielectric  susceptibility  and  dielectric  permittivity,  respectively,  of  water,  which 

is  the  background  medium. 

X  -  /C(^)  Dielectric  susceptibility  of  the  target. 

ko,  =  Velw/c  Wave  number  of  the  probing  microwave  in  the  background  medium. 

(o/ln  Microwave  frequency. 

P  =  P(x)  Induced  dielectric  polarization  in  the  target. 

Xo  Coordinate  of  the  n**  element  of  the  receiving  array. 

E«oh(Xo)  Electric  field  of  the  scattered  wave  measured  by  the  n^*  element  of  the 
receiving  array. 

Jo  Weighing  factor  to  be  applied  to  the  n**  element  of  the  array. 

Aw  =  Aw(x)  Vector  field  that  would  be  produced  by  a  set  of  current  elements  equivalent  to 
the  set  of  weighing  factors  (J,). 

The  theorem  described  above  may  be  considered  as  a  generalization  of  Lorentz  reciprocity 
theorem.-’^'^-’-^^  So,  multiplying  the  measured  scattered  field  at  each  of  the  array  element 
with  a  weighing  factor,  then  the  sum  of  the  products  is  c/e*.  times  the  sum  of  the  integrals  of 
[(l-Xm/x)P]  weighted  by,  respectively,  -V-Aw  and  ICn^Aw-.  The  weighing 
field  Aw  is  equal  to  the  vector  field  that  would  be  produced  by  a  set  of  current  distribution 
equal  to  the  weighing  factors,  {J„},  so  it  is  given  by; 


Aw(x)  = 


V  exp(ik«lx  -  x,l)  , 
^  —  Tw  r'l - J" 


(33) 


The  integral  on  the  left  hand  side  of  eq.  32  is  to  integrate  over  only  the  target,  therefore  it 
will  not  be  affected  by  the  values  of  Aw  outside  the  target.  This  gives  us  a  free  hand  to 
select  the  weighing  factors  {J„}.  To  retrieve  dielectric  property  of  the  target  at  a  focal  point, 
say  Xf ,  one  then  finds  a  set  of  weighing  factors  { J«)  such  that  the  corresponding  vector  field 
Aw  has  a  sharp  peak  at  X/  ,  and  negligibly  small  elsewhere  inside  the  target.  We  have  found 
that  this  is  possible  by  selecting  the  magnitude  of  the  vector  J,  to  be  the  phase-amplitude 
conjugation  corresponding  to  the  path  from  the  element  x,  to  the  desired  focal  point,  x/ : 

IJ„(x/)l  =  clx,  -  X/l-exp(-/lCnlx«  -  x/l)  (34) 
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As  to  the  direction  of  J„,  if  only  one  direction  of  the  scattered  field  is  measured  by  the  array 
elements,  then  may  be  chosen  to  be  in  that  direction  only,  so  that  the  sum  on  the  right 
hand  side  of  eq.  32  can  be  calculated  from  the  measured  field. 

To  summarize  our  approach  of  soft  focusing,  one  first  acquires  the  scattered  fields  at 
each  of  the  array  elements  at  {x,},  a  qualitative  dielectric  property  of  the  target  at  any  point, 
say  X/,  is  given  by: 


Qualitative  dielectric  property  at  x/  =  Z  J«-EjM«(Xn)  (35) 

n 


where  the  magnitude  of  J»  is  given  by  eq.  34,  and  its  direction  is  taken  to  be  that  of  the  field 
being  measured.  If  only  one  direction,  say,  y-direction,  of  Ewan  is  measured,  then  take  Jh  to 
be  only  in  y-direction,  so  that  only  the  measured  part  of  the  field  enters  the  above  equation. 
The  set  of  factors  {J„)  works  like  a  synthetic  soft  lens  for  focusing  the  scattering  data,  and 
the  weighing  field,  Aw(x),  gives  the  equivalent  field  pattern  of  the  soft  lens.  Therefore,  the 
quality  of  the  soft  lens  may  be  evaluated  by  analyzing  the  field  pattern  of  Ah-(x).  Note  that, 
given  any  lens  {J»},  the  corresponding  field  Ah.(x)  will  have  all  sorts  of  peaks  outside  some 
finite  region.  These  peaks  will  contribute  to  the  integral  on  the  left  hand  side  of  eq.  32 
unless  they  are  outside  the  target.  Therefore,  some  a  priori  knowledge  on  the  geometrical 
extent  of  the  target  is  a  necessary  condition  of  the  applicability  of  the  above  "soft  focusing" 
method. 

We  have  made  numerical  computation  and  modeling  to  analyze  the  sensitivity  of  the  above 
"soft  focusing"  technique.  The  computation  is  based  on  a  water-immersed  medical  imaging 
system  that  we  studied  for  the  Walter  Reed  Army  Institute  of  Research  (WRAIR).  The 
receiving  antenna  has  hexagonal  lattice  structure  with  127  waveguide-fed  antenna  elements, 
each  of  size  4  mm  x  7  mm  (see  Figures  la-Ic  of  Appendix  E  or  reference  5).  The  system 
operates  at  3  GHz  in  water.  These  results  show  a  3  dB  focusing  resolution  of  about  5  mm  in 
the  transverse  direction,  and  1 1  mm  in  the  longitudinal  direction  (see  Figures  2-5  and  Table  1 
of  Appendix  E  or  reference  5). 
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INTRODUCTION 

Acoustic  waves  generated  by  pulsed  microwaves  have  been 
cited  as  a  mechanism  for  microwave  hearing  tl],[2],[3]. 
More  recently  it  has  been  demonstrated  that  acoustic 
waves  are  transduced  in  dielectric  objects  simulating  the 
ocular  lens  when  exposed  to  pulsed  microwaves  [4]  .  The 
effect  has  also  been  cited  as  the  operant  mechanism  for 
cellular  damage  in  studies  of  the  murine  ocular  lens  in 
vitro  [5] .  Three  theories  on  the  microwave  auditory 
mechanism  have  been  reviewed  and  con^sared  by  Lin  [6], 
viz.,  the  radiation  pressure  theory,  the  electrostrictive 
theory,  and  the  thermoacoustic  theory.  Among  them  it  was 
found  that  only  the  thermoacoustic  mechanism  could 
produce  elastic  waves  of  magnitude  large  enough  to 
explain  the  experimental  observations.  Comparing  to  the 
magnitude  of  acoustic  pressure  generated  in  typical 
biological  tissues,  the  electrostrictive  effect  was  found 
to  be  about  two  orders  of  magnitude  smaller,  while  the 
radiation  pressure  three  orders  of  magnitude  smaller  [7] . 

Previous  theories  on  microwave  electrostrictive  process 
and  microwave  thermoacoustic  process  are  based  on  a 
linear  wave  equation  with  a  generating  function  derived 
from  the  microwave  pulses .  Pressure  waves  are  thereby 
generated  in  a  way  similar  to  a  forced  harmonic 
oscillator.  In  the  case  of  the  thermoacoustic  process, 
the  generating  function  is  derived  from  an  inhomogeneous 
heating  by  the  microwaves,  whereas  in  the  case  of  the 
electrostrictive  process,  it  is  derived  from  stress  due 
to  dielectric  polarization  [8], [9].  This  approach  fails 
to  consider  the  proper  balance  of  the  distribution  of  the 
absorbed  microwave  energy  among  the  internal  thermal 
energies  and  the  bulk  kinetic  energy,  and  the  effect  of 
thermal  or  dielectric  discontinuities  at  medium 
interface.  Also  neglected  is  the  coupling  between  the 
thermoacoustic  process  and  electrostrictive  process, 
which  may  result  from  dependency  of  dielectric 
permittivity  on  various  thermodynamical  coefficients. 


152 


Reprint  from  the  Annual  Report  of  the  1967  ICEC  Conference  on  Electrical  Insulation  and  Dielectric  Phenomena. 


Dielectric  Acoustic  Response  to  Microwave  Pulses 
Through  Thermoacoustlc  and  Electrostrlctive  Effects 


T.  C.  Guo  and  W.  W.  Guo 


such  as  mass  density,  temperature,  and  pressure;  this 
coupling  may  further  enhance  the  electrostrictive  effect. 

Earlier,  we  have  developed  a  thermodynamical  formulation 
of  the  process  of  microwave  absorption  and  its  coupling 
to  the  internal  thermal  energy  and  mechanical  bulk  energy 
in  the  dielectrics  [10]  .  It  was  shown  that  acoustic 
waves  may  be  generated  by  microwave  pulses  if  there  are 
discontinuities  in  thermodynamical  characteristics  in  the 
dielectrics,  or  discontinuities  in  microwave  absorptions. 
Numerical  results  for  a  one-dimensional  system  and  a 
spherical  system  with  water  as  the  absorbing  dielectric 
have  been  compared  consistently  with  experimental 
results,  which  included  measurement  of  pressure  waves  in 
a  cylinder  filled  with  an  aqueous  solution  of 
electrolytes  exposed  to  pulsed  microwaves  and  estimation 
of  the  coupling  efficiency  between  the  liquid  dielectric 
and  the  adjacent  air  [11] .  However,  these  studies  did 
not  include  electrostrictive  force  in  the  dielectric  and 
the  dielectric  interface.  In  this  paper,-  we  reformulate 
the  problem  to  include  electrostrictive  stress  and  show 
that  it  may  play  an  important  role  in  the  generation  of 
pressure  waves  by  microwave  pulses  in  dielectrics  with 
thermal  or  dielectric  discontinuities.  We  also  show 
that,  in  the  first-order  linear  approximation,  the  effect 
of  the  electrostrictive  force  on  the  acoustic  velocity  is 
negligible . 


THEORY  OF  MICROWAVE  ACOUSTICS 

Our  approach  to  derive  the  coupling  from  microwave  pulses 
to  acoustic  waves  is  to  start  from  the  laws  of  the 
conservation  of  mass,  the  conservation  of  momentum,  the 
conservation  of  energy,  and  the  thermodynamic  equation  of 
state.  When  a  dielectric  system  is  exposed  to  microwave 
radiation,  microwave  energy  is  absorbed  into  the 
dielectric  through  the  work  of  the  electric  field  on  the 
electric  polarization:  E-dP,  which  is  then  partially 
converted  into  the  internal  energy  Ui  and  partially  into 
the  bulk  kinetic  energy  Uk  via  thermal  expansion  and 
other  elastic  vibrations.  The  process  is  complex  and 
irreversible;  however,  we  shall  assume  that  the  time 
domain  which  concerns  us  is  much  greater  than  the  time  it 
tckes  to  reach  a  quasi-steady  state  and  that  every 
microscopically  small  part  of  the  system  is  in  local 
ec[uilibrium  so  that  thermodynamic  variables  may  be 
defined. 
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The  time  domain  of  our  interest  is  the  width  of  the 
microwave  pulse,  which  is  in  the  order  of  microseconds 
and  is  too  short  for  any  heat  conduction  and  convective 
flow  to  take  place.  On  the  other  hand,  the  time  cycle  of 
electric  field  is  in  the  range  of  subnanosecond,  which  is 
much  shorter  than  the  time  domain  of  our  concern. 
Therefore,  we  may  average  all  variables  over  a  microwave 
cycle.  The  resulting  time-average  of  E*dP  gives  the 
microwave  energy  absorbed  by  the  dielectric,  SQp ,  and  the 
electrostrictive  work  on  the  dielectric,  8Wp .  Thus,  the 
first  law  of  thermodynamics  gives: 

5Qp  +  5Wp  =  d(Ui  +  Uk)  +  8W  (1) 

where  8W  is  the  mechanical  work  done  by  the  dielectric; 
heat  flow  by  conduction  and  convection  have  been  ignored, 
as  explained  earlier.  Denote  by  pi j  the  stress  tensor  in 
the  dielectric,  and  by  Oi  j  the  electrostrictive  stress 
tensor  on  the  dielectric  by  the  electromagnetic  field, 
then  J||3jpijdV  and  JJj3jCijdV  are,  respectively,  the 
force  on  the  dielectric  volume  element  by  the  surrounding 
body  and  the  force  on  the  element  by  the  electromagnetic 
field.  The  total  stress  tensor  experienced  by  the 
dielectric  is  then  the  sum  of  pi;  and  di;,  which  will  be 
denoted  by  si  ; ; 


Si  j  =  Pi  j  +  ai  j 


If  e  is  the  dielectric  permittivity,  then  [12] 


(e  -  P(|^)t]5i, 


471 


EiE) 


(2) 


Averaging  over  the  microwave  cycle,  the  terms  EiE;  vanish 
except  for  i=j,  and  the  above  equation  give  : 


OiJ 


r 

L  167C 


[e 


9e 


(3) 


In  order  to  convert  eq.  (1)  into  a  differential  equation 
with  respect  to  space  and  time,  we  start  from  a  general 
form  of  the  equation  of  conservation.  Let  g  be  a 
themoclynamic  extensive  quantity,  Sg  be  the  rate  of 
production  of  g,  which  may  be  due  to  internal  production 
or  external  input,  and  v  be  the  bulk  velocity  function, 
then  the  equation  of  conservation  for  g  is: 
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atg  +  gV-v  =  Sg  (4) 

where  we  have  ignored  the  second  order  term  vVg  by 
assuming  that  both  g  and  v  are  small.  The  first  term  on 
the  left  represents  the  rate  of  increase  of  g  per  unit 
volume,  and  the  second  term  is  the  rate  of  outgoing  flow 
of  g.  Taking  g  to  be,  respectively,  the  mass  density, 
the  stress  tensor  Sij,  and  energy  density,  eq.  (4)  gives 


atp  +  pV-v  =  0  (5) 

p3cvi  -  d)Sij  =  0  (6) 

patUi  -  Sijdjvi  =  pp(x,t)  (7) 

where  P  is  the  specific  absorption  rate  (SAR)  of 
microwaves,  and  Ui  is  the  internal  energy  density.  Short 
of  an  analytic  expression  for  the  internal  energy,  we 
shall  express  atUi  in  terms  of  empirical  coefficients. 
In  terms  of  the  heat  capacities  Cp  and  Cv ,  isothermal 
compressibility  tensor  (Kt)ij,  and  isobaric  thermal 
expansion  coefficient  Pp,  it  may  be  shown  that: 


3t  Ui 


■-G 


Pp 


Pi  jdjVi' 

3tP 

p^ V-v  . 

P 

r(e-l) 

•(KT)jjatPij  +  dt  [ — g- — 


Note  that,  while  the  total  stress  on  the  dielectric  is 
sjj  -  Pi  j  +  Oij,  it  is  the  material  stress  pj  j  which  is 
responsible  for  the  material  strain,  and  therefore 
accountable  in  the  thermal  internal  energy.  Indeed,  in 
the  presence  of  an  electromagnetic  field,  the  material 
stress  automatically  adjusts  itself  in  reaction  to  the 
electrostrictlve  stress.  Therefore,  the  above  equation 
does  not  include  Oij.  To  obtain  atUi  in  terms  of  dtSij, 
we  replace,  in  the  above  equation,  pi j  by  sij-Oij  and  use 
eq.  (3)  for  CTi j ,  and  replace  V-v  by  -at  p/p.  It  gives; 


3rT.  _  _r  Cp1^_,^Slj)ajVl  Cv(KT)ij).,_ 

— d - 3; — 
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Eqs .  (5)-('.)  are  the  equations  of  the  system. 
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THE  LIHEAR  APPROXIMATION 


The  above  system  equations  may  be  mutually  substituted  to 
eliminate  the  terms  involving  ▼  and  p  and  thereby  yield  a 
wave  equation  for  si  j .  However,  it  is  rather  difficult 
to  eliminate  the  term  r(|Eil^9ivi)  (summing  over  i)  in 
eq.  (8),  which  is  highly  anisotropic.  One  may  select  a 
coordinate  system  such  that  |Ei 1  in  all  three  directions 
are  equal,  but  then  it  will  complicate  the  boundary 
conditions.  Nevertheless,  we  shall  make  a  crude 
approximation  that  |Ei | 1E| ^ /3 .  Then  £(|Ei|^3iVi)  = 

(1/3) lEl^V-v,  which,  by  eq.  (5),  is  -  (1/3) IE | ^ (dt p/p) . 
With  this  approximation,  the  electrostrlctive  tensor  Oj  j 
in  eq.  (3)  becomes  isotropic:  Oij  =  -CT-Sij,  where  o  is 
the  electrostrlctive  pressure  defined  as: 


o 


|E|" 

16jt 


>1 


O'l?’'  '  ] 


(9) 


Neglecting  terms  involving  second  order  products  of  P, 
Oij,  |E|,  and  derivatives  of  all  thermal  quantities,  we 
arrive  at  the  following  wave  equation  for  the  stress 
tensor  si j  : 
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where  the  subscript  o  signifies  that  the  quantity  is 
evaluated  at  its  ambient  equilibrium  value.  For  an 
isotropic  dielectric,  sij  =  -s&ij  and  (KT)i_;  =  Kt/3,  so 
2(kT)ii  (sura  over  /)  is  Kt,  the  above  wave  equation  is 
then  further  reduced  to 
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where  Uo  denotes  the  acoustic  velocity  in  the  dielectric 
in  the  absence  of  the  electromagnetic  waves.  To  estimate 
the  magnitude  of  the  change  of  acoustic  velocity  by 
electrostrictive  pressure  c,  we  consider  water  at  room 
temperature  as  the  dielectric,  for  which  p  =  1  gm/cm^, 
Cp  =  4.186x10''  erg/gm*C,  pp  =  2.8x10"'  *C‘',  and  E  =  80. 
We  also  assume  that  the  dielectric  permittivity  is 
linearly  proportional  to  the  mass  density,  so  Oe/3p)T  = 
e/p,  then  eq.  (9)  gives  o  =  (e  |E  ( ^ )  /  (247C)  .  Assume  a 
plane  electromagnetic  wave  of  intensity  1  =  1  kW/cm^ , 
then  its  field  magnitude  is  |E|^  =  8x1/ (cVe)  = 

0.94  erg/cm^ ,  where  c  is  the  speed  of  light  in  vacuum. 
Then  0  =  1  dyne/cm^,  and  the  correction  term  is 
(aPp)/(pCp)  =  6.7x10"'^,  which  is  insignificant. 

PRESSURE  NAVE  GENERATION  ‘ 


AS  illustrated  at  the  end  of  the  last  paragraph,  for  most 
dielectrics  and  for  microwave  intensity  up  to 
gigawatts/cm^ ,  the  term  a/p  on  both  sides  of  eq.  (11)  may 
be  neglected.  Eq.  (11)  may  then  be  rewritten  as 


V^s  -  ' 

11* 


s  - 


pp 

Cp 


3c 


(13) 

where  u  =  (Cp /pCv  kt  ]  ^ ^  is  the  speed  of  sound  in  the 
dielectric.  The  three  terms  inside  the  brackets  on  the 
right  hand  side  are  the  generating  forces  of  pressure 
waves.  The  first  term  represents  the  thermoelastic 
effect,  the  second  term  is  the  electrostrictive  force, 
and  the  third  term  comes  from  the  electromagnetic  energy 
of  dielectric  polarization.  We  shall  compare  the  sizes 
of  these  three  driving  forces  by  expressing  them  in  terms 
of  |E|^.  Let  a  be  the  attenuation  coefficient  of 
microwave  intensity  as  defined  by  the  equation  I  (x)  = 

Ioe"“*  .  Then  the  first  term  is  poP  =  ai  =  (acVe/87t)  |E1^. 
If  the  dielectric  permittivity  is  linear  in  p,  then 
Oe/3p)T  “  e/p,  so  eq.  (9)  gives  a  =  (e  |E|  ^ )  /  (247t)  and 
the  electrostrictive  term  gives  ( (peCvKi ) / (24rtPp ) ] 3t  IE | ^ . 
Denote  by  t,  the  rise  time  of  a  microwave  pulse,  then 
3tlE|^  ~  lEI^/tr.  Thus,  the  relative  magnitudes  of  these 
three  terms  are,  respectively. 


The  first  paragraph  in  this  section  was  inadvertently 
missed  from  the  manuscript  submitted  to  CEIDP. 
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One  sees  that  the  second  and  the  third  driving  forces  may 
be  larger  than  the  thermoelastic  term  if  the  microwave  is 
rapidly  modulated,  contrary  to  the  estimate  by  previous 
theory  that  the  electrostrictive  effect  is  two  orders  of 
magnitude  smaller  (61.  For  water  at  30  'C,  p  =  1  gm/cm^ , 
Cp  =  4.186x10’  erg/gm’C,  Kt  =  4.46x10'”  cm^/dyne,  Pp  = 
2.8x10'^  ‘C"’,  and,  at  3  GHz,  a  =  0.88  cm"’  and  e  ■=  80. 
The  relative  magnitudes  are,  respectively: 

,  Cv  0.75  nanosec  0.34  nanosec 
1/  *  r  r  Z  •  (lo) 

L.p  cr  tr 

One  may  conclude  that,  for  water,  a  rise  time  of 
1  nanosecond  is  the  marginal  modulation  rate,  beyond 
which  the  electrostrictive  term  and  the  polarization  term 
contribute  more  to  the  pressure  waves  than  the 
thermoelastic  term. 


In  the  absence  of  the  electrostrictive  force,  the  wave 
equation  for  the  pressure  waves  reduces  to  the  one 
developed  for  thermoacoustic  process  [10,11].  While  we 
have  not  yet  completed  calculations  on  the  pressure  waves 
generated  by  the  electrostrictive  effect,  pressure  waves 
from  the  thermoacoustic  effects  have  been  reported 
[10,11],  and  measurement  on  the  amplitude  and  pulse 
characteristics  have  also  been  made  [11]  .  Contrary  to 
previous  theory  that  thermoacoustic  pressure  waves  were 
generated  mainly  by  inhomogeneous  distribution  of 
microwaves  in  dielectrics  [6],  it  was  found  that  pressure 
waves  are  generated  whenever  there  is  a  discontinuity  in 
thermal  or  dielectric  parameters  in  the  medium. 
Equations  (10)  and  (11)  also  indicate  this  effect  for 
electrostrictive  pressure  waves.  With  an  air-water 
system  and  neglecting  the  electrostrictive  force,  it  was 
found  that  a  pea)c  SAR  of  4.5  IcW/gram  sufficed  to  generate 
a  pressure  wave  of  pea)t  pressure  equivalent  to  about  10% 
of  the  ambient  atmospheric  pressure,  which  is  the  upper 
limit  of  linear  approximation. 
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INTRODUCTION 

Acoustic  waves  generated  by  pulsed  microwaves 
have  been  cited  as  a  mechanism  for  microwave 
hearing  (1,2,3).  More  recently  it  has  been 
demonstrated  that  acoustic  waves  are  trans¬ 
duced  in  dielectric  objects  simulating  the 
ocular  lens  when  exposed  to  pulsed  microwaves 

(4)  .  The  effect  has  also  been  cited  as  the 
operant  mechanism  for  cellular  damage  in 
studies  of  the  murine  ocular  lens  in  vitro 

(5)  .  Three  theories  on  the  microwave  auditory 
mechanism  have  been  reviewed  and  compared  by 
Lin  (6),  viz.,  the  radiation  pressure  theory, 
the  electrostrictive  theory,  and  the  ther¬ 
moacoustic  theory.  Among  them  it  was  found 
that  only  the  thermoacoustic  mechanism  could 
produce  elastic  waves  of  magnitude  large 
enough  to  explain  the  experimental  observa¬ 
tions.  Comparing  to  the  magnitude  of  acoustic 
pressure  generated  in  typical  biological 
tissues,  the  electrostrictive  effect  was  found 
to  be  about  two  orders  of  magnitude  smaller, 
while  the  radiation  pressure  three  orders  of 
magnitude  smaller  (7)  . 

Previous  theories  on  microwave  electrostric¬ 
tive  process  and  microwave  thermcacoustic 
process  are  based  on  a  linear  wave  equation 
with  a  generating  function  derived  from  the 
microwave  pulses.  Pressure  waves  are  thereby 
generated  in  a  way  similar  to  a  forced  har¬ 
monic  oscillator.  In  the  case  of  the  ther¬ 
moacoustic  process,  the  generating  function  is 
derived  from  an  inhomogeneous  heating  by  the 
microwaves,  whereas  in  the  case  of  the  elec¬ 
trostrictive  process,  it  is  derived  from 
stress  due  to  dielectric  polarization  (8,9). 
This  approach  fails  to  consider  the  proper 
balance  of  the  distribution  of  the  absorbed 
microwave  energy  among  the  internal  thermal 
energies  and  the  built  Itinetic  energy,  and  the 
effect  of  thermal  or  dielectric  discontinu¬ 
ities  at  medium  Interface.  Also  neglected  is 
the  coupling  between  the  thermoacoustic  pro¬ 
cess  and  electrostrictive  process,  which  may 
result  from  dependency  of  dielectric  permit¬ 
tivity  on  various  thermodynamical  coeffi¬ 
cients,  such  as  mass  density,  temperature,  and 
pressure;  this  coupling  may  further  enhance 
the  electrostrictive  effect. 

Earlier,  we  have  developed  a  thermodynamical 
formulation  of  the  process  of  microwave  ab¬ 
sorption  and  its  coupling  to  the  Internal 
thermal  energy  and  mechanical  built  energy  in 
the  dielectrics  (10)  .  It  was  shown  that 
acoustic  waves  may  be  generated  by  microwave 
pulses  if  there  are  discontinuities  in  ther¬ 
modynamical  characteristics  in  the  dielec¬ 
trics,  or  discontinuities  in  microwave  absorp¬ 
tions.  Numerical  results  for  a  one-dimen¬ 
sional  system  and  a  spherical  system  with 
water  as  the  absorbing  dielectric  have  been 
compared  consistently  with  experimental  re¬ 
sults,  which  included  measurement  of  pressure 
waves  in  a  cylinder  filled  with  an  aqueous 
solution  of  electrolytes  exposed  to  pulsed 
microwaves  and  estimation  of  the  coupling 
efficiency  between  the  liquid  dielectric  and 


the  adjacent  air  (11)  .  However,  these  studies 
did  not  include  electrostrictive  force  in  the 
dielectric  and  the  dielectric  interface.  In 
this  paper,  we  reformulate  the  problem  to 
include  electrostrictive  stress  and  show  that 
it  may  play  an  important  role  in  the  genera¬ 
tion  of  pressure  waves  by  microwave  pulses  in 
dielectrics  with  thermal  or  dielectric  discon¬ 
tinuities.  We  also  show  that,  in  the  first- 
order  linear  approximation,  the  effect  of  the 
electrostrictive  force  on  the  acoustic  veloc¬ 
ity  is  negligible. 

THEORY  OF  MICROWAVE  ACOUSTICS 

Our  approach  to  derive  the  coupling  from 
microwave  pulses  to  acoustic  waves  is  to  start 
from  the  laws  of  the  conservation  of  mass,  the 
conservation  of  momentum,  the  conservation  of 
energy,  and  the  thermodynamic  equation  of 
state.  When  a  dielectric  system  is  exposed  to 
microwave  radiation,  microwave  energy  is 
absorbed  into  the  dielectric  through  the  worlc 
of  the  electric  field  on  the  electric  polari¬ 
zation:  E-dR,  which  is  then  partially  con¬ 
verted  into  the  internal  energy  Ui  and  par¬ 
tially  into  the  bul)c  jcinetic  energy  Uk  via 
thermal  expansion  and  other  elastic  vibra¬ 
tions.  The  process  is  complex  and  irrever¬ 
sible;  however,  we  shall  assume  ..hat  the  time 
domain  which  concerns  us  is  much  greater  than 
the  time  it  talces  to  reach  a  quasi-steady 
state  and  that  every  microscopically  small 
part  of  the  system  is  in  local  equilibrium  so 
that  thermodynamic  variables  may  be  defined. 

The  time  domain  of  our  interest  is  the  width 
of  the  microwave  pulse,  which  is  in  the  order 
of  microseconds  and  is  too  short  for  any  heat 
conduction  and  convective  flow  to  take  place. 
On  the  other  hand,  the  time  cycle  of  electric 
field  is  in  the  range  of  subnanosecond,  which 
is  much  shorter  than  the  time  domain  of  our 
concern.  Therefore,  we  may  average  all  vari¬ 
ables  over  a  microwave  cycle.  The  resulting 
time-average  of  E-dP  gives  the  microwave 
energy  absorbed  by  the  dielectric,  6Qp  ,  and 
the  electrostrictive  work  on  the  dielectric, 
6Wp .  Thus,  the  first  law  of  thermodynamics 
gives : 

6Qp  +  6Wp  =  d(Ui  +  Uk)  +  6w  (1) 

where  6w  is  the  mechanical  work  done  by  the 
dielectric;  heat  flow  by  conduction  and  con¬ 
vection  have  been  ignored,  as  explained  ear¬ 
lier.  Denote  by  pi ^  the  stress  tensor  in  the 
dielectric,  and  )oy  ojj  the  electrostrictive 
stress  tensor  on  the  dielectric  by  the  elec¬ 
tromagnetic  field,  then  JJj3jpijdV  and 
njd^OijdV  are,  respectively,  the  force  on  the 
dielectric  volume  element  by  the  surrounding 
body  and  the  force  on  the  element  by  the 
electromagnetic  field.  The  total  stress 
tensor  experienced  by  the  dielectric  is  then 
the  sum  of  pjj  and  Oij,  which  will  be  denoted 
by  Sij: 

Slj  =  PiJ  +  Oij 
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If  €  is  the  dielectric  permittivity,  then  (12) 

"  §5 

Averaging  over  the  microwave  cycle,  the  terms 
EiEj  vanish  except  for  i=j,  and  the  above 
equation  gives: 

-  P#^  1  -  ‘3> 

In  order  to  convert  eq.  (2)  into  a  differen¬ 
tial  equation  with  respect  to  space  and  time, 
we  start  from  a  general  form  of  the  equation 
of  conservation.  Let  g  be  a  thermodynamic 
extensive  quantity,  Sq  be  the  rate  of  produc¬ 
tion  of  g,  which  may  be  due  to  internal  pro¬ 
duction  or  external  input,  and  v  be  the  bulk 
velocity  function,  then  the  equation  of  con¬ 
servation  for  g  is: 


It  is  remarked  that,  while  the  total  stress  on 
the  dielectric  is  sij  “  pij  +  Oij,  it  is  the 
material  stress  pij  which  is  responsible  for 
the  material  strain,  and  therefore  accountable 
in  the  thermal  internal  energy.  Indeed,  in 
the  presence  of  an  electromagnetic  field,  the 
material  stress  automatically  adjusts  itself 
in  reaction  to  the  electrostrictive  stress  by 
the  field,  therefore,  the  above  equation  does 
not  include  Oij. 

We  would  like  to  express  the  partial  deriva¬ 
tives  of  Uth  with  respect  to  p  and  pij  in 
terms  of  known  empirical  quantities  such  as 
specific  heat  and  compressibility.  Utilizing 
the  Jacobian  method  of  relating  one  partial 
derivative  to  another,  one  may  derive  from  eq. 
(9)  the  following  equation: 


dtUth  -  [  -  ^  ]atp  +  [  ]dtP  (10) 


dtg  +  gV-v  =  S,  (4) 

where  we  have  ignored  the  second  order  term 
vVg  by  assuming  that  both  g  and  v  are  small. 
The  first  term  on  the  left  represents  the  rate 
of  Increase  of  g  per  unit  volume,  and  the 
second  term  is  the  rate  of  outgoing  flow  of  g. 
Taking  g  to  be,  respectively,  the  mass  den¬ 
sity,  the  stress  tensor  Sij,  and  energy  den¬ 
sity,  eq.  (5)  gives 


3tp  t  pV-v  =  0 

(5) 

p9tvj  -  djsij  =  0 

<6) 

h  -  Si  jdjVi  =  pP(3C,t) 

(7) 

where  P  is  the  specific  absorption  rate  (SAR) 
of  microwaves,  and  Ur  is  the  internal  energy 
density  of  the  dielectric  system. 

Equations  (5) -(7)  all  together  contain  four 
unknowns,  viz.,  p,  v,  s,  and  Ui  .  One  more 
equation  is  needed  to  complete  the  system, 
which  will  be  furnished  by  the  thermodynamic 
equation  of  state: 

Ui  =  Uth  (p,  p)  +  (i)  <|^)  lEI’  (8) 

where  the  first  term  on  the  right  hand  side  is 
the  thermodynamical  energy,  whereas  the  second 
term  is  the  electrical  energy  in  the  dielec¬ 
tric  that  arises  from  electric  polarization. 
Since  the  electric  energy  density  averaged 
over  a  microwave  cycle  is  (l/8jt)  Re  (P-E)  = 
(1/871)  Re  (e-1)  |E  1 ,  e  in  eq.  (8)  represents  only 
the  real  part  of  the  dielectric  permittivity; 
the  imaginary  part  of  e  contributes  to  pP  in 
eq.  (7)  ,  the  energy  absorption  by  the  dielec¬ 
tric.  As  there  is  no  analytic  expression  for 
the  thermodynamical  energy,  Uth,  we  shall  try 
to  express  it  in  terms  of  other  empirical 
coefficients,  such  as  density  p,  specific 
heats  Cv  and  Cp,  isobaric  thermal  expansion 
coefficient  Pp  f  and  isothermal  compressibility 
tensor  (Kt)jj.  As  a  thermodynamical  system, 
the  dielectric  has  two  Independent  thermody¬ 
namic  variables,  which  may  be  any  pair  among 
the  variables  p,  T,  p,  and  Uth,  and  all  other 
variables  may  be  considered  as  functions  of 
the  pair.  Since  p  and  p  are  the  variables  in 
eqs.  (5)- (7),  they  are  selected  as  the  in¬ 
dependent  variables.  A  differential  of  Uth 
may  be  expressed  as: 

dUth  =  (dpUth)  dp  +  (d  Uth)pdpij  (9) 

F*  j  F-*  J 


where  Cv  and  Cp  are  the  specific  heat  (per 
mass)  at  constant  volume  and  constant  pres¬ 
sure,  respectively;  Pp  is  the  isobaric  thermal 
expansion  coefficient,  and  Kt  the  isothermal 
compressibility.  It  is  straightforward  to 
generalize  the  above  result  to  include  the 
stress  tensor.  Noting  that  the  sign  of  pi j  is 
generally  defined  such  that,  for  an  isotropic 
case.  Pi J  =  -p5i j ,  where  p  is  positive  when  it 
is  a  pressure,  and  negative  when  it  is  a 
tension.  Eq.  (10)  then  takes  the  form: 


^tOth  =  ^  ^  ]dtp  - 

(11) 

where  Cp  is  given  by: 


Cp  =  [  dtUth 


and  (KT)ij  is  the  compressibility  tensor 
defined  as 


In  deriving  the  above  equation,  we  have  used 
the  generalized  first  law  of  thermodynamics 
for  mechanically  anisotropic  media: 

6Q  =  dUth  -  dV. 

The  electromagnetic  terms  are  not  included 
because  the  thermodynamical  quantities,  Cv , 
Cp,  Pp,  and  are  generally  measured  in 
the  absence  of  electromagnetic  fields.  To 
understand  the  above  equation,  we  consider  8Q 
as  the  total  heat  input  to  a  material  element 
of  unit  mass  during  time  dt,  during  which  the 
total  energy  is  increased  loy  dU.  The  total 
work  ^one  by  the  element  is  equal  to 
-9-(V'p)dt  times  the  volume,  1/p.  On  the 
other  hand,  from  eq.  (5)  ,  the  increment  in 
volume  during  dt  is  dV  =  d(l/p)  =  (l/p)V-vdt. 
Therefore  the  work  done  associated  with  a 
volume  increment  dV  is  -[V- (v-p) /V-v]dV,  which 
may  be  separated  into  two  terms: 

- (Pi jdjvi /V-v) dV  -  (vidjpi j/V-v) dV, 

of  which  the  second  term  cancels  with  the  bulk 
kinetic  energy  Un,  so  only  the  first  term  is 
accounted  for  in  the  above  equation. 
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As  explained  earlier,  it  is  the  material 
stress  pij  which  is  responsible  for  the 
material  strain,  and  therefore  accountable  in 
the  thermal  internal  energy  Uth.  On  the  other 
hand,  the  acoustic  wave  equation  is  the  equa¬ 
tion  for  the  total  stress  Sij.  To  obtain 
3tUth  in  terms  of  dtsij,  we  replace,  in  the 
above  equation,  pij  by  sij-aij  and  substitute 
(3)  for  oij.  Eq.  (11)  then  gives: 


.  r  lEP  -,ae>  1  e  lEil^diVi  1  dtp 

^  L  TW  ■  P<5p>"  ’  ■  55 - 

+  ^(KT)iiat[  11^  (E  -  p(||)T  1  -  gllEil^  ] 


By  eq.  (5),  the  factor  V-v  in  the  denominator 
may  be  replaced  by  -dtp/p,  the  above  equation 
then  becomes; 


^  mlULL  -  ^(KiJijatSij 

*  [  155“  ■  P'lp’^  '  ]‘^  ^  alp' 

.  ^(KT)iiat[||^  (E  -  p(|f)T  1  -  ] 

(12) 

We  may  now  substitute  eq.  (12)  for  the  term 
Uth  in  eq.  (8)  to  obtain  an  expression  for 
dtUi  .  With  some  rearrangement  of  terms,  one 
obtains : 

■  t-  iife  •  If  -  'ff'-i 


(E-1) lEl^l 


+  ^(Ki)iiat[  j||-  (£  -  P(|^)t  )  -  jflEil^  ] 


Eqs .  (5) -(7)  and  (13)  are  the  equations  of  the 
system.  To  summarize,  we  list  these  equations 
below  (summation  over  repeated  indices) : 


dtp  +  pV'V  =  0 

pdtvi  -  djSij  =  0 

pdtUi  -  Sijdjvi  =  pP(x,t) 


[■ 


^  |E|^  ,E  .dE,  ,  1,  „ 

T55p  'p  ■  <5p>^’  J^’^p 


.  i^si,  - 


Cv(KT)ii 


sip' ] 

t  ^(KT)liat[  ||L  [E  -  P(||)t  1  -  ^lEll*  ] 


THE  LINEAR  APPROXIMATION 

The  above  system  equations  may  be  mutually 
substituted  to  eliminate  the  terms  involving  v 
and  p  and  thereby  yield  a  wave  equation  for 
Sij.  However,  it  is  rather  difficult  to 
eliminate  the  term  Z(lEil^3iVi)  (summing  over 
i)  in  eq.  (17),  which  is  highly  anisotropic. 
One  may  select  a  coordinate  system  such  that 
lEil  in  all  three  directions  are  equal,  but 
then  it  will  complicate  the  boundary  condi¬ 
tions.  Nevertheless,  we  shall  make  a  crude 
approximation  that  I  E.i  I  ^  =  IE  I  ^ /3 .  Then 

Z(IEil^divi)  “  (1/3)  IE  I ^ V-v,  which,  by  eq. 

(14),  is  -(1/3) IE|2 (dtp/p) .  With  this  approx¬ 
imation,  the  electrostrictive  tensor  Oi  j  in 
eq.  (3)  becomes  isotropic:  ai  j  -  -a-Sij,  where 
a  is  the  electrostrictive  pressure  defined  as: 

■>  ■  -  [  If  -  p'H"  '  ]  »•' 

Neglecting  terms  involving  second  order  pro¬ 
ducts  of  P,  Oij,  IE|,  and  derivatives  of  all 
thermal  quantities,  we  arrive  at  the  following 
wave  equation  for  the  stress  tensor  Sij: 

aj[  ^  djsij  ]  -  3t[  ^  5  ^(Kt)  jjdtsij] 


“  [  5^  ^  f 

X  [dipt  [fe(KT)il]>^a-di^[<i^]] 

where  the  subscript  o  signifies  that  the 
quantity  is  evaluated  at  its  ambient  equi¬ 
librium  value.  For  an  isotropic  dielectric. 
Si)  =  -siij  and  (KT)ij  =  Kt/3,  so  L(KT)ii  (sum 
over  i)  is  ict,  the  above  wave  equation  is  then 
further  reduced  to 

X  [p„diPt  [P^]^ai.a-di»[li^]] 

Equation  (20)  is  a  typical  .linear  wave  equa¬ 
tion  with  the  velocity 

U  =  fl  +  56£ll/2r  _Cp  11/2  ^  r  ^  opp11/2 

L  pCpJo  L  pCvKr  Jo  L  pCpJo  ° 


where  Uo  denotes  the  acoustic  velocity  in  the 
dielectric  in  the  absence  of  the  electromag¬ 
netic  waves.  To  estimate  the  magnitude  of  the 
change  of  acoustic  velocity  by  electrostric¬ 
tive  pressure  o,  we  consider  water  at  room 
temperature  as  the  dielectric,  for  which  p  = 
1  gm/cm’,  Cp  =  4.186x10''  erg/gm*C,  Pp 

2.8x10"^  ’C"',  and  E  =  80.  We  also  assume 
that  the  dielectric  permittivity  is  linearly 
proportional  to  the  mass  density,  so  (aE/dp)T 
=  E/p,  then  eq.  (18)  gives  O  =  (£  IE  I  ^  )  /  (24it)  . 
Assume  a  plane  electromagnetic  wave  of  inten¬ 
sity  1=1  kW/cm*,  then  its  field  magnitude  is 
|E1‘  =  8kI/(cV£)  =  0.94  erg/cm’,  where  c  is 

the  speed  of  light  in  vacuum.  Then  a  = 
1  dyne/cm* ,  and  the  correction  term  is 
(0Pp)/(pCp)  =  6.7x10“'^,  which  is  insig¬ 

nificant  . 
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As  illustrated  at  the  end  of  the  last 
paragraph,  for  most  dielectrics  and  for 
microwave  intensity  up  to  gigawatts/cra^ ,  the 
term  a/p  on  both  sides  of  eg.  (20)  may  be 
neglected.  Eg.  (20)  may  then  be  rewritten  as 


V2s  - 

(22) 


where  u  =  (pCvKr /Cp J ^ ^ ^  is  the  speed  of  sound 
In  the  dielectric.  The  three  terms  inside  the 
brac)iets  on  the  right  hand  side  are  the 
generating  forces  of  pressure  waves.  The 
first  term  represents  the  thermoelastic  ef¬ 
fect,  the  second  term  is  the  electrostrictive 
force,  and  the  third  term  comes  from  the 
electromagnetic  energy  of  dielectric  polariza¬ 
tion.  We  shall  compare  the  sizes  of  these 
three  driving  forces  by  expressing  them  in 
terms  of  |E|^.  Let  o  be  the  attenuation 
coefficient  of  microwave  intensity  as  defined 
by  the  eguation  I(x)  =  loe"®*.  Then  the  first 
term  is  poP  =  al  =  {ac'lt/8n)  \%\^  .  If  the 
dielectric  permittivity  is  linear  in  p,  then 
(de/3p)T  =  £/p,  so  eg.  (18)  gives  o  = 
(e  IE  1  ^  )  /  (24ii)  and  the  electrostrictive  term 
gives  [  (PECvKt  )  /  (24ll&p)  1  dc  |E|  ^  .  Denote  by  tr 
the  rise  time  of  a  microwave  pulse,  then 
dtlE|2  -  |E|2/tr.  Thus,  the  relative  mag¬ 
nitudes  of  the  electrostrictive  and  polariza¬ 
tion  terms  with  respect  to  the  thermoelastic 
term  are,  respectively. 


Electrostrictive  effect  _  pCvKtV£  1 
Thermoeiastic  effect  “  3caPp  ' tr ' 

(23) 

Polarization  effect  _  Ve . 1 
Thermoeiastic  effect  ~  co’tT' 


One  sees  that  the  second  and  the  third  driving 
forces  may  be  larger  than  the  thermoeiastic 
term  if  the  microwave  is  rapidly  modulated, 
contrary  to  the  estimate  by  previous  theory 
that  the  electrostrictive  effect  was  two 
orders  of  magnitude  smaller  (6) .  For  water  at 
30  *0,  p  =  1  gra/cm^,  C»  =  0.987-Cp  =  4.132x10’' 
erg/gm'C,  Kt  =  4.46x10"*'  cra^/dyne,  Pp  = 

2.8x10"^  "C-*,  and,  at  3  GHz,  a  =  0.88  cm"* 
and  e  =  80,  the  above  ratios  are  then,  re¬ 
spectively: 


Electrostrictive  effect  _  0.75  nanosec 

Thermoeiastic  effect  tr  ' 

(24) 

Polarization  effect  ^  0.34  nanosec 

Thermoeiastic  effect  ”  tr 


One  may  conclude  that,  for  water,  a  rise  time 
of  1  nanosecond  is  the  marginal  modulation 
rate,  beyond  which  the  electrostrictive  term 
and  the  polarization  term  contribute  more  to 
the  pressure  waves  than  the  thermoeiastic 
term. 

In  the  absence  of  the  electrostrictive  force, 
the  wave  eguation  for  the  pressure  waves 
reduces  to  the  one  developed  for  thermoacous¬ 
tic  process  (10,11).  While  we  have  not  yet 
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completed  calculations  on  the  pressure  waves 
generated  by  the  electrostrictive  effect, 
pressure  waves  from  the  thermoacoustic  effects 
have  been  reported  (10,11),  and  measurement  on 
the  amplitude  and  pulse  characteristics  have 
also  been  made  (11)  .  Contrary  to  previous 
theory  that  thermoacoustic  pressure  waves  were 
generated  mainly  by  inhomogeneous  distribution 
of  microwaves  in  dielectrics  (6),  it  was  found 
that  pressure  waves  are  generated  whenever 
there  is  a  discontinuity  in  thermal  or  dielec¬ 
tric  parameters  in  the  medium.  Eguations  (19) 
and  (20)  also  indicate  this  effect  for  elec¬ 
trostrictive  pressure  waves.  With  an  air- 
water  system  and  neglecting  the  electrostric¬ 
tive  force,  it  was  found  that  a  pealc  SAR  of 
4.5  )cW/gram  sufficed  to  generate  a  pressure 
wave  of  pealc  pressure  10*  pascal  (i.e.,  10^ 

dynes/cm^),  which  is  eguivalent  to  about  10% 
of  the  ambient  atmospheric  pressure.  This 
level  of  peak  pressure  may  be  considered  as 
the  upper  limit  of  linear  approximation.  As 
to  the  ratio  of  electrostrictive  effect  to 
thermoacoustic  effect,  it  depends  on  the  rise 
time  of  the  microwave  pulse,  being  greater 
than  one  for  rise  time  shorter  than  a  nanosec¬ 
ond,  and  smaller  than  one  for  rise  time  longer 
than  a  nanosecond  This  result  is  in  contrary 
to  previous  estimate  based  on  a  simpler  theo¬ 
ry,  which  gave  a  ratio  of  the  order  of  10"^, 
and  independent  of  the  rise  time  of  microwave 
pulses . 
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ABSTRACT 

A  formulation  for  scattering  of  electromagnetic  fields  by  arbitrary 
three-dimensional  dielectric  objects  is  presented.  The  formulation  also 
provides  an  inverse  solution,  viz.,  determination  of  the  dielectric  profile  of 
scatterers  through  measurement  of  the  scattered  fields  in  a  finite  region 
outside  of  the  target.  Uniqueness  of  the  inverse  scattering  solution  is 
assured  owing  to  the  dissipativlty  of  the  propagation  medium. 

I .  INTRODUCTION 

Analytic  expressions  of  scattering  of  electromagnetic  fields  may  be 
obtained  only  in  simple  cases,  such  as  the  Mie  scattering  from  a  uniform 
dielectric  sphere,^  or  scattering  from  concentric  dielectric  spheres,  which 
may  be  treated  as  extension  of  Mie  scattering.  For  arbitrary  dielectric 
objects,  numerical  approach  must  be  taken.  Even  so,  there  are  many 
theoretical  difficulties  which  require  certain  approximation  schemes.  The 
simplest  among  them  is  the  Born  approximation,  which  takes  the  first  order 
term  in  the  expansion  of  the  scattered  field  by  iteration.^  A  better  approach 
was  developed  by  Riclmiond  which  uses  the  moment  method, however,  it  applies 
to  only  two-dimensional  objects,  viz.,  objects  with  cylindrical  symmetry.  The 
moment  method  has  also  been  used  to  evaluate  the  field  inside  an  arbitrary 
three-dimensional  dielectric  body, but  it  has  been  found  that  the  numerical 
solutions  tend  to  diverge  with  respect  to  the  subdivision  of  cells. ^  Most 
recently,  a  testing  procedure  was  developed  to  analyze  the  numerical  stability 
and  reduce  the  computation  time.*  These  approaches  approximated  the  field 
inside  the  scatterer  by  a  linear  combination  of  some  basis  functions  and  then 
numerically  solve  the  linear  coefficients.  In  this  paper,  we  reformulate 
Richmond's  moment  method  in  an  operator  form  and  extend  it  to  three 
dimensions.  We  then  develops  a  method  to  calculate  the  integrals  of  the 
Green's  function  in  analytic  form,  thereby  provide  a  formula  for  the 
scattering  field  from  three-dimensional  scattering  objects.  The  formulation 
provides  a  solution  of  the  field  anywhere,  inside  or  outside  the  scatterer. 
Furthermore,  it  may  be  inverted  to  provide  an  inverse  scattering  formulation, 
from  which  one  may  develops  a  formula  to  obtain  the  three-dimensional 
dielectric  profile  from  the  measurement  of  the  scattered  field. 

II.  THE  SCATTERING  FORMULATION 

Consider  an  arbitrary  dielectric  body  placed  in  a  homogeneous  background 
medium.  Let  Em  and  fin  denote,  respectively,  the  dielectric  permittivity  and 
the  magnetic  permeability  of  the  background  medium,  and  let  E  (x)  and  ii(x) 
denote  those  inside  the  dielectric  object.  Thus  the  quantity  e (x) -Em  vanishes 
outside  the  scatterer.  We  shall  limit  our  scattering  problem  to  non-magnetic 
object  such  that  |i(x)  >  lim  for  all  x.  Consider  now  a  plane  wave  incident  upon 
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the  dielectric  object.  In  the  absence  of  free  charge  and  current,  the 
Maxwell's  equations  of  the  system  may  be  written  in  the  form:  M|f>  =  0,  where 
M  represents  the  Maxwell's  differential  operators  and  f  the  field  quantities. 
The  total  field  may  be  expressed  as  the  sum  of  the  scattered  field  and  the 
incident  field:  f  =  f*®*  +  f<^*.  We  also  separate  the  operator  M  into  two 
parts:  M  =  Hn-S,  with  Hn  representing  the  Maxwell's  operator  in  a  homogeneous 
background  and  -S  the  operator  due  to  excess  dielectric  permittivity,  E(x)-era. 
Then  the  Maxwell's  equations  may  be  written  as  Mni|f*^’>  +  Mra|f‘®’>  =  S|f>. 
Since  f*^’  satisfies  the  Maxwell's  equations  in  the  homogeneous  medium,  so 
Mm|f*^’>  =  0.  The  equation  then  becomes:  Minlf'“’>  =  Slf>.  Operating  both 
sides  by  the  inverse  of  Mm,  denoted  by  Gm,  one  then  gets: 

If'®' >  =  G„Slf>.  (1) 


It  is  remarked  that  the  Maxwell's  operator.  Mm,  in  general,  has  no  inverse 
because  it  has  a  null  space.  However,  we  are  only  concerned  with  the 
scattered  waves  that  contain  the  factor  expO’kmr),  which,  owing  to  a  positive 
imaginary  part  of  km,  vanishes  at  infinity  faster  than  any  power  of  r. 
Therefore,  in  the  subspace  of  scattered  waves.  Mm  has  an  inverse,  which  is  the 
Green's  function  Gm .  Here  we  also  remark  that,  owing  to  the  existence  of  the 
inverse  of  the  Maxwell's  operator  in  the  dissipative  medium,  the  inverse 
scattering  problem  will  have  a  unique  solution. 


Thus  the  scattering  problem  has  been  reduced  to  a  source-field  problem  in 
the  homogeneous  medium,  with  the  source  Slf>.  Adding  |f‘'’>  to  both  sides  of 
eq.  (1)  and  rearranging  terms,  the  equation  gives  (1  -  GmS) | f >  =  lf‘''>,  or, 

lf>  =  (1  -  GmS)-'  |f<">  (2) 

Eq.  (2)  gives  the  total  field  in  terms  of  the  incident  field.  Specific 
expressions  of  Gm  and  S  depend  on  the  choice  of  |  f>,  which  can  be  a  one- 
dimensional  potential  |<p>,  or  a  three-dimensional  electric  field  1E>,  or  a 
four-dimensional  potential  |A>  =  |ip,A>.  For  any  choice  of  |f>,  the  operators 
are  derivable  from  the  Maxwell's  equations.  Once  eq.  (1)  is  obtained,  one  may 
re-factorlze  the  product  GmS  to  any  pair  of  operators,  (Gm'.S'),  such  that 
Gm'S'  =  GmS.  So,  for  the  case  I  f>  being  the  electric  field  1E>,  we  first 

write  the  Maxwell's  equations  in  the  form:  Mm'|f'®'>  =  S'|f>,  where,  in  the 

Gaussian  system  of  electromagnetic  unit,  the  operators  Mm 'and  S'  are  as 
follows : 

Mm'  =  [V  +  km "3  •?  (3) 

S'  =  -  .L  (W  +  km^-y](e  -  Em)  (4) 

Cm 


where  a  time-dependency  of  has  been  assumed.  From  eq.  (3)  one  obtains 

the  x-representation  of  the  inverse  of  Mm': 


_  1  exp(ikm |x  -  x' I )  V 
47t  |x  -  x'  I 


(5) 


Since  [W-  +  km^  •?]  commutes  with  Mm',  it  commutes  with  its  inverse,  Gm ' . 

Therefore,  from  eq.  (4),  Gm's'  =  [W-  +  km^ -ijCm '  [- (£  -  Em)/Em].  Defining  Gm 
as  [W*  +  km^'i]Gm',  we  may  re-factorize  Gm'S'  to  the  product  of  the  following 
pair  of  (Gm,S): 
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Gm  (X,X') 


expdkm  |x  -  x'  I ) 
lx  -  x' 1 


(6) 


s  =  -  4-  [e(x)  -  e„] (7) 

bm 

The  operator  S  vanishes  outside  the  scatterer  and  Is  diagonal  In  the  sense 
that  S(x,x')  =  S(x)S(x-x').  In  the  next  section,  we  shall  convert  the  above 
Integral  equations  to  matrix  equations  by  digitizing  the  operators  Gm  and  S. 

III.  DIGITIZING  THE  SCATTERING  EQUATION 


To  calculate  the  inverse  of  (1  -  GmS)  in  eq.  (2) ,  we  divide  the  space 
into  small  cells  so  that  the  integral  operators  Gm  and  S  may  be  converted  Into 
matrix  operators.  We  first  divide  the  entire  space  into  subspaces,  Vo 
(a=0, 1, 2, . . . ) ,  with  Vo  being  the  region  occupied  by  the  scatterer.  We  further 
divide  each  of  Vo  into  N  small  cells,  labeled  by  tj,  with  j  =  1,  2,  ...,  N. 

Assuming  that  each  cell  Tj  is  small,  the  field  f*“’  and  the  dielectric 
permittivity  e  inside  each  cell  may  be  represented  by  their  respective  mean 
values  inside  the  cells.  This  means  that  the  dimension  of  the  cells  must  be 
smaller  than  a  fraction,  say  1/4,  of  the  electromagnetic  wavelength  inside  the 
scattering  object.  Denoting  by  Xj  the  center  of  the  cell  Xj,  eq.  (1)  may  then 
be  converted  into  summation  of  integrals  over  cells  Xj  as  below: 

<Xlf*”>  “  S  f  <XlGmlx'>dx'  S(X))<Xj|f>  (8) 

^  ^  x'e’tjcVo 

where  we  have  used  the  diagonality  of  the  operator  S;  the  integral  |dx'  is 
carried  over  only  the  space  Vo  since  S(x')  vanishes  for  all  x'  outside  of  Vo. 

Let  fa'®'  denote  a  N-dlmensional  vector  (of  which  each  element  is  a 
3-dimensional  electric  field)  representing  the  scattered  field  in  the  N  cells 
in  the  region  Va,  and  fo  and  fo ' "  ,  respectively,  the  total  field  and  the 
incident  field  in  the  N  cells  in  the  scattering  body  Vo  .  Also  denote  by  Ga 
and  S  the  N  x  N  square  matrices  with  the  matrix  elements  defined  as  below: 

(Ga)i]  =  J  G(Xi,  x')  dx'  with  XjEXjcVa  (9) 

x'exjcVo 

“1 

SjJ  =  —  le(Xi)  -  Em]-1  5j]  with  Xj,3EXj,)CVo  (10) 

Cm 

Then  the  eq.  (8)  may  be  abbreviated  as: 

|fa‘”>  =  GaS|fo>  (a=  0,1,2,  ...)  (11) 

Note  that,  for  Sij,  both  Xi  and  Xj  belong  to  Vo,  the  scattering  object, 
whereas  for  (Ga)i],  Xj  belongs  to  Vo  and  Xi  belongs  to  Vo,  which  may  be  Inside 
(if  a=0)  or  outside  (if  ou^O)  the  scatterer.  Eq.  (11)  represents  the  part  of 
eq.  (1)  for  the  subspace  Va .  For  a  =  0,  we  derive  the  equivalence  of  eq.  (2) 
from  eq.  (11),  as  we  did  in  deriving  eq.  (2)  from  eq.  (1) .  We  get: 

|fo>  =  (1  -  GoS)-' lfo'‘'  >  (12) 
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To  get  the  scattered  field  Inside  or  outside  the  scatterer,  we  substitut . 
eq.  (12)  into  eq.  (11).  With  some  simple  manipulation,  it  gives: 

=  G„(S-1  -  (a  =  0,  1,  2,  ...).  (13) 

The  a]30ve  equation  gives  the  scattered  field  in  any  region,  Vo,  in  terms  of 
the  incident  field  and  dielectric  permittivity  inside  the  scatterer.  In  the 
next  section,  we  shall  evaluate,  in  analytic  form,  the  matrix  elemenLs  of  Go 
and  Go  in  eq.  (13),  which  are  the  Integrals  Gi j  defined  in  eq.  (9). 

IV.  CALCULATION  OF  THE  GREEN'S  FUNCTION  INTEGRALS 

The  two  matrices.  Go  and  Go,  in  eq.  (13)  are  )tlnetic  quantities,  in  the 
sense  that  they  depend  only  on  the  bac)cground  medium.  The  matrix  elements  of 
these  two  matrices  are  defined  in  eq.  (9)  as  integrals  of  the  Green's  function 
over  the  cells  Inside  Vo,  and  with  Xj  being  the  center  of  the  cell  ti,  which 
is  inside  Vo  for  the  matrix  Go,  or  Inside  Vo  for  the  matrix  Go.  To  derive  the 
Integrals,  one  substitutes  eq.  (6)  for  the  integrand  of  eq.  (9),  the  latter 
then  gives: 


(Go ) j  ) 


e'ticV 


jV,- 


exp(i]Cm  IXj  -  at'  I ) 
Ixj  -  X' I 


dx' 


(14) 


where  Vj  represents  the  gradient  operator  with  respect  to  xj  .  Since  Xj  in  the 
integrand  always  comes  with  x'  in  the  form  cf  xj-x',  we  may  replace  Vj  with 
the  gradient  operator  with  respect  to  the  vector  Xj-x.'  Assuming  that  the 
volume  of  the  integral,  tj,  is  sufficiently  small  such  that  it  may  be 
represented  by  a  sphere  of  radius  a,  centered  at  Xj ,  then  the  above  integral 
depends  on  Xi -xj ,  a,  and  )Cm  only.  Denoting  |Xi-x'|  by  |x|,  the  above  equation 
then  gives: 


(Go) 1 ) 


exp ( ikm |x| ) 
lx| 


dx 


(15) 


where  the  origin  of  the  coordinates  is  xj ,  and  t  is  a  sphere  of  radiuc  a 
centered  at  xj  .  Since  (V^  +  km^  )  •  [exp(ikni  (x| )  /  |x|  ]  =  -4n6  (x)  ,  so 


(Go)ij 


exp(ikro |x| ) 
lx| 


dx 


+ 


5 (X) dx 


(16) 


The  second  integral  on  the  right  hand  side  gives  unity  when  the  sphere  t 
includes  the  origin,  i.e.,  when  «=j,  and  vanishes  otherwise.  As  to  the  first 
integrals,  it  may  be  evaluated  with  the  help  of  the  Identity,  Vg(r)  =  ±g'(r), 
and  the  generalized  Gauss's  theorem,  which  relates  a  volume  integral  to  a 
surface  Integral.  The  above  equation  then  gives: 


(Go)  i  j 


6ii 


d  ,exp(ikmr) 
’dr^  r 


The  evaluation  of  the  above  surface  integral  is  straightforward;  the  result 
depends  on  whether  the  sphere  of  Integration,  t,  contains  the  origin  (i.e., 
whether  i  =  j.)  Thus  we  have,  for  o  =  0,  1,  2,  ...; 
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i  =  j:  (Go),,  =  1  II 

i  *  i:  {Ga)n  =  (co3(Jcma)  -  sin  ()c„a)  ] 

I Xi  1  Km“ 


(17) 


I  S(1  - 


1  1 
ikn|xx-xj|  («k„  |Xi-Xj  I )  ^ 

*i-*l  ,1  _  3  _ 


IXj-Xjl  |Xj-Xj 


ikmIXj-Xjl  (ikmIXj-Xjl) 


r)  ]  (18) 


With  the  above  analytic  expressions  for  the  matrix  elements  lGa)i],  and 
the  expression  for  the  matrix  element  Si  j  given  by  eq.  (10),  one  may  then 
calculate  the  scattered  field  in  any  region,  say,  Va,  using  eq.  (13) . 


IV.  SUMMARY  AND  CONCLUSION 


We  have  developed  a  three-dimensional  scattering  formulation  from  which 
one  may  calculate  the  scattered  vector  field  anywhere,  inside  or  outside  the 
scattering  object;  the  scatterer  may  be  of  arbitrary  three-dimensional  shape 
with  arbitrary  dielectric  profile.  There  are  two  practical  problems  that  need 
further  investigation.  One  is  the  stability  of  the  inversion  of  the  matrix 
(S''-Go)  in  eq.  (13).  Another  is  the  size  of  the  matrices  Go  and  Ga,  viz., 
the  number  of  cells  into  which  the  scattering  object  is  divided.  We  shall 
discuss  these  two  problems  in  the  following  paragraphs. 

On  the  stability  of  matrix  inversion,  the  question  is:  How  stable  is  the 
inverse  of  (S' '-Go)  if  there  is  an  error  in  (S' '-Go)?  This  subject  may  be 
investigated  by  analyzing  the  so  called  condition  number  of  the  matrix 
(S' '-Go),  which  is  a  measure  of  singularity  of  the  matrix.®  The  singularity 
of  a  matrix  may  be  reduced  by  eliminating  symmetry  of  the  matrix,  viz., 
eliminating  the  interdependency  between  the  columns  or  the  rows.  Therefore, 
owing  to  the  term  S' ' ,  there  is  not  much  stability  problem  if  the  dielectric 
profile  of  the  scatterer  is  highly  asymmetric.  For  a  symmetric  dielectric 
body,  it  may  be  necessary  to  reduce  the  symmetry  by  other  means,  such  as 
making  the  division  into  N  cells  asymmetrically. 

As  to  the  size  of  the  matrices,  it  is  limited  by  two  conditions.  One  is 
the  desired  resolution  or  degree  of  accuracy,  and  another  is  the  approximation 
that  we  used  to  digitize  the  original  Integral  equation.  The  first  limitation 
is  a  subjective  one  which  needs  not  be  remarked.  As  to  the  second  limitation, 
recall  that  in  deriving  eq.  (8),  it  was  assumed  that  the  values  of  the  total 
field  and  the  dielectric  permittivity  Inside  each  cell  of  the  scattering 
object  do  not  vary  much.  For  the  field  to  remain  constant  inside  each  cell, 
the  cel'l  must  be  smaller  than  a  quarter  wavelength  inside  the  scattering 
dielectric.  In  this  aspect,  the  computation  time  increases  with  the  cube  of 
the  frequency  of  the  scattering  waves.  It  is  worth  noting  that,  as  some 
Investigators  have  found,  higher  symmetry  in  the  division  of  cells  may 
facilitate  the  computation  by  consideration  using  group  algebra;'"  this 
argument  applies  only  to  the  computation  involving  the  matr-'ces  Ga  since  they 
do  not  involve  the  dielectric  properties  of  the  scatterer,  wnlch  may  be  quite 
asymmetric.  On  the  other  hand,  as  we  mentioned  above,  higher  symmetry  may 
increase  the  condition  number  of  the  matrices  and  thereby  reduce  the  accuracy 
of  the  result.  For  scatterers  with  some  symmetric  properties,  it  is  possible, 
however,  to  utilize  both  arguments  to  benefit  the  computation.  One  may  start 
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by  using  a  higher  synmietrlc  division  of  cells  and  then  block-diagonalize  the 
matrices  Ga  and  (S'^-Go)  of  eq.  (13)  to  irreducible  forms.  The  computing  time 
is  reduced  by  the  diagonalizing  operation  while,  since  each  diagonal  block  is 
highly  asymmetric,  the  stability  of  the  data  may  be  preserved.  However, 
rigorous  proof  and  further  investigation  on  this  argument  are  needed. 

Finally  we  remark  that  the  formulation  described  in  this  paper  may  be 
inverted  to  provide  an  inverse  scattering  formula ^ ^  Briefly,  the  approach 
is  to  make  as  many  scatterings  with  independent  incident  waves  as  there  are 
number  of  cells  in  the  target.  For  each  scattering,  we  have  a  pair  of  fa '  “ ' 
and  fo ’ ^  ’  that  satisfy  eq.  (13).  Lining  up  the  column  vectors  representing 
fa'”  and  fo ' ^  ,  respectively,  we  obtain  two  N  x  N  square  matrices,  and 

F**>.  Then  eq.  (13)  gives;  Fa'”  =  Ga  (S' '  -  Go ) '  *  Fo  '  ”  .  The  equation  may 
then  be  inverted  to  yield  an  expression  of  S'*  in  terms  of  the  F'”’  . 
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Abst  ract 

The  complex  permittivities  of  three-O imens i ona i  inhomogeneous  Oiologicai  poOies  can  pe 
extracted  from  microwave  scattering  data  by  inverse  scattering  approach.  A  wa t e r -  i  nTne r s ed 
microwave  system  is  used  to  contract  the  wavelength  to  mill imeter  range  and  to  enhance 
impedance  matching  with  the  biological  body.  Contraction  of  the  wavelength  increases  the 
I  mage  resolut  lon.  wh lie  i mpedance  ma  t  ch i ng  promotes  the  m i c  r  owa ve  penetrat  lon.  Scattered 
fields  are  measured  using  an  array  of  127  o i po i e  elements  ana  a  total  size  of  a pp r o > i ma t e i ^ 
15cm  X  18cm  with  operating  frequency  at  3  GHz.  Two  inverse  scattering  approaches  have  been 
developed.  One  approach .  wh i ch  has  been  published  earlier.’-^  utilizes  an  inverse 
scattering  theorem  which  may  be  considered  as  a  generalization  of  the  Lorentz  reciprocity- 

theorem  to  dissipative  media.  The  other  approach,  which  is  presented  in  this  arric'e.  tafe-y 

scattering  measurement  by  an  array  with  various  directions  of  incident  wa-ve;  the  wave 
equation  is  converted  to  a  matrix  equation  by  dividing  tne  dielectric  body  into  a  number  o* 
cells,  the  dielectric  data  is  then  obtained  by  inverting  the  matrix  equation.  in  botn 

approaches,  uniqueness  is  assured  owing  to  the  d i ss i pa t i v i t y  of  the  propagation  medium, 

I n  t  r  oduc  t  I  on 

M  I  c  r  owa  ve  for  med  i  ca  i  i  mag  ing3-‘.&  provides  niany  ad  van  t  a  ges  over  i_  u  r  r  e  n  t  technique'.-.  ’  rie 
underlying  idea  here  is  to  reconstruct  a  t h r ee - d i mens i ona i  dielectric  image  based  on  the 
scattering  data.  The  technique  is  relatively  safe  since  only  non-ionizing  laaiation  is 

used .  Compa  r  ing  to  all  existing  techniques,  it  provides  the  best  contrast  be  twee  ri  s  o  f t 

tissues  in  terms  of  pa  t  ho  i  og  i  ca  i  information  since  it  is  the  dielectric  constant  t'-'dt  < 
imaged.  Simply  by  analyzing  tne  dielectric  image.  it  can  also  be  used  t  i.i  provide 

temperature  mapping  in  soft  tissues  with  a  resolution  of  one-tenth  of  a  degree  Ce  >  'us. 
Economical ly,  the  cost  of  a  microwave  imaging  system  is  uxuecteo  to  be  much  less  than  those 
using  the  existirig  imaging  techniques.  As  to  the  resolution.  it  depenqs  on  tne  nurriper  f 
elements  in  the  detection  array  and  the  image  reconstruction  method.  On  a  system  of  127 
array  elements  that  iS  being  developed,  we  have  achieved  a  transverse  reso  ution  -  ,3  u6 
width)  of  about  5.3  mm ,  which  is  about  haif-waveiength  in  water,  and  an  axiai  resc'iution  of 

about  11.7  mm,  which  is  about  one  wavelength  i  ri  watei  .  a  •  s  .  r 

in  order  to  reduce  reflection  of  the  incident  microwave  from  the  target,  boir 
tiansmitting  and  detecting  arrays  are  submerged  in  water.  Submerging  the  arrays  it  water 
not  only  provides  impedance  matching,  but  a i so  contracts  the  wavelength  and  tnereby  improves 
the  resolution.  Another  advantage  of  submerging  tne  arrays  in  water  is  e  i m i  iia t  i c n  of 

multipath  Signals  (viz..  Signals  scattered  from  surrounding  equipments  outside  of  the 
target.)  and  thereby  improving  s i gna i  - 1 o - no i se  ratio.  There  are  a  few  methods  to  further- 
improve  the  resolution,  all  of  them  require  aoditionai  scattering  data  collection.  However, 
increasing  the  amount  of  data  does  not  always  improve  tne  imaging  resoution  since, 
depending  on  the  reconstruction  methods,  redundancy  ana  i n t e r - depe nde nc v  of  reconstruction 
a  1 y o  r  I  i hm  may  negate  the  additional  inf  orma  t  i on  i npiu  t  .  The  pr  i ma  ry  concern  of  this  attic'e 
is  on  the  algorithm  of  image  reconstruction  from  the  meaSurertieot  s  of  scattered  fieid. 

The  objective  of  microwave  imagery  is  to  determine  the  dielectric  characteristics  of  a 
target  from  the  microwave  field  pattern  exterior  of  the  target.  For  far-fieid  and 
intermediate  zones,  the  traditional  approaches  are  the  Fraunhofer  and  Frenei  diffractions 
based  on  the  theory  of  Fourier  t r ans i orma t i on .  in  many  applications,  particularly  in 
medical  imagery,  the  target  must  be  in  the  vicinity  of  the  microwave  detector,  then  these 
approaches  are  not  applicable  since  neither  the  field  pattern  of  the  target  nor  the  field 
pattern  of  the  receiving  antenna  may  be  characterized  by  a  Fourier  transform  of  the 
structure  of  the  target  or  the  antenna.  Though  this  aspect  c omp i  icates  the  ma themat  ' ca i 
problem  of  near-field  imageries,  development  of  some  near-fieid  approach  may  provide  a 
pcissibi I ity  of  three-dimensional  focusing  wh i cn  is  not  achievable  in  either  Fraunhofer  or 
Frenei  diffractions. 

An  important  theoretical  problem  in  imageries  concerns  the  uniqueness  of  the  target 
reconstruction.  The  dielectric  structure  of  the  target  is  part  of  the  source  of  the  field 
exterior  to  the  target.  in  particular,  apart  from  undesirable  noises,  it  is  the  sole  source 
of  the  scaitcred  field.  The  Maxwell  s  equations  guarantee  a  unique  solution  of  the  field 
once  1  he  source  and  the  boundary  condition  are  specified.  However,  the  reverse  is  not  true. 
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indeed,  a  simple  consideration  ot  tne  field  due  to  charges  on  a  conducting  sphere  snows  tnat 
the  same  field  exterior  to  the  sphere  may  De  given  Dy  many  possible  distributions  within  the 
sphere,  for  example,  a  distribution  of  surface  charges  on  a  sphere  of  sma l ler  radius,  or  a 
point  charge  at  the  center  of  the  sphere.  Therefore,  theoretically,  knowing  the  field 
everywhere  extei lOr  to  an  imaging  target  does  not  guarantee  a  unique  solution  of  the 
dielectric  structure  within  the  target.  In  this  paper  we  shall  show  that  if  the  background 
medium  is  dissipative,  then  the  inverse  scattering  has  a  unique  solution.  One  may  then 
argue  that,  for  theoretical  calculation,  even  in  a  non-d i ss i pat i ve  medium,  one  can  add  a 
small  dissipative  part  in  the  dielectric  constant  to  start  with  and  later  set  it  to  zero. 
The  problem  is  that,  even  though  the  solution  is  theoretical ly  unique,  to  reconstruct  the 
target  requires  a  large  number  of  independent  data.  The  sma I ler  is  the  imaginary  part  of 
the  dielectric  constant,  the  more  redundancy  in  the  measured  data.  On  the  other  hand,  the 
stronger  the  dissipation  of  the  background  medium  is,  the  weaker  the  scattered  signal  is, 
and  therefore  the  smaller  is  the  s i gna i -to-no i se  ratio.  Thus,  we  are  faced  with  a  practcai 
problem  of  optimizing  the  d i ss i pat  i  v i ty  of  the  background  medium. 


Another  practical  aspect  further  complicates  the  issue:  in  reality,  one  can  measure  the 
field  only  at  a  limited  number  of  points  exterior  to  the  target.  The  problem  is  a  typical 
inverse  scattering  problem,  viz.,  how  much  information  on  the  dielectric  structure  of  the 
target  can  one  infer  fr om  a  set  of  limited  data  of  the  field  exterior  to  the  target^ 
Previously,  we  developed  an  inverse  scattering  theorem’ ■ »  and  a  method  of  phase  and 
amplitude  conjugation  techniques^ ■  to  reconstruct  tne  target  dielectric  profile,  with  which 


a  transverse  resolution  of  \/2  and  an  axial  resolution  of  x  were  obtained,  i 
we  present  another  method  which  is  based  on  discretely  samp i i ng  the  Green's 
then  converting  it  to  a  matrix  (referred  to  as  the  Green's  matrix).  With  this 
resolution  can  be  improved  unlimitedly;  the  resolution  i mp roves  linearly  with 
sanipting  points  of  the  scattered  field  anywhere  outside  of  the  target. 


e ,  with  wh I ch 
n  this  paper  , 
f  unct I  on  and 
approach,  the 
the  number  of 


The  scattering  formulation 

in  this  section  we  describe  tne  logic  of  the  formulation  and  present  the  scattering 
formula.  First,  we  mention  a  few  words  about  the  notation.  Since  we  are  dealing  with 
scattering  of  electromagnetic  fields  Dy  three-dimensional  objects,  spatial  coordinates  are 
t h r ee - d i mens lonai  vectors  wh iie  field  quantities  may  be  three-di mens lonal,  such  as  the 
electric  field,  or  f our -d imens i ona i ,  such  as  the  charge-current  density.  To  distinguish 
three-dimensional  quantities  from  f our -d imens i ona I  ones,  we  denote  the  former  with  bold  face 
letters  while  use  normal  letters  for  the  latter.  we  shall  also  use  Dirac  bra-ket.  say.  if,-, 
to  denote  a  field  quantity  f,  and  <xif>  to  denote  the  function  f(K),  which  is  f  in  the  x- 
representat i on .  Thus.  an  operator,  say.  the  Green's  function  G.  takes  the  form 
< X  I G I  X '  > =G( X , X '  ;  in  the  x-  representation.  unless  otherwise  stated,  all  products  between 
vectors  and  between  operators  and  vectors  are  vector  products,  namely,  the  products  are  to 
sum  over  either  the  3  or  the  4  components  without  the  t  notation;  we  reserve  the  E  notation 
for  summa  tion  over  the  numbe  r  of  ceils  in  discrete  samp  ling  of  data.  On  the  eiectr  omagne  t i c 
equat lOhs ,  we  shall  use  the  Gauss i an  system  of  eiectr  omagne  tie  unit:  conve  rsion  of  equatiuns 
between  the  Gaussian  system  and  the  mkS  system  may  be  found  at  the  end  of  the  text  book  by 
Jack  son . * 


Cons ider  an  arbitrary  dielectric  body,  wh ich  will  be  our  scattering  object,  occupying  a 
region  V,,  in  a  homogeneous  background  medium.  Let  and  Vm  denote  respectively  the 
dielectric  pe  rm ittivity  and  the  ma  gne  t i c  permea  b iiity  of  the  background  med i um .  and  let  <  <  x ) 
and  uCx)  denote  those  inside  V„ .  Thus  the  quantity  t(x)-(~.  ■ s  lora'ized  'n  the  sense  that 
It  vanishes  outs i de  Vo •  We  shall  limit  our  system  to  non-magne t ic  object  such  that  nx)  = 
Um  for  all  X  inside  Vo • 


Consider  now  a  p 
current,  the  Maxwe I  I 


lane  wave  incident  upon  the  target.  in  the  absence  of  free  charge  and 
s  equations  ot  the  system  may  be  written  in  the  fol lowing  form: 


where  M ,  wh ich  is  a  function  of  <<x),  represents  the  Maxwe ll's  differential  operators  and  f 
the  field  quantities.  The  total  field  may  be  expressed  as  the  sum  of  the  scattered  field 
and  the  incident  field:  f  =  f‘"’  +  fx'‘.  in  order  to  express  the  scattered  field  in  terms 
of  the  incident  field  and  the  dielectric  distribution  of  the  target,  we  separate  the 
operator  M  into  two  parts:  VW,  and  -S.  The  first  part,  IvUi,  represents  the  Maxwell's  operator 
in  a  homogeneous  background  medium  and  is  independent  of  (Cx),  The  second  part. 


represents  the  operator  due  to  excess  dielectric  permittivity,  «(x)-<m.  Then  eq .  (i)  may  be 
written  as  NW,  I  f  *  ‘  >  -*•  NUilf‘*’>  *  Slf>.  Since  f‘'»  satisfies  the  Maxwell's  equations  in  the 
homogeneous  medium,  so  lAn  I  f  ‘  >  *  0.  Eq,  (1)  then  reduces  to 
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Hr,  •  t  ‘  ^  > 


S  »  f  > 


<•  3  ) 


Thus  the  scattering  problem  has  been  reduced  to  a  source-field  problem  in  the  homogeneous 
medium,  where  the  scattered  field  corresponds  to  the  field  due  to  a  source  distribution 

^lf>.  However,  the  quantity  if>  on  the  right  hand  side  st i i i  contains  the  unknown  field 

To  r emove  the  scattering  fieid  tr om  the  right  hand  side  of  eq.  O).  we  operate  both 
sides  of  the  equation  by  the  inverse  of  NV«.  Note  that,  since  the  operator  fAr^  has  a  nu  i  i 
space,  it  in  general  has  no  inverse.  However,  since  we  are  only  concerned  with  the  outgoing 
fields  that  contain  the  factor  exp(\k#r»r>,  wntch,  owing  to  the  imaginary  part  of  krr^.  vanish 
at  infinity  faster  than  any  power  of  r.  we  may  then  limit  ourselves  to  the  subspace  of  such 
fields.  In  such  subspace.  NV«  has  an  inverse,  which  is  the  Green  5  function  Grr* .  uince 
IS  an  outgoing  scattered  wave,  so  GrtM«if‘**>  =  Then  eq.  ^3^  gives; 

I  f  *  -  >  >  =  Gr^S  If'.  ^  > 

Adding  if‘'»>  to  both  sides,  it  then  becomes  1  f  >  =  Gr«Si  f>  or, 

(  1  -  ■)if>  =  lf‘'>'  cbi 

it  can  be  proven  that  the  operator  (i-Grr>S>  or^  the  left  hand  side  has  an  inverse.  Therefore 

the  above  equation  gives 

1  f  >  =  Cl-  Grr^S  )  -  M  f  ‘  t  t>  > 

Subtract  the  incident  field  from  both  sides  of  the  above  equation,  one  then  obtains 

the  scattered  field 

I  f  *  -  >  >  -  t  (  1  -  Gr«S  )  '  ’  -  1  )  I  f  ‘  .  (  7  ) 

Eq .  (7>  gives  the  scattered  field  in  terms  of  the  incident  f>eid.  in  the  next  section  we 

shall  develop  a  technique  to  evaluate  the  inverse  of  (i-GmS>. 

So  far  we  have  described  the  scatter«ng  formulation  in  general  terms.  Specific 

expressions  of  these  equations  depend  on  the  cno'ce  of  the  quantity  if-;  if'  can  be  the  one- 

dtmensionai  scalar  potential  >«>.  or  the  three-d imens • ona i  electric  field  iE>.  or  the  four- 

dimensional  potential  IA>  *  For  any  choice  of  the  fieio  if',  the  operators  are 

derived  from  the  Maxwell's  equations.  in  be  1 ow  we  list  the  results  for  the  case  where  if' 
IS  IA>,  and  tor  the  case  where  if>  is  iE>.  Details  of  ttie  derivation  witi  pe  given  in  a 
forthcoming  paper. 

We  shall  assume  a  time-dependency  of  e-'»-«-  and  replace  d/at  with  -\w.  if  the  fieid  if' 
is  the  4 -d imens 1 ona )  vector  potential,  we  snail  use  the  following  notation  to  represent 
matrix  operation: 

H  _ 

I 


where  x  denotes  the  1 -d i mens i ona I  time  component  of  a  a - d i mens i ona I  vector,  and  (— — - ) 

denotes  the  3-d imens i ona l  space  component.  Then  tne  operators  M.,  and  S  may  be  represented 
as  below: 


r~  1 

-  —  w  t  < 

f  ^  rr» 

c 

IWUmf  <  -  I 


;  t  -frr,) 

.-(rr.)  T 


Note  also  that  the  operator  w  on  the  left  of  inside  the  matrix  not  only  operates  on 

hut  also  operates  on  whatever  quantity  to  be  operated  by  S.  Since  NW,  is  a 
differential  operator.  its  inverse,  G..,.  is  an  integral  operator.  In  the  X- r  ep  i"  ese  nt  a  t  i  on  , 
Gm  may  be  explicit  y  expressed  as  below: 


GwfK.K' ) 


<Xt  Gm I K'  > 


J _  exp(  1  kr.,1  X  -  X  ' 

4n  I  X  -  X  ■  I 


The  above  expression  means  that,  in  the  x-representat i on ,  Gm  is  an  integral  operator: 


<xlG,.,l  t  > 


/  G..,f  X  ,  X' )  f  (  X  '  )  dx  ' 
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Note  that  the  aet initions  ot  operators  G„  ana  S  are  not  unique  in  the  sense  that,  once 
eq .  <■  4 )  is  obtaineq,  one  may  re-factorize  the  product  G,wS  to  any  pair  ot  operators, 
CGm'.S),  such  that  G„'S'  =  G^S .  This  is  exactly  what  we  shall  Oo  for  the  case  where  the 
field  if>  is  the  electric  field  IE>.  We  shall  first  derive  a  pair  denoted  by  Gm  and  S', 
where  G„, '  corresponds  to  the  Maxwell's  operator  MW,',  then  re-factorize  the  pair  to  another 
pair  Gw,  and  S,  The  re- f actor i zat i on  will  not  affect  the  results  of  scattered  fields  or 
inverse  scattering  image  reconstruction.  but  simplifies  the  scattering  and  inverse 
scattering  formulation.  Thus,  for  the  case  where  if>  is  the  electric  field,  IE>,  we  first 
derive  from  the  Maxwell's  equations  lAw'and  S': 


MW,'  =  ♦  l(w,“  1  •  1 


I  VW  •  •»  Kw,^  •  I  1  <  <  -  <w, ) 


From  eq.  C11)  we  obtain  the  x- represent a t i on  of  the  inverse  of  MW,'; 


GW,'  (  K.  X  '  ) 


exp(  I  k.T'ix  -  i 
IX  -  X ' I 


Since  tww-  f  ■  T  1  i  n  eq  .  ( 

inverse,  Gw,'.  Therefore, 

IW-  *  Kw,“  •  T  1  Gw,  ■  1  -  f  f  -  <w,>/<w,i 

-T  ■  (  <  -  <w.)/<w,,  respectively; 

Gw,<  X  ,  X  ■  >  = 


(1^)  comnutes  with  MW,'  =  1  +  km*  1  •  T  ,  so  it  coiTrnutes  with  its 

,  from  eq.  <12),  Gw,'5'  =  -Gw.'lWW'-*'  kw,*  ■  T  1  (  <  -  tmJ/tw, 

tw,i.  Then  we  define  the  pair  G,„  and  S  as  I  wy  •  -i-  k„,*-llGw,'  and 


I  T  expc  1  kw,  I  X  -  X  '  I  ) 


lx  -  X '  1 


-  (  <  (  X)  -  <rr,)  •  I 

<fT> 


The  operator  S  in  eq .  t  1  tj  )  is  diagonal  in  the  sense  that,  in  the  x- r ep r esent a t i on , 
Slx.x'l  «  SCx)4(x-x  ).  It  also  vanishes  outside  of  the  target  region  .  It  is  the 
diagonal  elements  of  this  operator  that  we  wish  to  recover  from  the  scattering  data.  in  the 
next  section,  we  shall  covert  the  above  integral  equations  to  matrix  equations  by  digitizing 
the  operators  Gw,  and  S. 


In  this  section,  we  describe  a  numerical  method  to  calculate  the  inverse  of  -  Gw,5 ) . 
and  thus  the  scattered  field  f‘»’  using  eq .  (7),  for  the  case  when  f'»'  is  the  electric 
field,  E‘“’.  in  this  case.  Gw,  and  S  are  given  by  eqs .  <14)  and  (15).  Our  approach  will  be 
to  divide  the  space  into  small  cells  so  that  the  integral  operator  Grr,  may  be  converted  to  a 
matrix  operator.  To  this  end,  we  divide  the  entire  space  into  subspaces,  and  label  each 
subspace  by  v,, ,  with  being  the  entire  region  occupied  by  the  target.  We  further  divide 

each  of  into  N  sma ll  cells,  labeled  by  t,.  with  j  =  1,  2 .  N ■  We  assume  that  each 

cell  r  j  IS  sma  M  enough  that  the  field  f**’  and  the  dielectric  permittivity  r  inside  each 
cell  may  be  represented  by  their  respect! ve  mean  values  inside  the  cell.  This  means  that 
the  dimension  of  the  cell  must  be  smaller  than  a  fraction,  say  1/4.  of  the  microwave 
wavelength  inside  the  target.  We  shall  denote  by  Xj  the  center  of  the  cell  r^.  Then 
eq .  1 4 )  may  be  converted  into  summa  tion  of  integrals  over  ceils  tj  as  Pel ow ; 


<x I f ‘ ‘ > 


1  1  <XI(jw,IK'><X'ISIX" 

XX*’ »  f  >aM '  ax” 

=  j  J  <Xlljw,IX'>'S(X')6(X'-X"><X"l 

f  > dx ' dx ” 

1  <Xl<3w,IX'  >  ■  S(X'  )<x 

■  1  f  >ax'  =  z 

I  r  < X 1  (jw,  1  X  '  > 5(  X  ' )  < X  '  1  f  >ax '  1 

^  fV„ 

t 

L  X ' <T , c Vo  J 

1  1  <xi(3w,ix'  >dx' 

^  1  X  '  <  T  J  C  Vo 

j  S(X  J  XX  J  1  f  > 

(16) 

where  the  integral  /dx'  is  carried  over  only  the  space  Vo  since  Sex')  vanishes  for  all  x' 
outside  of  Vo. 

Denote  by  G, j  and  S, j  .  respectively,  the  following  quantities: 

G,  ,  «  I  Gw.(x»,  X')0X'  (17) 

X ' <  T  j  C  Vo 

S,  j  »  =  -  —  IcCX,)  -  Cw,!  6,J  (18) 
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By  approximating  the  integrated  region  tj  with  a  sphere,  the  integral  j  of  eg.  <i7)  may  Pe 
calculated  analytically;  the  analytic  form  is  presented  at  the  end  of  this  section.  With 
the  definition  of  Gv j  and  Svj.  eq.  <t6)  may  De  written  as: 


In  the  matrix  form,  the  above  equation  is: 


Gil  Gia  Gi3 

Ga 1  Ga  a  Ga  3 
G3 1  Gaa  Gas 


Gm 1  G^a 

Gm^  1  ,  1 


Gin  Gi.n-»i 
OaN  Ga  .  1 

Gsn  G3  .  N-*  1 


Gnn  Gn  .  N-*  1 
Gn>«'  1  .  N 


Si  0  0 

O  Sa  0 
0  0  Sa 


1 


where  subscripts  1  to  N  denote  cells  inside  Vo .  subscripts  N+1  to  2N  denote  those  inside  v,. 
etc.  We  may  break  the  above  equation  into  equations  representing  the  operation  for  each 
subspace  V*.  Thus,  for  the  space  Vo  occupied  by  the  target,  the  above  equation  gives; 


G11  Gia  Gta 

Ga 1  Ga  a  Ga  3 

Gai  Gaa  Gas 


Gn 1  Gna  • 


Gin  Si  0  0 

GaN  0  Sa  G 

G3N  0  0  S3 


1 

1 1  f  ^  I 


and.  for  the  subspaces  ca  *  1,  2.  eq.  (20)  gives: 


Goin*  1  .  1  G««n«  1  .  a  OotN*  1  .  3 
GoiN^a  .  1  G«,N*a  .  a  GoiN<»a  .  3 
Go«n*-3  .  1  Go,n*3  .  a  Go«n-3 . 3 


Go<N-kN  .  1  Go»N-»N  .  a  G«n*-n  .  3 


1  .  » 

OoiN-*-a  .  * 
0«(N«  3  .  > 


GwN-*-  N  .  I 


Si  0  0 

0  Sa  0 

0  0  $3 


The  aDove  two  equations  may  t)e  aooreviated  as: 


where  and  lfo>  represent,  respectively.  the  scattered  and  total  fields  in  Vo,  the 
region  occupied  by  the  target:  ifo,‘»»>  and  I f„>  represent  those  in  v*  («  =  1,  2,  ...>,  a 
subspace  of  any  volume  outside  of  the  target.  In  the  above  two  equations,  the  matrix 
elements  of  G„ ,  G,. ,  and  S  are  as  defined  in  eqs .  tl7>  and  (10)  with  the  following 
assignments  for  x.  and  x.: 


(Go) . j 


I  Gcx. .  X' : 


with  X.  C  T  ,  cVc 


(Go.)  .  , 


J  G(x. .  X' : 


with  X. «  T . cVo 


-  [<(X,)  -  6.,  with  Xx.jeT..jCVo 

<f*1 


The  values  of  the  matrix  elements  S, ,  in  eq .  C27)  are  the  quantities  to  be  reconstructed 
from  the  scattering  data;  the  values  of  the  matrix  elements  (Go>,j  and  (Go,),,  may  be 
calculated  from  the  analytic  forms  presented  at  the  end  of  this  section.  Note  that  eq .  (23) 
is  similar  to  eq .  (4).  Indeed,  it  represents  the  part  of  eq.  (4)  for  the  subspace  Vo,  the 
region  occupied  by  the  scattering  target.  We  may  then  derive  the  equivalence  of  eqs.  (6) 
and  (7)  from  eq .  (23),  as  we  did  in  deriving  eqs.  (6)  and  (7)  from  eq.  (4).  We  then  get: 
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lfo>  =  Cl  -  GoS)  -  ’  I  t  >  c28> 

=  (Cl  -  G„S5  ’  (29) 

To  get  the  scattereci  field  outside  the  target,  we  substitute  eq.  c28)  into  eq .  (24).  With 

some  simple  manipulation,  it  gives: 

lfo.'»‘>  *  G„CS-’  -  Go)"’lf„‘'>>  Ca  =  1.2....).  (30) 


The  above  equation  gives  the  scattered  field  in  any  region,  V„ ,  outside  of  the  scattering 
target  in  terms  of  the  incident  field  and  dielectric  permittivity  inside  the  target.  in  the 
next  section,  we  shall  invert  this  equation  to  obtain  the  dielectric  permittivity,  S,  inside 
the  target  in  terms  of  the  scattered  field  outside  the  target. 


We  shall  summarize  the  scattering  approach  described  in  this  section.  One  first  divides 
the  targe t  i n t c  a  number  of  cells,  say.  N  cells.  Then  one  calculates  the  incident  field  at 
the  center  of  each  cell  and  constructs  a  N-di mens i ona l  vector,  fo‘'’  cof  wh i ch  each  el emen t 
may  still  be  a  3-d i mens i ona l  vector,  such  as  an  electric  field).  One  also  constructs  a 
N  X  N  diagonal  matrix.  S,  of  which  each  diagonal  element  is  given  by  the  average  value  of 
the  dielectric  permittivity,  «.  in  eacn  cell  of  the  target.  Then  one  calctilates  integrals 
tGo)>j  =  XGt X, , X' ) dK' .  with  X,  being  the  center  of  the  ceil  and  integrating  x'  over  the 
cell  Tj;  both  T,  and  rj  are  inside  the  target.  Thus  one  constructs  the  N  x  N  square  matrix. 
Go,  with  matrix  elements  being  the  integrals  (Go)»j.  Finally,  to  obtain  the  scattered  field 
anywhere  outside  the  target,  say,  in  a  region  Vo.,  one  divides  the  region  into  N  cel  is  and 
calculates  the  integrals  <Go.).j  =  XG< . X ' >dx ' .  with  x,  being  the  center  of  the  cell  r,  in 
Vo,  and  integrating  x'  over  the  cell  Tj  i ns i de  the  target  region,  Vo .  One  thus  constructs 
the  N  X  N  square  matrix.  Go,,  with  matrix  elements  being  the  integrals  (Go.>,j.  The  scattered 
fields  inside  Vo.  can  then  be  obtained  from  eq .  (30).  with  the  ith  element  of  the  vector 
fo,‘">  representing  the  scattered  field  m  the  cell  t,  of  Vo,. 


Finally,  we  present  results  of  the  integrals  <Go),j  and  (G„),j 
consider  the  case  where  the  field  of  interest  is  the  electric  field, 
by  eq.  (14).  The  expressions  for  the  integrals  in  eqs .  (25)  and  (26) 


in  analytic  forms.  We 
then  G(x.x')  is  given 
are  then; 


i  =  j  :  (Go  )  ,  )  or  (Go.)  ,  , 


1  *  J  ■■  (Go)  ,  j  or  (Go)  ,  j 


I  I  1 


2 

3 


q  ^exp( ' kma) 
da  a 


(31) 


a 


exp('ko.ix>-Xi  I  ) 
I  X, -X , I 


•  t  cos  (  l(o,a  > 


1 

s I n( K„a ) ) 

'  iT>a 


•  I  1(1 


_ 1 _ 

1  kroIXx  -X,  I 


_ 1 

(  1  Ko,  I  Xx  -X  ,  I  )  “ 


) 


Xx -Xj  Xx -X, 

I  Xx  -X  J  I  I Xx  -X J  I 


1 _ 

1  Ko,  I  Xx  -X  ,  I 


( 1 k„IXx -Xj I )= 


(  32  ) 


where  we  have  assumed  that  each  cell  is  of  equal  volume,  and  that  the  volume  of  the  ceils 
are  small  and  may  be  approximated  by  a  sphere  of  radius  a.  Mote  that  the  expressions  for 
( Go )xj  and  for  ( Go,  >xj  are  of  the  same  form;  the  oniy  difference  is  the  value  of  Xx:  for  Go. 
Xx  IS  inside  the  scattering  target,  whereas  for  G* .  Xx  is  outside  the  target:  for  both,  xi 
IS  inside  the  target.  With  the  above  expressions  for  (G„)xj  and  (Go.)xj,  eqs.  (29)  and  (30) 
provide  a  formula  for  the  scattered  field  anywhere,  inside  or  outside  of  the  scattering 
target.  Previously,  complete  formula  for  scattered  fields  has  been  developed  only  for 
two-dimensional  targets. Formula  for  evaluating  fields  Inside  arbitrary 
three-dimensional  dielectric  bodies  has  been  derived,’’ ■’ »  but  it  has  been  found  that  the 
numerical  solutions  tend  to  diverge  with  respect  to  the  subdivision  of  cells.’)*  Most 
recently,  a  testing  procedure  was  developed  to  analyze  the  numerical  stability  and  reduce 
the  computation  time.’-*  However,  all  these  approaches  approximated  the  field  inside  the 
scatterer  by  a  linear  ciomb i na t i on  of  some  basis  functions  and  then  nume rically  solved  the 
linear  coefficients.  This  is  the  first  ti me  a  comp i e t e  f o rmu la  is  developed  for  the  vector 
scattered  fields  from  an  arbitrary  three-dimensional  target. 

Inverse  scattering  -  Green's  function  approach 

The  scattering  formulation  presented  in  the  last  two  sections  provides  a  method  to 
calculate  the  scattered  field  of  an  arbitrary  target.  In  this  section.  we  shall  show  that 
the  result  of  the  formulation  may  be  inverted  to  provide  a  method  to  evaluate  the  dielectric 
permittivity  distribution  of  an  arbitrary  target. 
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Cons i der 
direction  of 
accord i ng I y . 
i ndependent  . 


independent  vectors.  lfoi‘‘' 
wave  vector  k  j  ( j  *  1 .  2 , 

,  the  incident  matrix: 


eq.  (30).  If  the  incident  field  fo'"’  'S  changed,  say,  Oy  changing  the 

the  incident  wave.  then  the  resulting  scattered  field  wMi  oe  changed 
Two  incident  fields  with  different  incident  directions  are  linearly 
So,  if  the  target  is  divided  into  N  cells,  we  may  produce  N  linearly 

I  > . fow' ‘ ’ t .  each  of  which  corresponds  to  an  incident 

N).  and  thereby  construct  a  non-singular  N  x  N  square  matrix 


fo..  f. 


So.  if  the  incident  waves 
the  incident  matrix  Fo‘'’ 


>  are  the  plane  waves  exp<ikj-x).  then  the  matrix  elements  of 


exp < \ h  j 


the  center  of  the  ith  cell. 


The  dielectric  profile  of  the  target  will  be  obtained 
scattered  fields  by  inverting  eq.  (30).  The  measurements 
made  in  any  region  outside  the  target,  say,  in  the  region  v, 
defined  at  the  beginning  of  the  last  section.  Since  Vo 
divide  V  into  N  cells  and  measure  the  scattered  field  at  the 
incident  field.  Denoting  by  fj**’  the  N-d imens i ona i  vect 


fields 
F«-»  : 


in  the  N  cel  Is  for  the 


srget  will  be  obtained  from  the  measurements  of  the 
.  (30).  The  measurements  of  the  scattered  fields  may  be 

get.  say,  in  the  region  V,  which  may  be  one  of  the  v*. ' s 
last  section.  Since  Vo  is  divided  into  N  cel  Is.  we  also 
the  scattered  field  at  the  center  of  each  cell  for  each 
•>  the  N-d imens 1 ona I  vector  that  represents  the  scattered 
i  nc i dent  wave .  we  then  construct  a  N  x  N  scattered  ma  t  r i x . 


fwi  fwa  fw3  f  MM  J 

With  the  N  incident  waves  and  the  N  scattered  waves  lined  up  as  two  N  x  N  square  matrices, 
eq .  ( 30  )  gives: 

F‘->  =  G(S-’  -  G„)-’Fo‘'’  '36) 

where  we  have  dropped  the  subscript  <x  in  the  matrix  G  since  no  reference  was  made  to  any 
particular  Vo.  for  sampling  the  scattered  field.  In  the  above  equation,  each  column  of  the 
matrix  F‘«'  represents  the  scattered  wave  due  to  an  incident  wave  represented  by  the 
corresponding  column  in  F„‘ ' ’ ,  Each  element  in  a  column  of  F'»>  represents  the  scattered 
field  at  a  cell  in  V,  whereas  each  element  in  a  column  of  Fo''*  represents  the  incident 
field  of  the  corresponding  cell  inside  the  target  Vo-  Since  Fo‘'>  is  a  non-singuiar  matrix, 
we  may  invert  it  and  bring  eq.  (36)  to  the  following  form; 


which  then  gives 


Eq.  (38)  is  the  final  result  of  this  inverse  scattering  approach;  the  inverse  of  each 
diagonal  element  of  the  matrix  on  the  left  hand  side  gives  the  negative  relative  excess 
dielectric  permittivity,  -(<-<„)/<„,  in  the  corresponding  cell  of  the  target. 


The  inverse  scattering  approach  may  be  summarized  below.  One  first  divides  the  target 
into  a  number  of  cells,  say.  N  cells.  Then  one  produces  N  incident  fields,  all  of  which 
have  the  same  amplitude  and  frequency,  but  each  has  a  different  incident  direction.  For 
each  incident  field,  one  calculates  the  incident  field  at  the  center  of  each  cell  inside  the 
target  region  and  constructs  a  N-d imens iona I  vector,  fo' ' ’  (of  which  each  element  is  a  3- 
dimensional  electric  field  vector).  Lining  up  all  incident  vectors  for  the  incident  waves 
as  column  vectors,  one  then  forms  a  N  x  N  square  matrix,  F„''*,  referred  to  as  the  incident 
matrix.  Also,  for  each  incident  field,  one  measures  the  scattered  field  in  each  cell  in  any 
region  V.  which  is  outside  the  target;  the  measurement  region,  v.  and  its  division  into  N 


36  /  SPie  Vol  767  Medics!  Imaging  (1987) 


7 


cells  must  De  unchanged  for  aM  incident  fields.  Tne  scattered  ma  t  r i x ,  f ‘ ,  may  then  De 
constructed  from  the  measurement  ot  the  scattered  fields,  in  the  same  way  as  the  incident 
matrix  Fo^'*  was  constructed  from  the  incident  fields.  Then  one  calculates  the  integrals 
cGo)xj  =  JGC Xv , X ' ) dx '  given  by  eqs.  <3i)  and  <32>.  with  Xt  being  the  center  of  the  cell  t, 
inside  the  target  and  integrating  x'  over  the  cell  Tj.  also  inside  the  target.  Thus  one 
constructs  the  N  x  N  square  matrix.  Go,  referred  to  as  the  target  Green's  matrix,  with 
matrix  elements  being  the  integrals  <Go>xi.  Similarly,  one  calculates  the  integrals  (G)xj  = 
/GC X, . X ' ) dx '  from  eqs.  (31)  and  (32).  with  Xv  being  the  center  of  the  ce M  in  the 

measurement  region,  V,  and  integrating  x'  over  the  cell  Tj  inside  the  target  region.  Vo • 
One  thus  constructs  the  N  x  N  square  matrix,  G.  referred  to  as  the  cross  Green's  matrix, 
with  matrix  elements  being  the  integrals  (G)xj.  Finally,  one  substitutes  these  matrices 
into  eq .  (38)  and  calculates  the  matrix  5“'.  ideally,  the  result  will  be  a  diagonal  matrix, 

of  which  the  negative  of  the  inverse  of  each  diagonal  element  Is  the  relative  excess 
dielectric  permittivity,  in  the  corresponding  cell  of  the  target  divided  by  the 

background  dielectric  constant,  as  described  by  eq .  (18). 

Cone  I  us  I  on 

We  have  developed  a  three-d imens i ona i  scattering  formulation  from  which  one  may  calculate 
the  scattered  vector  field  anywhere  outside  the  scattering  target;  the  scattering  target  may 
be  of  arbitrary  three-d imens i ona i  shape  with  arbitrary  dielectric  profile.  We  also 
developed  an  inverse  scattering  formulation  by  inverting  the  scattering  equation.  Tne 
inverse  scattering  formulation  provides  a  formula  to  calculate  the  three-dimensional 
dielectric  permittivity  inside  the  scattering  target  from  measured  data  of  the  scattered 
field  in  a  chosen  region  outside  the  target.  There  are  two  crucial  aspects  that  enable  our 
formulation.  One  is  the  d i s s i pa t  i  v  i  t y  of  the  background  medium  that  makes  the  scattering 
problem  invertible  and  provides  the  uniqueness  of  the  inverse  solution.  Another  is  the 
conversion  of  the  integral  equation  into  matrix  equation  by  digitizing  the  space.  Tne 
d  I  ss  I  pat  I V 1 ty  aspect  is  the  theoretical  basis  of  the  formulation.  whereas  the  space 
digitizing  provides  practical  support.  Though  the  formulation  provides  a  neat  equation  for 
inverse  scattering  calculation,  two  practical  probi ems  st'M  require  further  investigation. 
One  is  the  stability  of  the  ma  trix  inversion  of  the  scattered  ma  trix  whici.  is  defined 

in  eq.  (35).  Another  is  the  s « ze  of  the  matrices  that  are  involved  in  the  inverse 
scattering  f  o  rmu la.  viz.,  the  number  of  ceils  into  wh ich  the  target  is  divided.  we  sha  i  i 
d  I  scuss  these  two  problems  in  the  following  paragraphs. 

On  the  stability  of  ma  trix  inversion,  the  question  is:  How  stable  is  the  inverse  of 
if  there  is  an  error  in  This  subject  may  be  investigated  by  analyzing  the  so  cai'ea 

condition  number  of  the  ma  trix  F  < • *  .  wh i ch  i s  a  measure  of  singularity  of  the  ma  t  r i x .  ^ »  The 
Singularity  of  a  matrix  may  be  reduced  by  eliminating  s yriTne try  of  the  matrix,  viz., 
el iminating  the  interdependency  between  the  columns  or  the  rows.  Each  column  of  tno 
scattered  matrix.  F<»>.  corresponds  to  an  incident  wave  while  each  row  corresponds  to  a  cet 
in  the  scattering  region,  v.  Note  that.  for  all  incident  waves,  the  region  v  ana  its 
subdivision  of  ceils  remain  the  same.  One  way  to  increase  asymmetry  of  the  matrix  is  to 
select  the  incident  directions  in  a  random  fashion.  such  as  using  three- d imens i ona i 
directions  with  unequal  angular  spacings.  Similarly.  random  division  of  the  scattered 

region  into  the  N  cells  may  also  increase  asymmetry.  The  latter  means  that  random  placement 
of  the  elements  of  the  receiving  antenna  array  is  better  than  symmetric  placement. 

As  to  the  size  of  the  matrices,  it  is  limited  by  two  conditions.  One  'S  the  desired 
resolution,  and  another  is  the  approximation  that  we  used  to  digitize  the  original  iritegrai 
equation.  The  first  limitation  is  a  subjective  one  which  needs  not  be  remarked.  As  to  the 
second  limitation,  recall  that  in  derivi ng  eq.  (16).  it  was  assumed  that  the  values  of  the 
total  field  and  the  dielectric  permittivity  inside  each  ceil  in  the  target  region  do  not 

vary  much .  For  the  field  to  rema in  constant  inside  each  cell,  the  ceil  must  be  sma i ' er  than 

a  quarter  wavelength  i ns : de  the  target.  For  a  3  GHz  incident  wave,  the  wavelength  inside  a 
typical  biological  target  is  about  12  mm.  so  the  dimension  of  each  ce i i  must  be  sma i 'er  tnar 
3  mm.  This  will  require  an  enormous  computat ion  task  for  a  typical  biological  target  of, 
say.  30  cm  in  size.  One  way  to  reduce  the  size  of  the  matrices  is  to  use  a  lower  frequency 

source.  Using  a  1  GHz  source  will  reduce  the  size  of  the  matrix  by  a  factor  of  27.  and 

thereby  reduce  the  number  of  matrix  elements  by  a  factor  of  729.  Another  approach  15  to 
divide  the  target  into  several  regions  for  image  recons truci*  i on  by  shielding  the  microwave 
radiation  to  all ow  only  limited  exposure.  Theor et ically,  a  division  into  2  regions  reduces 
the  matrix  elements  of  each  region  by  a  factor  of  4,  thereby  a  net  gain  of  2  in  computation 
reduct  ion.  However ,  f  r  om  the  eng i neer 1 ng  point  of  v i ew .  limiting  the  exposure  requires 
additional  technical  consideration  on  i nter -reg i ona I  interference. 
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Any  total  field,  i.e.,  an  incident  and  its  scattered  fields  f‘'*+f‘“'. 

The  incident  field. 

The  scattered  field. 

The  value  of  f  in  the  jth  cell. 

The  value  of  f‘'»  in  the  jtn  ce i  .  . 

The  value  of  f‘«*  in  the  jtn  cell.  in  some  occasions  the  subscript  j  refers  to 
the  jth  incident  direction  (see  foj‘‘’  below),  then  fj‘“’  represents  the  M- 
dimensional  column  vector  derived  from  measurements  of  the  corresponding  scattered 
field  in  any  region  V  outside  of  the  target  which  is  divided  into  N  cel  Is  for  data 
samp  I  I ng . 

The  N-d imens I ona I  vector  representing  the  sampling  of  the  total  field  in  Vo  • 

The  N-d imens i ona I  vector  representing  the  sampling  of  the  incident  field  in  v^ . 

The  N-d imens I ona I  vector  representing  the  sampling  of  the  scattered  field  in  Vo ■ 
The  N-d imens i ona I  column  vector  f„  for  the  jth  incident  direction;  incident  fields 
are  sampled  for  a  total  of  N  different  cases  <j  -  1,  2.  ...  N).  each  with  a 
different  incident  direction  but  otherwise  same  target  setup. 

The  N-O I  mens ional  vector  represent! ng  t  he  samp  ling  of  the  total  field  in  Vo. . 

The  N-d imens I ona I  vector  representing  the  sampling  of  the  incident  field  in  Vo. . 

The  N-O imens I ona I  vector  representing  the  sampling  of  the  scattered  field  in  Vo.. 
The  Green's  matrix  with  matrix  elements  defined  as  ;g,o( x , , x ' ) dx '  with  x,  being  the 
center  of  the  cell  r,  and  dx'  integrated  over  the  ceil  t,. 

The  matrix  element  of  G,  defined  as  XGo,( x, , x ' ) dx '  with  x,  being  the  center  of  the 
cell  T.  and  dx'  integrated  over  the  cell  t,. 

Green's  operator  of  Ht,  (GrrJAr^  *  fAt,Gim  *  1  in  the  space  of  outgoing  scattering 

fields. 

Another  expression  of  C5m  based  on  a  different  factorization  of  the  product  G>t,'S'  » 
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The  Green's  function  of  NU,,  viz.  Go,  in  the  x-representat  i  on . 

)  The  Green's  function  of  lA,, '  ,  viz.  Go,'  in  the  x-representat  lon. 

The  N  X  N  square  submatrix  of  G  with  the  matrix  elements  G, j ,  where  both  \  and  j 
refer  to  ce I  I s  i n  Vo . 

2,  ...)  The  N  X  N  square  submatrix  of  G  with  the  matrix  elements  G, , , 

where  i  refers  to  the  cell  r,  in  v*  and  j  refers  to  cell  r,  m  Vo. 

The  wave  vector  in  the  background  medium. 

The  jth  incident  wave  vector  in  the  background  medium  (ikji  =  iko,i):  if  the  target 
is  divided  into  N  cells,  then  N  scattering  waves  are  used,  each  with  a  different 
incident  direction. 

Maxwell  operator  with  background  medium  and  scattering  target. 

Maxwell  operator  in  a  homogeneous  background  medium. 

Another  expression  of  NV.,  based  on  a  different  factorization  of  the  product  Go,'S'  = 
Go,S.  Where  Gw'  and  Gw  are  inverse  of  NU'  and  NU, .  respectively,  m  the  subspace  of 
all  outgoing  scattered  waves. 

Total  number  of  cells  into  which  the  target,  or  a  region  outside  of  the  target,  is 
divided  for  data  sampling. 

Scattering  part  of  M  (-S  =  M  -  NU,).-  it  depends  on  t  and  u.  For  the  case  where  f 
IS  the  electric  field  E.  then  S  =  - ( < (x) -<wl r cw . 

Another  expression  of  S  based  on  a  different  factorization  of  the  product  Gw'S'  = 
G^S  . 

Negative  of  the  relative  excess  value  of  <  m  the  target  with  respect  to  the 
background  medium.  - t « (x j ) -<„ 1 / <w.  It  is  the  value  of  S  in  the  jth  ceil  for  the 
case  where  f  is  the  electric  field  E. 

A  three-dimensional  vector  representing  a  point  in  space. 

The  vector  representing  the  center  of  the  jth  cell. 

The  region  occupied  by  the  target;  also  represents  its  volume.  V„  is  divided  into 
N  cells  for  data  sampling. 

(a  X  1,  2,  ...)  A  region  outside  of  the  target  in  which  the  scattered  field  is 

sampled.  V„  is  divided  into  N  ceils  for  data  sampling. 

A  function  of  x  representing  the  distribution  of  dielectric  permittivity  in  space. 
Dielectric  permittivity  of  the  homogeneous  background  medium. 

The  electric  scalar  potential. 

A  function  of  x  representing  the  distribution  of  magnetic  permeability  in  space. 
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Magnet  *  c  permeaD i  i 
The  j  t  h  ce i  i  i n  Vo 


tv  of  the  homogeneous  DacKg^ound  medium, 
or  Vo.;  also  represents  its  volume. 
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ABSTRACT 


Theoretical  and  experimental  studies  on  microwave  pulse 
induced  pressure  waves  are  reported.  A  thermodynamical  formula¬ 
tion  of  microwave  interaction  with  dielectrics  is  summarized.  It  is 
shown  that  acoustic  waves  may  be  generated  by  pulsed  microwaves 
even  in  the  absence  of  inhomogeneity  of  microwave  absorption, 
owing  to  discontinuities  of^thermodynamical  variables  and  micro- 
wave  exposure  conditions  across  the  dielectric  interface.  The 
formulation  is  applied  to  a  spherical  system  and  some  numerical 
results  are  presented.  Experimental  results  include  measurement 
of  pressure  waves  in  a  cylinder  filled  with  an  aqueous  solution  of 
electrolytes  exposed  to  pulsed  microwaves  and  estimation  of  the 
coupling  efficiency  between  the  liquid  dielectric  and  the  adjacent 
air.  Pressure  scaling  with  energy  per  pulse  is  also  demonstrated. 

INTRODUCTION 


Previous  derivations  of  the  thermoelastic  mechanism  induced 
by  microwave  pulses  in  dielectrics  were  based  on  the  linear  elastic 
equation  with  an  external  excitation  derived  from  the  thermal 
expansion  that  was  produced  by  an  Inhomogeneous  absorption  of 
microwave  energy  [1].  This  approach  did  not  consider  the  proper 
balance  of  the  distribution  of  the  absorbed  microwave  energy 
among  the  internal  thermal  energies  and  the  bulk  kinetic  energy, 
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and  thereby  failed  to  define  all  thermal  coefficients  adequately. 
Furthermore,  all  analyses  have  been  limited  to  either  isolated 
systems  or  completely  constrained  systems,  both  of  which  exclude 
the  possibility  of  coupling  to  the  external  medium. 

To  correct  these  deficiencies,  we  have  approached  the  sub¬ 
ject  from  thermodynamical  considerations  and  established  a  thor¬ 
ough  formulation  of  the  macroscopic  electromagnetic  interactions 
in  dielectrics  [2].  Couplings  across  dielectric  interfaces  are  in¬ 
cluded  through  discontinuities  of  th^  thermodynamical  variables 
and  conditions  of  microwave  exposure.  Results  of  these  couplings 
may  provide  a  method  for  experimental  evaluation  of  the  micro- 
wave  generated  pressure  waves  in  small  dielectric  objects  for 
which  direct  measurement  is  difficult.  In  this  paper,  the  theoreti¬ 
cal  formulation  is  briefly  summarized.  It  is  then  applied  to  a 
spherical  water-air  system  and  some  dynamical  data  are  thereby 
calculated.  Failure  of  the  linear  approach  under  high  peak  power 
exposures  is  predicted. 

These  results  have  important  implications  for  biosystems 
exposed  to  pulsed  fields.  For  example,  ocular  lenses  exposured  to 
pulse  modulated  microwave  fields  exhibit  histopathological  changes 
which  are  quantitatively  and  qualitatively  different  from  CW 
exposures  of  equal  average  power  [3].  A  leading  mechanism  of 
cellular  injury  which  may  be  invoked  by  pulse  modulation  is  ther- 
moetastic  transduction.  A  microwave  pulse,  upon  absorption  in 
water  dominated  dielectrics,  produces  a  transient  pressure  wave. 
This  transduction  does  not  neccessarily  depend  upon  microwave 
gradients  within  the  dielectric  [2j.  The  transduced  pressure  wave 
then  propagates  withitvthe  dielectric  on  the  basis  of  mechanical 
(elastic)  impedance  subject  to  mechanical  and  thermodynamical 
boundary  conditions.  In  order  to  investigate  the  role  of  such 
pressure  waves,  a  series  of  experiments  have  been  performed  to 
study  the  effect  of  pulse  parameters  as  well  as  temperature.  A 
central  theme  in  these  studies  is  the  role  of  energy  per  pulse  as 
related  to  peak  pressures  induced  by  thermoelastic  transductions. 

Pressure  measurements  were  made  in  saline  filled  vessels 
which  maintain  the  geometric  section  and  material  properties  of 
the  chamber  in  which  earlier  ocular  lens  exposures  took  place  [4]. 
In  addition,  measurements  were  also  made  in  vessels  of  similar 
geometric  section  but  different  material  properties.  Pressure 
measurements  were  made  outside  of  the  waveguide  where  expo¬ 
sures  took  place.  Pressure  measurements  were  made  both  in  the 
saline  medium  and  in  the  air  at  the  air/liquid  interface,  in  order  to 
investigate  the  coupling  efficiency  between  the  medium  in  which 
the  transduction  takes  place  (saline)  and  the  adjacent  non- 
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absorbing  medium  (air).  Theoretical  calculations  for  a  one¬ 
dimensional  case  comparable  to  the  vessel  geometry  have  been 
published  previously  [2], 

THE  NON-LINEAR  WAVE  EQUATIONS 


The  objective  is  to  derive  the  governing  equations  for  the 
thermodynamical  and  mechanical  variables,  which  include  mass 
density  p,  stress  tensor  p,  thermodynamic  internal  energy  U|,  bulk 
kinetic  energy  U|^,  bulk  velocity  v,  temperature  T,  and  microwave 
specific  absorption  rate  (SAR)  P.  Since  the  details  of  the 
formulation  have  been  published  earlier  [2],  here  we  just  give  a 
brief  summary.  The  time  domain  of  our  concern  is  the  acoustic 
relaxation  time,  which,  as  we  shall  show,  is  of  the  order  of 

milliseconds.  We  assume  that  this  is  much  greater  than  the 

microwave  pulse  width,  which  in  turn  is  much  greater  than  the 
period  of  the  microwave  carrier  and  the  time  to  reach  a  local 
equilibrium.  We  also  assume  that,  as  it  will  be  verified  later,  the 
acoustic  relaxation  time  is  in  general  much  smaller  than  the 

relaxation  time  to  reach  a  thermal  steady-state  by  heat 

conduction.  Based  on  these  assumptions,  we  shall  ignore  heat 
conduction,  heat  transfer  by  convective  flow,  and  any  dependency 
of  U|  on  electric  polarization.  To  reduce  complexity,  we  shall  also 
ignore  the  internal  energy  associated  with  the  anisotropicity  of  the 
strain  tensor.  It  can  then  be  shown  that  the  thermoacoustic  effect 
may  essentially  be  derived  from  the  laws  of  the  conservation  of 
mass,  the  conservation  of  momentum,  the  conservation  of  energy. 


and  the  thermodynamic  equation  of  state: 

dtP^+  pV*v  =  0  (1) 

pdtVj  +  3jPjj  =  0  (2) 

Pi  j  3jVj  +  pdjU,  =  pP  (3) 

^t*^!  =  (  Vp^JP;  j  3j  v, /V*v]djp 


where  d^  =  3j  +  y*V  is  the  convective  derivative,  and  Cp  are 
the  specific  heat  (per  mass)  at  constant  volume  and  constant 
pressure,  respectively,  8p  is  the  isobaric  thermal  expansion  coeffi¬ 
cient,  and  Kj  the  isothermal  compressibility.  The  last  equation, 
which  may  be  substituted  by  another  thermodynamic  equation  of 
state,  is  derived  from  the  general  expression  Uj  =  U[p,pJ  by  using 
the  Jacobian  method  of  relating  one  partial  derivative  to  another 
and  applying  the  first  law  of  thermodynamics. 
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By  mutual  substitution  among  equations  (1)  through  (4),  one 
may  obtain  the  wave  equations  for  any  of  the  thermodynamic 
quantities.  In  particular,  if  the  strain  tensor  3{Uj  is  isotropic,  one 
obtains  the  following  wave  equation  for  the  stress  tensor: 

=  3,llP<,/p)''^tBp/SjPl  (5) 

Similar  wave  equations  for  the  njass  density  p  or  the  velocity 
V  may  also  be  obtained.  These  equations  must  be  supported  by  the 
respective  boundary  conditions,  which  are  implicated  in  eqs.  (1  }- 
(3).  All  boundary  conditions  may  be  obtained  by  integrating  each 
term  of  these  equations  across  a  thin  layer  of  the  boundary 
interface  and  applying  the  generalized  Stoke's  theorem.  However, 
since  the  boundary  surface  is  not  stationary,  one  must  first 
transform  these  equations  to  the  Lagrangian  specification  to  obtain 
the  boundary  conditions  [2],  The  wave  equations  and  boundary 
conditions  are  greatly  simplied  in  the  one  dimensional  case.  Note 
that  in  the  one-dimensional  case,  the  factor  in  ep*(5)  is 

replaced  by  (.Pg/pJ“^  •  i^or  our  calculations  later,  we  list  the 
result  for  the  pressure  and  the  bulk  velocity: 

with  the  boundary  conditions  that 

p(x,0  and  Pq’^  3jjP(x,t) 
are  continuous  across  the  dielectric  interface,  and, 

with  the  boundary  conditions  that 

v(x,t)  and  (pCp3^v  -  Pq  SpP )  /  (  P^C^  icj ) 
are  continuous  across  the  interface. 

These  wave  equations  show  that  acoustic  waves  are 
generated  either  by  an  inhom ogeneity  or  temporal  variation  of  the 
product  of  the  microwave  SAR  and  some  thermodynamical 
variables. 
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LINEAR  APPROXIMATION  IN  A  SPHERICAL  MODEL 


To  illustrate  the  implications  of  the  wave  equations  and  their 
respective  boundary  conditions,  we  consider  a  dielectric  sphere  of 
radius  a  surrounded  by  air  at  1  atmospheric  pressure.  A  square 
pulse  of  microwave  of  duration  t  and  amplitude  P^  is  incident  upon 
the  sphere.  The  sphere  is  assumed  to  be  small  so  that  the 
microwave  absorption  exhibits  no  spatial  variation.  In  the  linear 
approximation  eq.  (5)  gives,  for  the  pressure  waves  in  the 
frequency  domain;  * 


p^(r,<u)  =  ^ 

P 


{- 


ipC^i  O _ 

“lyllc^a) 


-1 


(8) 


and 


P2(r,u))  =  -•(— ^J-1  TTjirpry 

io(Jliij _ 

T^TlTTiy-— 


{- 


jo(h»)  -  h^'lkjay  io'(h») 


(9) 


where  jf,(x)  and  h^(x)  are  respectively  the  spherical  Bessel  and 
Hankel  functions  of  the  flTst  kind  of  order  n,  the  subscripts  1  and  2 
label  the  dielectric  and  the  air  media,  respectively,  c-)  ^  denotes 
the  corresponding  value  of  /Cp/  (  pC^  Xy)  in  either  medium,  which 
is  the  velocity  of  the  acoustic  wave  in  the  medium,  tan(J)  = 
(pc)2/(pc)i  is  the  ratio  of  the  acoustic  impedances  of  the  air  to  the 
liquid,  ki  2  ”  acoustic  wavenumber  in  either 

medium,  and  P(a))  is  the  Fourier  transform  of  the  square  pulse  SAR, 
which  is  given  by: 

P(«1  >  -  p„  ID)  no) 


The  factor  (u)  -  iO)  in  the  denominator  represents  a  singularity  at 
the  origin  immediately  below  the  real  axis  of  the  complex  (j*-plane. 
Multiplying  P'1  ^  2(  p“)  e"'*^^  and  integrating  over  (uusing  the 

contour  integration  technique,  eqs.  (8)  and.  (9)  then  give  the 
pressure  waves  in  the  time  domain: 
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«  _  0$  P  [  <o_  I 

^-rVh  exp(-i'^n^) 

n-0  pin 

(C2/C1  -  tan4>)  /  (iOj^a/ci  )^  +  tanij)  c s c^ ( w^a/c ^  ) 


and 


0?  _  p8  P ( (i)  It 

P2('''t)  =  L"^  t-rVJ,  ir^  exp(-ia)^t) 
n-o  pin 

^  [i  -  C2/(<o„a)J  h^(ai„r/c2)/ho(ui„a/c2) 
(c^/c^  -  tanij)) /((Uj^a/c-i  +  tarnj)  c  s  ( lo^^a/c  ^  ) 


where  oj^  are  the  frequency  eigenvalues  which  are  given  by  the 
solutions  of  the  equation 

[1  -  ((oa/cj)  cot(iJJa/c  1)  ]  tanifi  =  ~  •(‘“•‘/ci)  (13) 

For  a  water-air  system,  tanij)  is  about  3  x  10~^,  the  above 
equation  may  then  be  solved  by  making  the  Taylor  expansion  for  oj 
as  a  funcion  of  tan<j).  One  gets,  to  the  first  order  in  tan<|). 


(n/O) 


(C2/C,  )  /(mr)  +  i 

nTT  -  - = - - - - - -  taniji 

1  +  {c2/c^  )^/(mr)2 


(n=0) 


I2I 


to  to  t  9 

-H  -sf  -l)tan^J 


(14) 


which  gives  an  acoustic  relaxation  time  of  t  =  1  /lm(t^  )  =  11.2 
milliseconds  for  the  fundamental  mode  in  a  water  sphere  with  2a  = 
1  cm.  This  is  much  smaller  than  the  time  for  the  sphere  to  reach  a 
thermal  steady  state  by  heat  conduction,  which  is  about  19.4 
seconds  [5].  Since  the  acoustic  relaxation  time  is  proportional  to 
the  radius  a  while  the  conductive  relaxation  time  is  proportional  to 
a^,  the  two  relaxation  times  become  closer  as  the  size  of  the 
dielectric  decreases,  and  they  become  equal  when  the  radius  is 
about  2.9  microns,  beyond  which  thermal  conduction  becomes 
Im  portant. 

To  obtain  the  acoustic  energy  coupled  to  the  air,  one  first 
calculates  the  solution  for  the  velocity  in  the  frequency  domain. 
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which  is  givenby  [-i/(<*ip)J^P^ 


The  total  acoustic  energy  coupled  into  the  air  may  be  obtained  by 
calculating  the  work  done  by  the  dielectric  on  the  air  at  the 
interface.  To  the  first  order,  the  result  is: 

a> 

^air  “  v^(a,u)=0) 

=  /IIT  Pot^il 

Comparing  this  to  the  total  microwave  absorption  by  the  liquid, 
Ejbs  "  /Tif (4Tra^/3)PiP(a)=0)  ,  one  obtains  the  coupling  efficiency 
of  the  absorbed  microwave  energy  to  the  acoustic  energy  in  the 
air: 

"lair  =  ^air^^abs  =  Po^-'^Jr 

Take  the  dielectric  to  be  water,  so  8  =  2.8  x  1 0““*  ®C~^  at  30 

®C  [6],  Cp  =  4.186x10^  erg/gm-°C,  and  p  =  1  gm/ cm^ .  With  the 
equilibrium  pressure  of  the  air,  being  1  atm  =  1.01  x  10^ 
dyne/cm^,  the  coupling  efficiency  is  then  rijjj.  =  6.7  x  10~^.  This 
value,  as  we  shall  show  in  the  next  section,  appears  to  agree  with 
the  results  from  pressure  measurements. 

PRESSURE  MEASUREMENTS 
Material  and  Methods 

Pressure  measurements  were  made  in  both  air  and  saline 
media  for  two  types  of  vessels  placed  in  the  WR  975  exposure 
system  (cf.  Fig.  1).  The  first  vessel  duplicated  the  cross  section 
and  material  properties  of  the  lens  exposure  chamber.  That  is,  it 
was  made  of  square  quartz  tubing  ca.  one  inch  on  each  side  with  a 
1/8  inch  wall  thickness.  The  second  vessel  was  circular  in  cross 
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Fig.  1.  Functional  block  diagram  of  the  WR  975 
exposure  system. 

section,  one  inch  in  outer  diameter,  1/8  inch  wall  thickness,  and 
made  of  acrylic  plastic.  Both  were  filled  with  phosphate  b offered 
saline  (PBS).  These  vessels  were  placed  into  the  WR  975  via  the 
same  waveguide  below  cut-off  window  as  the  actual  lenticular 
exposure  chamber.  The  two  vessels  extended  beyond  the  wave¬ 
guide  to  provide  a  PBS  column  for  coupling  of  the  acoustic  wave 
into  regions  where  the  hydrophone  would  not  be  subject  to  direct 
influence  of  the  electromagnetic  field.  Each  vessel  was  separately 
matched  to  the  microwave  source  with  the  triple  stub  tuner. 
Typical  return  losses  were  ca.  -15  to  -20  dB  for  either  vessel. 

Pressure  measurements  made  in  the  liquid  were  performed 
with  the  equipment  illustrated  in  Fig.  2.  The  transducer  was  of 
Bruel  and  Kjaer  manufacture,  the  type  8103  miniature  hydrophone. 
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Fig.  2.  Functional  block  diagram  of  the  acoustic 
measurement  system. 
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This  transducer  is  of  the 
piezoelectric  type  as 
shown  in  Fig.  3a.  It  is  9 
mm  in  outer  diameter 
with  an  integral  water¬ 
proof  connection  to  6 
meter  of  double  shielded 
cable.  The  type  8103 
has  a  200  KHz  band¬ 
width  with  a  directivity 
response  as  shown  in 
Fig.  3b.  Since  the  hy¬ 
drophone  has  serious 
lobing  when  used  with 
axial  incidence,  an  open 
mesh  (ca  2  by  2  mm)  fi¬ 
berglass  screen  was 
loosely  wrapped  around 
the  transducer  to  center 
it  in  the  vessel's  cross 
section.  The  hydrophone 
was  calibrated  with  a  B 
&  K  type  4223  hydro¬ 
phone  calibrator.  Ail 
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Fig.  3a.  The  hydrophone  construc¬ 
tion  and  dim  ensions. 
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Fig.  3b.  The  axial  and  transverse  patterns  of  the  hydrophone. 
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water  based  measurements  are  expressed  in  dB  relative  to  0  dB  =  1 
micro  Pascal.  The  hydrophone  was  interfaced  to  a  charge 
amplifier  also  of  B  &  K  manufacture,  the  type  2650.  This 
instrument  allows  insertion  of  the  hydrophone  calibration  factor 
into  the  signal  conditioning  network  just  prior  to  the  output 
amplifier.  When  used  in  the  charge  mode,  the  2650  input  amplifier 
is  bypassed  and  the  first  stage  becomes  a  band-pass  filter.  Prior  to 
the  2650,  the  hydrophone  cable  was  terminated  by  an  Erie  low  pass 
filter  with  30  dB  rejection  at  1  GHz.  Hydrophone  wave-forms  were 
recorded  by  means  of  a  Hewlett  ^Packard  1740  osciloscope 
connected  to  the  output  of  the  B  &  K  2650.  Photographs  were 
made  with  a  Hewlett  Packard  1976  osciloscope  camera  and 
Polaroid  type  107  film.  The  photographic  exposure  duration  was 
one  second,  thereby  each  photograph  contains  a  total  of  5  pulses 
superim  posed. 


Pressure  measurements  in  air  were  made  with  a  Bruel  and 
Kjaer  sound  pressure  level  (SPL)  meter,  the  type  2230.  The 
microphone  was  a  0.5  inch  condenser  type  of  B  &  K  manufacture, 
the  type  4155.  The  SPL  meter  and  microphone  were  shielded  from 
direct  influence  of  the  electromagnetic  field  and  ambient  noise  by 
a  1  m^  enclosure  of  microwave  absorber  Eccosorb^  EC-77.  This 
combination  was  calibrated  with  a  microphone  calibrator  the  B  &  K 
type  4230  for  an  average  SPL  of  1  Pascal.  All  air  based  pressure 
measurements  are  made  in  dB  relative  to  0  dB  =  20  micro  Pascal. 
The  2230  must  be  used  with  the  detector  in  the  peak  responding 
mode  with  display  of  the  maximum  value  using  broadband  (lin) 
frequency  response  and  fast  time  weighting.  Electrical  analogs  of 
the  pressuie  wave  forms  are  available  from  the  2.5  mm  AC  out 
minijack  on  the  side  of  2230.  Waveform  display  and  recording  was 
performed  using  the  same  system  described  for  the  hydrophone 
measurements.  Under  these  conditions,  the  rise  time  of  the  system 
is  ca  30  microseconds. 


The  microwave  pulse  waveform  (see  Fig.  4)  was  measured 
from  the  forward  power  arm  of  the  Dielectric  Communications  WR 
975  directional  coupler  by  means  of  a  Hewlett-Packard  directional 
detector  HP  type  423A.  All  RF  pulses  were  monitored  for  rise 
time  (ca  0.5  microseconds)  and  stated  pulse  widths  are  RF  pulse 
widths  rather  than  video  pulse  widths,  i.e.,  the  video  rise  time  is 
excluded.  The  standard  RF  pulse  was  10  microseconds  in  duration 
with  a  peak  transmitted  power  of  10  KW  at  918  MHz.  This 
produced  an  energy  density  of  ca  326  micro  Joules/cm^. 
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Fig.  4.  A  detected  microwave  pulse  is  shown  with  5 
psec  per  horizontal  division.  The  negative 
detector  (HP  type  423A)  produces  a  down¬ 
ward  pulse. 

Both  air  and  PBS  pressure  measurements  were  tested  for 
direct  electromagnetic  field  effects.  In  the  case  of  air 
measurement  this  was  accomplished  by  recording  the  SPL  with 
blockage  of  the  microphone  air  path  by  ca  20  dB  with  foam 
rubber.  Also  the  SPL  meter  body  was  covered  with  EC-77,  in  the 
case  of  the  PBS  based  measurements,  the  hydrophone  transducer 
was  removed  from  the  PBS  coiumn  with  all  other  factors  left 
constant,  thus  blocking  the  water  path  by  ca  60  dB. 

Results 

The  pressure  measurements  recorded  with  the  hydrophone  are 
shown  in  Fig.  S  as  a  functren  of  position  in  the  vessel.  The  peak 


Fig.  5.  Peak  pressures  as  a  function  of 
distance  outside  of  the  waveguide. 


I  Voll  •  lOOO 


416 


Guo  e(  al. 


pressures  recorded  were  independent  of  the  vessel  material 
composition  for  equal  net  energy  absorbed  per  unit  volume  of 
liquid.  Note  that  no  pressure  measurements  were  made  inside  the 
wave  guide.  The  pressure  maximum  (exclusive  of  the  wave  guide 
interior)  was  at  a  distance  of  about  1  cm  up  from  the  flange  of  the 
wave  guide  below  cut-off  window.  At  this  location  the  pressure 
had  a  peak  value  of  ca  10  kilopascals.  The  peak  pressure  declined 
nonmonotonically  with  distance  to  about  1  kilopascal  at  the 
air/water  interface. 

t 

The  pressure  waveform  recorded  with  the  acoustical 
measurement  system  is  shown  in  Fig.  6.  The  pressure  wave  rise 
time  in  the  PBS  is  estimated  in  Fig.  7  to  be  ca  30  microseconds. 
This  particular  figure  shows  the  pressure  rise  time  for  the  20 
microsecond  RF  pulse,  but  that  for  the  10  microsecond  RF  pulse 
was  substantially  identical  (probably  due  to  the  rise  time 
limitations  of  the  modulator  and  hydrophone  systems).  The 
pressure  fall  time  to  the  abcissa  is  ca  30  milliseconds.  This  is 
followed  by  a  slow  (ca  50  milliseconds)  pressure  wave  undershoot  of 
modest  amplitude  before  returning  to  the  baseline. 


Fig.  6.  The  pressure  waveform  is  shown  above  with 
100  pascals  per  vertical  division  and  20  msec 
per  horizontal  division. 


Pressure  scaling  with  energy  per  pulse  is  demonstrated  for 
the  standard  puise  and  doubling  the  energy  per  pulse  by  doubling 
the  pulse  width  from  10  to  20  microseconds  (see  Fig.  8).  Similarly, 
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linear  scaling  of  the  peak  pressure  is  demonstrated  in  Fig.  9  for 
doubling  the  energy  per  pulse  by  doubling  the  peak  power  leaving 
the  pulse  width  constant  at  10  microseconds.  Such  pressure  scaling 
is  consistent  with  the  observation  of  increased  histopathological 
damage  with  increased  energy  per  pulse  for  fixed  average  power 
and  total  absorbed  energy  [7]. 


Fig.  9.  Pressure  scaling  with  energy  per  RF-pulse  is 
shown  for  the  case  of  doubling  the  peak  power 
transmision  from  lOKW  to  20KW.  The  axes  and 
energy  densities  are  those  shown  for  Fig.  8. 


The  air  pulse  is  shown  in  Fig.  10.  The  displayed  maximum 
SPL  in  air  was  90.4  dB,  with  standard  deviation  2.26  dB.  This 
suggests  a  ca  64  dB  coupling  factor  from  liquid  pressure  to  air 
pressure  across  the  interface.  This  confirms  the  predicted  coupling 
efficiency  presented  earlier  in  this  report.  The  rise  time  appears 
to  be  ca  5  microseconds  (also  probably  limited  by  the  SPL  meter's 
band  width)  and  a  fall  time  of  .5  milliseconds. 


Fig.  10.  The  pressure  pulse  in  the  air,  with  500  psec 
per  horizontal  division  and  0.2  volt  per 
vertical  division. 
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REMARKS  AND  DISCUSSIONS 


We  have  shown  that  acoustic  waves  arc  generated  if  there  is 
any  spatial  discontinuity  in  the  dielectric  medium,  or  temporal 
discontinuity  in  the  microwave  absorption.  Coupling  of  the 
acoustic  energy  from  an  absorbing  medium  to  a  nonabsorbing 
medium  has  been  formulated  through  the  thermodynamic  equations 
and  the  boundary  conditions.  The  result  may  be  useful  for  the 
evaluation  of  the  pressure  waves  inside  small  biological  objects  by 
measuring  the  pressure  waves  or  other  thermodynamical  quantities 
in  the  surrounding  medium. 

The  result  of  the  linear  approximation  indicates  that  the 
percentage  of  the  absorbed  microwave  energy  being  coupled  to  a 
non-absorbing  medium  is  proportional  to  the  pressure  in  the 
nonabsorbing  medium  and,  apart  from  that,  is  independent  of  its 
other  thermodynamic  properties.  It  is  also  independent  of  the 
microwave  pulse  width,  so  that  total  acoustic  energy  coupled  to 
the  air  is  proportional  to  the  total  microwave  absorption  cross 
section.  The  predicted  value  of  6.7  x  10~^  for  the  coupling 
efficiency  (cf.  eq.  (18))  is  exactly  the  same  as  the  value  derived 
previously  for  the  one-dimensional  case  (2),  which  appears  to  agree 
with  the  experimental  result  of  the  pressure  measurements 
described  in  the  last  section.  Indeed,  comparing  eqs.  (17)  and  (18) 
to  the  corresponding  equations  for  the  one-dimensional  case,  we 
arrive  at  the  conclusion  that  the  coupling  efficiency  is  independent 
of  the  geometry  of  the  system. 

The  formulation  indicates  that  the  microwave  induced  ther¬ 
moacoustic  effect  is  in  gei^eral  nonlinear.  At  low  radiation  den¬ 
sity,  one  may  make  the  linear  approxim  ation  by  ignoring  all  second 
and  higher  order  terms.  The  validity  of  the  linear  approximation 
can  be  evaluated  by  calculating  the  generated  pressure  wave  and 
comparing  it  to  the  initial  equilibrium  value.  Taking  an  optimal 
pulse  width  of  t  =  it/wq,  the  amplitude  for  Pi(r,t)  in  eq.(11) 
corresponding  to  the  fundamental  frequency  component  is  equal  to 
(2a/Tr)(  c8p/Cp)  iPg.  Using  the  values  of  Bp  and  Cp  for  30  ®C  as 
cited  below  eq.  (18)  and  assuming  1  cm  diameter  of  a  water  sphere 
surrounded  by  air  of  1  atm  pressure,  a  peak  SAR  of  =  30  kW/gm 
results  in  a  pressure  wave  of  amplitude  equal  to  0.1  bar  in  the 
water,  which  is  10%  of  the  initial  equilibrium  pressure.  This 
fraction  of  variation  in  pressure  will  also  result  in  at  least  equal 
fraction  of  error  in  the  linear  approximation.  For  the  one 
dimensional  case  [2]  with  4  cm  length  of  water,  which  maybe  used 
to  approximate  the  dielectric  vessel  described  in  the  previous 
section,  the  same  pressure  amplitude  is  generated  with  a  peak  SAR 
of  3  kW/gm.  The  above  estimate  is  based  only  on  the  fundamental 
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harmonic.  Eq.  (11)  (as  well  as  eq.  52  of  ref.  2)  shows  that  the 
pressure  amplitudes  of  the  higher  harmonics  are  approximately 
inversely  proportional  to  the  frequency,  so  the  sum  does  not 
converge  rapidly.  Therefore  the  total  pressure  amplitude  could  be 
much  greater  than  that  estimated  with  the  fundamental  harmonic 
alone.  This  is  indeed  evidenced  by  the  pressure  measurements 
described  in  the  preceding  section. 

The  experimental  portion  of  this  work  is  subject  to  some 
limitations  which  merit  discussion,  fhe  first  of  these  is  the  fact 
that  no  measurements  were  made  inside  of  the  waveguide.  The 
reason  for  this  exclusion  was  to  reduce  artifact  due  lo  direct 
electromagnetic  coupling.  The  consequence,  however,  is  that  we 
have  no  measurement  of  pressure  at  the  location  where  the  ocular 
lens  specimen  was  placed.  We  presume  that  pressures  at  that 
location  may  be  higher  than  those  recorded  outside  of  the 
waveguide.  Also  the  acoustic  mode  distribution  in  the  vessel  is 
complex  as  judged  by  the  nonmonotonic  pressure  attneuation  with 
increasing  axial  distance,  and  as  evidenced  by  the  marked  pressure 
gradients  in  the  transverse  plane  of  the  vessel.  Work  is  presently 
in  progress  to  develop  a  subminiature  hydrophone  based  upon 
optical  transduction  mechanisms  which  will  be  free  of  direct 
electromagnetic  influence.  This  probe  will  be  used  not  only  inside 
the  waveguide,  but  also  in  vivo  for  direct  pressure  measurements. 

Another  instrumental  limitation  of  importance  is  the 
bandwidth  of  the  transducers.  This  imposes  a  lower  limit  in  the 
pressure  wave  rise  time  estimates.  Nevertheless,  the  recorded 
pressure  rise  time  in  PBS  of  30  microseconds  is  remarkably  fast 
when  compared  to  pressure  rise  times  in  physiologic  situations.  It 
is  also  very  fast  compared  to  more  common  overpressure  situations 
due  to  air  conducted  impulsive  noise  or  blast  effects.  Of  course, 
air  conducted  impulse  noise  is  poorly  coupled  to  internal  organs 
such  as  the  ocular  lens  where  no  air  path  exists  (this  stands  in 
contrast  to  the  situation  of  the  so-called  hollow  viscus  or  air  filled 
space  such  as  the  lungs.)  We  speculate  that  the  pressure  rise  time 
may  play  a  role  in  cellular  damage  that  is  separable  from  the  peak 
pressure.  This  factor  may  not  be  examined  in  the  present  series  of 
experiments  since  the  exposure  system  is  a  power  oscillator  with 
fixed  rise  time  of  ca  0.5  microseconds. 

The  relationship  of  this  work  to  that  of  Olsen  and  Lin  [8] 
must  also  be  discussed.  The  pressure  measurements  reported  in 
that  paper  were  directed  primarily  at  a  determination  of  the 
acoustic  modes  of  the  head  of  various  experimental  animals,  rather 
than  absolute  pressure  and  air  coupling  measurements  as  was  the 
case  in  our  work.  This  is  evidenced  by  the  fact  that  Olsen  and  Lin 
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made  no  attempt  to  locate  a  pressure  maximum.  In  fact,  coupling 
was  intentionally  minimized  to  reduce  artifact.  Furthermore,  no 
cognizance  was  taken  of  the  spatial  response  pattern  of  the 
piezoelectric  disk  transducer  or  its  directivity  with  respect  to 
presumed  acoustic  generator  centers.  Nevertheless,  a  calibration 
curve  was  provided  along  with  the  transducer  output  signals.  The 
calibration  curve  is  somewhat  difficult  to  interpret  jlue  to  axis 
mislabeling,  but  if  we  assume  that  sensitivities  (decibels 
referenced  to  one  volt  per  micropascal)  to  be  preceded  by  a 
negative  sign,  and  that  opposite  label  is  Pascals  per  millivolt,  then 
the  pressure  traces  may  be  quantified.  For  example,  the  trace  for 
rat  2  in  figure  five  (ref.  8)  demonstrates  a  peak  of  about  200 
Pascals.  Since  Olsen  and  Lin  used  the  same  pulse  width  but  ca  3.3 
times  lower  peak  power,  their  result  would  scale  to  ca  650  Pascals 
(neglecting  the  probable  underestimation  of  a  pressure  peak).  Our 
results  show  a  peak  of  ca  1  kilopascal  towards  the  distal  end  of  the 
tube  and  much  higher  values  as  the  flange  of  the  waveguide  below 
cut-off  window  is  approached  (cf.  Fig.  5  of  this  paper).  Based  on 
the  differences  in  geometry  and  the  artifact  problem,  the 
measurements  of  Olsen  and  Lin  are  broadly  consistent  with  ours, 
save  those  near  the  window  flange.  The  measurements  reported 
here,  although  confirming  the  coupling  between  the  absorbing 
(saline)  and  nonabsorbing  (air)  media,  can  not  be  compared  to  our 
theoretical  predictions  in  the  case  of  the  very  high  pressure  values 
recorded  near  the  window  flange.  On  one  hand  this  may  be  due  to 
the  fact  that  our  numerical  estimate  is  derived  from  only  the 
contribution  of  the  fundamental  harmonic.  On  the  other  hand, 
until  independent  confirmation  of  these  high  values  can  be  made 
with  another  transducer  of  better  intrinsic  decoupling  from  the 
incident  field,  the  highest  values  must  remain  suspect. 

The  theory  and  experimental  results  presented  here  may 
represent  the  microwave  analog  of  photoacoustic  spectroscopy 
[9].  Some  earlier  works  with  microwave  (X  band)  sources  for 
electron  paramagnetic  resonance  were  reported  by  Diebold  and 
McFadden  [10].  They  made  neither  measurements  of  coupling 
efficiency  nor  pressure  measurements  in  the  liquid  phase.  Also, 
the  Rosenc waig-Gersho  theory  [9]  appears  to  be  derived  from 
White  [11]  wherein  the  thermoacoustic  elastic  transduction  process 
is  assumed  to  derive  from  spatial  inhomogeneity  of  microwave 
specific  ab sorptio n  rate. 
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Recent  Developments  in  Microwave 
Medical  Imagery — Phase  and  Amplitude 
Conjugations  and  the  Inverse  Scattering 

Theorem 

Theodore  C.  Guo,*  Wendy  W.  Guo,*  and  Lawrence  E.  Larsen^ 


A  theoretical  analysis  of  the  local  field  of  microwave  lattice  radiation  source  and  an  in¬ 
verse  scattering  theorem  are  presented.  The  results  are  applied  to  a  water-immersed 
microwave  array  system  for  medical  imaging.  It  is  shown  that,  using  a  technique  of 
phase  and  amplitude  conjugations,  a  satisfactory  three-dimensional  focusing  for  a  tar¬ 
get  located  in  the  neighborhood  of  the  array  may  be  achieved.  The  focusing  resolutions 
for  transverse  and  longitudinal  directions  are  approximately  X/2  and  X,  respectively, 
where  X  is  the  wavelength  in  the  dielectric.  By  increasing  the  element  spacing  of  the 
array,  the  resolutions  can  be  as  good  as  5.3  mm  and  11. 7  mm,  respectively,  at  the  operat¬ 
ing  frequency  of  3  GHz.  Combining  the  technique  of  phase  and  amplitude  conjugations 
with  the  inverse  scattering  theorem,  the  array  will  be  able  to  provide  a  three-dimen¬ 
sional  imaging  system  with  satisfactory  resolution. 


1.  INTRODUCTION 

The  use  of  imagery  methods  for  microwave  dosimetric 
analysis  has  many  advantages.  Chief  among  these  are  the 
non  invasive  nature  of  the  data  collection  and  the  ease  with 
which  imagery  as  a  form  of  data  display  can  be  related  to  the 
medical  traditions  of  anatomy  and  pathology.  In  fact,  a 
considerable  literature  is  being  developed  or.  the  subject  of 
medical  applications  of  microwave  technology  in,  e.g.,  the 
cardiovascular  system  as  described  elsewhere  in  this  volume 
by  Lin  [I)-(13|.  Nevertheless,  medical  imagery  with  mi¬ 
crowave  radiation  faces  a  number  of  technical  and  theoret¬ 
ical  problems.  Historically,  the  most  prominent  technical 
problem  has  been  the  dilemma  cast  by  the  contradictory 
requirements  of  resolution  and  propagation  loss  with  respect 
to  the  choice  of  the  frequency  of  operation.  When  imagery 
is  formed  by  sampling  the  scattered  fields  in  air,  system 
operation  at  frequencies  below  ca.  10  GHz  produces  unac¬ 
ceptably  low  spatial  resolution  whereas  higher  frequencies 
of  operation  impose  excessive  propagation  loss  [14].  Simul¬ 
taneously,  the  use  of  air  coupling  in  system  design  introduces 
two  additional  undesirable  features:  multipath  contamina¬ 
tion  from  propagation  paths  exterior  to  the  target,  including 
lateral  or  Beverage  wave  propagation  at  the  air-dielectric 
interface  (15|;  and  poor  power  transfer  between  the  target 
and  the  antennas  used  for  data  collection  [16].  These  prob- 
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lems  are  solved  by  operation  of  the  antennas  in  a  medium 
of  high  permittivity  and  loss  tangent  comparable  to  the  bi 
ologic  target  under  study.  The  medium  of  choice  is  water 
since  it  provides  a  resolution  enhancement  by  wave  length 
contraction  of  nearly  a  factor  of  9  at  3  GHz  [3]  while  simul¬ 
taneously  the  impedance  match  is  improved  and  reflection 
is  reduced  in  comparison  to  air  coupling.  Water  coupled 
arrays  offer  all  the  advantages  of  water  r<  upled  elements 
plus  the  enormous  increase  in  data  collection  speed  made 
possible  by  electronic  scanning  and  focusing.  Additional 
discussion  is  presented  in  a  companion  paper  by  Foti  et  al 
elsewhere  in  this  volume. 

In  spite  of  these  operational  advantages,  water  coupling 
introduces  its  own  complications.  The  worst  of  these  is  the 
fact  that  the  anten:'as  must  be  placed  in  close  proximity  to 
the  subject  to  reduce  coupling  losses.  Thus,  the  antennas 
must  operate  within  their  reactive  or  local  zone.  This  is 
especially  troublesome  in  the  case  of  array  antennas  since 
Fresnel  or  Fraunhofer  diffraction  theory  cannot  be  used  to 
effectively  focus  or  stear  the  array.  We  address  this  problem 
in  terms  of  a  new  method  for  array  focusing  which  accom¬ 
modates  not  oniy  local-zone  operation  with  3-dimensional 
focus  control,  but  also  element-to-element  variation  in 
network  parameters.  This  method  is  based  on  a  phase- 
amplitude  conjugation  of  fields  sampled  by  the  array  when 
illuminated  by  a  half-space  omnidirectional  radiator. 

The  major  theoretical  problem  in  medical  application  of 
microwave  imagery  is  recovery  of  resolution  in  the  direction 
if  propagation  for  forward  scattering  (bistatic)  based  sys¬ 
tems.  In  medical  imagery  from  radiologic  disciplines,  this 
is  known  as  the  tomography  problem;  in  microwave  and 
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Fig.  la.  I  .attice  structure  of  the  planar  array  of  12T  elements. 


electromagnetic  propagation  disciplines,  this  is  the  inverse 
scattering  problem.  .4  brief  review  of  the  inverse  scattering 
problem  and  its  historical  solutions  is  presented  elsewhere 
in  this  volume  by  Boerner  and  Chan.  We  express  the  prob¬ 
lem  chiefly  in  terms  of  axial  resolution  since  usable  resolu¬ 
tion  is  obtainable  in  the  transverse  plane  for  isolated  organs. 
Obviously,  operation  in  uivo  would  require  solution  of  the 
inverse  scattering  problem  in  any  case.  We  present  a  new 
inverse  scattering  theorem  which  may  be  considered  as  a 
generalization  of  the  Lorentz  reciprocity  theorem  to  the  case 
of  lossy  media  ll9l-(211.  It  is  applied  to  forward  scattered 
fields  in  combination  with  the  phase-amplitude  conjugation 
in  a  model  data  collection  system  which  simulates  the  DART 
(Dosimetric  Analysis  by  Radiofrequency  Tomography) 
Mark  4  system  under  development  at  the  Walter  Reed  Army 
Institute  of  Research  [2),  13),  '6l-(8l,  Il2l-[13l. 

The  next  section  provides  a  brief  description  of  the  DART 
system.  It  is  followed  by  a  theoretical  analysis  of  the  local 
field  of  the  receiving  array  in  Section  3.  Section  4  introduces 
the  topic  of  phase-amplitude  conjugation  and  the  application 
to  three-dimensional  focusing.  Section  5  presents  the  new 
inverse  scattering  theorem  and  its  use  in  combination  with 
the  phase-amplitude  conjugation  method  for  image  recon¬ 
struction.  The  proof  of  the  new  inverse  scattering  theorem 
is  given  in  Section  6  and  concluding  remarks  appear  in 
Section  7. 


2.  DESCRIPTION  OF  THE  SYSTEM 

The  system  is  composed  of  two  antenna  arrays,  one  for 
transmission  and  another  for  reception,  submerged  in  a  cy¬ 
lindrical  water  container  of  about  3  feet  in  diameter  and  3 
feet  high.  Both  antenna  arrays  are  of  hexagonal  shape.  The 
elements  are  placed  in  a  brick-staggered  arrangement,  cor¬ 
responding  to  a  planar  lattice  with  one  lattice  vector  at  59° 
from  another  and  0.97  times  the  length  (Fig.  la).  The  re¬ 


ceiving  array  is  composed  of  127  elements  whereas  the 
transmitting  array  contains  151  elements.  Each  element  is 
a  short,  water-filled  waveguide.  The  cross  section,  4  mm  X 
7  mm,  is  that  of  a  degenerate  ellipse  (see  Fig.  lb).  The  feed 
structure  consists  of  waveguide- tu-coaxial  adapter  with  an 
insulated  end  feed  which  is  shorted  to  the  broad  wall  of  the 
waveguide  as  shown  in  Fig.  Ic.  This  element  differs  from  the 
one  described  by  Foti  et  al  (elsewhere  in  this  volume)  in  that 
it  is  more  amenable  to  series  production  in  a  monobloc  array 
by  numerically  controlled  milling  machines.  It  is  designed 
for  fixed  tuning  (VSWR  <  1.5)  over  a  iGHz  band  centered 
on  3GHz. 

At  an  axial  distance  of  5  cm  or  farther,  the  underwater 
field  pattern  of  each  element  in  the  forward  direction  is 
similar  to  that  of  a  dipole.  Both  transmitting  and  receiving 
antennas  are  mounted  in  adjustable  frames,  facing  each 
other  for  forward  scattering  imagery.  The  target  is  to  be 
placed  between  the  two  antennas.  The  axial  distance  be¬ 
tween  the  antennas  and  the  target  may  be-adjusted  from  as 
close  as  5  cm  to  a  distance  of  about  35  cm.  Other  engineering 
details  on  water  coupled  antennas  for  medical  microwave 
imagery  may  be  found  in  Ref.  [3]. 

In  order  to  compensate  for  the  differences  of  the  distance 
from  each  of  the  elements  to  the  target,  a  method  of  phase 
and  amplitude  conjugations  is  used.  That  is,  a  factor  which 
includes  both  phase  and  amplitude  is  applied  to  each  ele¬ 
ment  depending  on  its  distance  to  each  focal  point  (see 
Section  4,  cf.  Elq.  (25)].  Instead  of  applying  the  conjugations 
to  the  transmitting  array,  which  would  require  RF  attenu¬ 
ators  and  phase  shifters,  the  conjugations  are  applied  to  the 
receiving  elements  in  the  form  of  off-line  data  processing. 


Fig.  Ib.  Subarray  In  a  stainless  steel  monoblot 
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i.e..  by  multiplying  the  received  complex  field  amplitude  for 
each  element  by  a  complex  factor  that  corresponds  to  the 
conjugation  of  the  phase  and  amplitude  of  the  scattering 
parameter  S21  [17|  measured  for  each  element  illuminated 
by  an  omnidirectional  source  (this  will  be  described  in  detail 
in  Section  5).  On  the  other  hand,  the  phasing  of  the  emitting 
array  is  designed  to  produce  a  near  plane-wave.  Alterna¬ 
tively,  sequentially  overlapping  subarrays  may  be  energized 
to  provide  illumination  of  selected  areas  of  the  target. 

In  order  to  describe  the  application  of  this  technique  to 
a  water  coupled  microwave  imaging  system,  a  brief  digres¬ 
sion  into  the  design  of  a  multiplex  receiver  is  necessary.  The 
receiver  consists  of  127  open  ended  waveguide  elements  with 
127  coax  lines  are  routed  to  23  6P1T  diode  switches  in  a  3- 
tier  reverse  corporate  power  divider  network.  This  network 
provides  the  switching  to  connect  each  receiver  array  ele¬ 
ment  to  a  harmonic  converter.  Two  low  noise  amplifiers  of 
20  dB  gain  are  inserted  under  operator  selection  prior  to  the 
harmonic  converter  to  compensate  for  path  losses  through 
the  coupling  medium.  The  local  oscillator  for  the  harmonic 
converter  is  derived  via  a  directional  coupler  from  a  digital 
synthesizer  which  serves  as  the  signal  source  for  the  trans¬ 
mitting  array.  The  RF  port  of  the  harmonic  converter  is  at¬ 
tached  to  each  element  of  the  receiver  array  via  the  switch 
matrix.  The  IF  port  returns  the  down -converted  receiver 
signal  to  a  complex  ratiometer  which  compares  another 
sample  of  the  transmitted  signal  with  the  IF  signal  from  the 
receiver. 

The  amplitude  and  phase  conjugation,  therefore,  includes 
not  only  element-to-element  variations  in  the  array  geom¬ 


etry,  but  also  path  length  and  insertion  loss  variations  in  the 
switch  matrix  enroute  to  the  harmonic  converter.  The 
complex  ratiometer  provides  this  measurement  by  esti¬ 
mating  B21  for  each  element  over  the  entire  RF  signal  path 
from  the  source  to  the  RF  port  of  the  harmonic  converter. 

The  source  used  for  the  S21  measurement  should  ideally 


Fig.  Id.  A  water-immersed  lens  antenna  for  receiver  calibrations 
is  shown.  This  derived  the  measured  data  for  phase  and  amplitude 
roniugations  to  accomplish  three-dimensional  focusing. 
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Fig.  le.  Pattern  cuts  from  the  dielectric  lens  calibration  antenna. 


consist  of  a  2ir  steradian  omnidirectional  radiator.  One  re¬ 
alization  of  such  a  calibration  source  is  shown  in.  Fig.  Id. 
This  design  is  based  upon  a  dielectric  lens.  It  provides  a  3-dB 
beam  width  of  130®  in  azimuth  and  in  elevation  as  shown  in 
Fig.  le.  The  measured  amplitude  for  each  element,  inclusive 
of  its  path  to  the  harmonic  converter,  for  a  given  ptosition  of 
the  calibration  source  provides  the  amplitude  taper  needed 
to  compensate  for  path  losses  through  the  coupling  medium 
and  insertion  losses  in  the  switch  matrix.  That  is,  the  needed 
amplitude  taper  is  the  inverse  of  the  measured  amplitude 
taper  for  that  position  of  the  calibration  source.  The  needed 
phase  taper  to  insure  coherent  addition  at  each  position  of 
the  calibration  source  is  the  conjugate  of  the  measured  phase 
taper.  All  focal  spot  positions  are  provided  by  translation  of 
the  calibration  source. 

The  calibration  data  set  derived  as  described  above  is  then 
applied  to  the  Ssi  measurements  in  the  presence  of  a  target 
for  each  element  fm  each  focal  point  In  this  way,  the  forward 
scattered  fields  are  scanned  by  a  sensitive  volume  element. 
The  receiver  array  focusing  takes  place  off-line.  The  array 
is  focused  only  in  the  coupling  medium,  not  in  the  target. 
The  forward  scattered  fields  from  the  target  are  differenced 
from  the  beam  pattern  of  the  illuminator  recorded  at  the 
same  plane  in  the  absence  of  the  target. 


3  THE  LOCAL  FIELD  ANALYSIS 

In  this  section,  the  formula  for  the  electromagnetic  field  in 


the  neighborhood  of  a  lattice  radiation  source  is  derived.  The 
formula  will  be  used  to  calculate  the  field  pattern  of  the 
antenna  array  and  the  beam  characteristics.  The  antenna 
array  is  treated  as  a  localized  distribution  of  charge  and 
current  in  a  lattice  structure.  The  following  three  assump¬ 
tions  are  made; 

Assumption  1:  If  df  and  di  represent,  respectively,  the  Ion 
gitudinal  and  transverse  dimensions  of  each  array  element 
with  respect  to  the  direction  of  the  point  of  observation,  and 
In  the  distance  from  the  nth  element  to  the  observation 
point,  it  is  assumed  that,  for  every  element,  the  following 
magnitude  comparisons  are  valid: 

rn  »  dr  and  di,  X  »  dr.  X  »  d-,/r„  ( 1 1 

where  X  is  the  wave  length  of  the  microwave  signal  in  water. 
Under  this  assumption,  the  field  due  to  each  radiating  ele¬ 
ment  may  be  approached  by  the  dipole  approximation. 

Assumption  2:  Mutual  couplings  between  the  radiating  el¬ 
ements  are  included  in  the  local  field  formula  to  the  extent 
that  the  effect  of  ail  coupling  is  assumed  to  he  identical  in 
every  element.  In  other  words,  the  difference  between  the 
peripheral  elements  and  the  interior  elements  with  regard 
to  the  effect  of  mutual  coupling  is  assumed  to  be  negli 
gible. 

Assumption  3:  Mutual  coupling  between  the  radiating  ele 
ments  is  linear  with  respect  to  the  phase  and  amplitude  of 
the  power  input  to  the  elements.  Measured  mutual  coupling 
using  a  Hewlett  Packard  8542C  Automated  Network  ,-\n  i 
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lyzer  proved  this  assumption  to  be  true  and  demonstrated 
that  mutual  coupling  is  less  than  -.lO  dB. 

It  is  remarked  that  Assumption  1  is  of  a  quantitative  na¬ 
ture,  in  the  sense  that  its  degree  of  satisfaction  depends  on 
the  degree  of  quantitative  precision  needed  for  the  field 
pattern.  Although  the  broad  dimension  of  each  array  ele¬ 
ment  is  about  half  a  wavelength  in  the  coupling  medium,  it 
is  the  actual  current  distribution  that  determines  the  size 
of  the  source.  Since  the  dominant  mode  (TEio)  of  the  electric 
field  in  the  aperture  with  respect  to  boresight  angle  is  known 
to  be  sinusoidal,  the  effective  size  of  the  current  distribution 
is  shorter  than  the  broad  dimension  of  the  guide  due  to 
center  weighting.  Our  calibration  measurement  shows  that 
the  field  of  each  individual  array  element  resembles  that  of 
a  dipole,  indicating  that  this  assumption  is  valid  for  the 
system.  It  will  be  made  clear  where  this  assumption,  as  well 
as  the  other  two  assumptions,  enter  into  the  derivation,  so 
that  the  percentage  error  of  the  derived  quantities  may  be 
determined. 

To  derive  the  local  field  formula,  consider  a  localized 
charge  density,  f>.  and  current  density,  •!,  distributed  in  a 
space  region  For  an  antenna  in  a  free  space,  V'  indicates 
the  space  occupied  by  the  antenna  array,  as  well  as  its  ac¬ 
cessories.  Without  losing  any  generality,  monochromatic 
time-variation  is  assumed,  so  that 

ptx.tl  =  p(x)e”'*' 

and  (2) 

-/(x.tl  =  ./(xle"'*’. 

Accordingly,  all  other  field  quantities  resulting  from  p  and 
also  vary  with  time  monochromatically.  Any  other  time- 
variation  can  always  be  obtained  by  superposition  of 
monochromatic  waves.  From  Maxwell's  equation,  the  vector 
potential  at  any  point  outside  of  V  is  given  by,  in  the 
Gaussian  system  of  units  [18|, 

7  r  1  ,,  exptiklx  -  x'l)  , 

A(x)  =  I  -J(r) - ; - ; - dx  .  (3) 

Jvf  |x-x| 

from  which  one  obtains  the  magnetic  induction  B  and  elec¬ 
tric  field  strength  E: 


B(x)  =  VXA(x)  (4) 

and 

6(x)  =  p  V  X  B(x),  (5) 

k 

where  k  is  the  magnitude  of  the  wave  vector  in  the  medium. 
For  water,  k  is  a  complex  quantity  ki  -t-  ikz  where  k|  is  equal 
to  2x  times  the  inverse  of  the  wave  length  and  k2  is  the  in¬ 
verse  of  the  distance  over  which  the  field  is  attenuated  by 
a  factor  of  e  =  2.72  (equivalent  to  a  power  loss  of  8.7  dB).  At 
an  operating  frequency  of  3  GHz  in  water,  the  values  of  k| 
and  k2  are 

kj  =5.5  cm"'  and  kz  =  0.44  cm"'.  (6) 


The  current-charge  volume,  V,  is  divided  into  a  number 
of  subvolumes,  each  denoted  by  V„,  which  represents  the 
space  occupied  by  the  nth  radiating  element.  Let  Xn  be  the 
center  of  Vn  and  denote  by  Jn  and  p„,  respectively,  the  cur¬ 
rent  density  and  the  charge  density  in  Vn  with  respect  to  its 
center,  then 

Jn(x)  =  +  Xp)  and  Pn(x)  =  p(x  +  Xn).  (7) 


Clearly,  Equation  3  may  also  be  written  as: 


N  ^  1 

A(X)  =  L'  - 

n=  1  Vn  C 


1  -  exp(ik|x  -  x' 

J(x)  — [7 — 

X  -  X 


■dx'. 


(8) 


Making  the  change  of  variable  x'  -  Xn  -t  x  .  where  the  new 
x'  is  a  vector  from  the  center  of  each  element  to  the  volume 
dx',  which  is  identical  for  every  n  and  thus  may  simply  be 
regarded  as  a  vector  in  V,.  using  Eq.  7.  Eq.  8  becomes 


A(x 


n=  1  '  l 


1  -  exp  ik] (x  -  Xn)  ~  x' 
"  Jn(x') 77 77 — 

C  (X  -  Xn>  -  X 


dx'.  i9i 


Note  that  each  integral  in  the  right  hand  side  of  Eq.  9  is  the 
same  as  that  in  Eq.  3.  except  that  x  is  replaced  by  x  -  Xn  and 
the  space  of  integration  is  over  only  the  center  element,  V,. 
instead  of  the  entire  array,  V.  Therefore,  even  though  the 
observation  point  is  in  the  neighborhood  of  the  array,  as  long 
as  |x  -  Xnl  =  rn  is  much  greater  than  the  size  of  each  array 
element,  d.  one  may  expand  the  integrand  in  Eq,  9  in  powers 
of  x'/rp  (note  that  |x'|  ^  d). 

So,  denote  by  x/  and  Xt',  respectively,  the  longitudinal  and 
transverse  components  of  x'.  i.e..  the  projections  of  x'  in  the 
directions  parallel  and  perpendicular  to  x  -  Xp.  then 


Xc'  =  "^  |x'  •  (X  -  .Xp)|(x  -  .Xp) 


1 10) 


and 


Define 


Xt'  = - -  |x'  X  (X  -  Xp)|  X  (X  -  Xp),  111) 

rp^ 


=  |Xf  X,  =  |x, 


(12) 


Then  one  has  the  identity 

|x  -  Xp  -  x'l  =  [(Tp  -  Xe')'i  -)-  Xt'^l'^'^ 
or  (13) 


^  ^  (xt'/rp)'^  |i-"2 

(l-x//rn)2 


Under  Assumption  1,  Xt'/tp  and  xo'/rp  are  both  small  quan¬ 
tities.  So  expand  |x  -  Xp  -  x'|  and  oi,l;“r  functions  of  it  in 
powers  of  n-e'hn  and  Xt'/rp,  then  the  following  series  result: 
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and 


As  to  the  factor  expiik I X  -  Xn  -  x'| ),  its  expansion  depends 
not  only  on  the  relative  magnitude  of  rn  and  d,  but  also  on 
the  magnitudes  of  krn  and  kd.  If  the  real  part  of  kd  is  small, 
which  is  valid  under  Assumption  1,  then,  except  for  the  first 
term,  every  term  on  the  right  hand  side  of  Eq.  14  is  much 
smaller  than  1.  The  exponential  of  the  series  then  gives 


exp(ik|x  -  Xn  -  x'|)  =  explikrnXl  -  ikx/>'  +  ...).  (16) 


Combining  this  expansion  with  Eq,  1,'i.  the  integrand  of  Eq. 
9  becomes 

1  exp( ik| X  -  Xp  -  x'l  I 

-  .I„lx  l - 


|x  -  Xn  -  x' 


1  exp(ikrn)  ,  ,,,, 

= - Jn(x  ) 

c  fn 


+ 


—  -  ikxe'  +  .  .  . 


(17) 


Substituting  the  leading  term  on  the  right  hand  side  of  Eq. 
17  into  Eq.  9,  an  approximation  for  the  vector  potential  is 
obtained  as 


A(x) 


^  I  expdktn) 

n=  1  C  Tp 


X, 


.]n(x')dx'. 


(18) 


It  can  be  shown  that  the  integral  in  the  right  hand  side  is 
proportional  to  the  total  dipole  moment,  Pn,  of  the  array 
element  Vp: 


i 


X, 


jn(x  )dx' 


i 


J(x')dx', 


(19) 


where  the  dipole,  Pn.  is  defined  as  the  moment  of  the  charge 
distribution  of  the  nth  radiating  element  with  respect  to  its 
center: 


Pn  = 


(20) 


Therefore,  the  vector  potential  may  also  be  written  as 


.A(x) 


N 

-ik  ?■ 


Pn 


exp(-ik|x  -  Xn| ) 


(21) 


which  is  the  field  due  to  N  radiating  dipoles. 

Assumptions  2  and  3  are  now  applied  to  Eq,  21.  Noting 
that  all  elements  have  the  same  geometry,  the  only  factors 
that  could  contribute  to  different  values  of  Pn  for  different 
elements  are  the  input  power  and  phase  and  the  differences 
in  the  current-charge  distributions  due  to  mutual  coupling. 
Under  Assumption  2,  the  last  factor  is  assumed  to  be  negli¬ 
gible,  and,  under  Assumption  3,  Pn  must  be  proportional  to 
the  input  phase  and  amplitude  factor.  Therefore  pInCn  may 
be  substituted  for  Pn  in  Eq.  21 ,  where  p  is  the  dipole  moment 
for  each  radiating  element  at  a  certain  standard  input.  is 
the  illumination  factor  for  the  nth  element,  and  Cn  is  a 
complex  factor  representing  the  phase  and  amplitude  con¬ 
jugations.  In  is  used  as  a  controlling  factor  to  modify  the 
mainbeam  shape.  Equation  21  then  becomes 


N 

A(x)  =  -ikp  X 

n=  1 


(22) 


The  electric  and  magnetic  fields  may  be  obtained  from  the 
above  equation  by  applying  Eqs,  4  and  o  on  A(x).  Again 
applying  Assumption  1,  the  results  are 


-  ,,  N  X  -  Xn  exptiklx  -  Xnl  ) 

B(x)  =  -k-^pX  i:  InCnTT - Xt  |.  .  I  (23) 

n=l  |X-Xni  |x-Xn| 


and 


_  N  /  X  -  Xr.  \ 

E(x)  =  -k2  InCn  px— ^ 
n=l  I  |x-Xn|/ 

X  X  ~  Xn  exptik|x  -  Xn|  ) 

|x-Xnl  |x-Xn| 


(24) 


The  definition  of  the  quantities  in  Equations  22  through 
24  are  summarized  below" 


A  =  the  vector  potential,  in  the  Gaussian  system  ol 
units 

B  =  the  magnetic  induction,  in  the  Gaussian  system  of 
units 

Cn  =  the  complex  number  representing  the  phase  and 
amplitude  conjugation  for  the  nth  radiating  ele¬ 
ment  (see  next  section) 

E  =  the  electric  field  intensity,  in  the  Gaussian  system 
of  units 

1,,  =  the  illumination  factor  for  the  nth  element:  this 
factor  is  used  to  control  the  beam  shape  (see  next 
section) 

k  =  the  complex  number  representing  the  magnitude  of 
the  wave  vector  of  the  radiation  in  water;  the 
values  of  its  real  and  imaginary  parts  for  a  3  GHz 
radiatii)n  are  given  in  Equation  6 

n  =  a  subscript  denoting  the  nth  radiating  element 

N  =  the  total  number  of  radiating  elements  in  the 
array 

p  =  the  dipole  moment  of  each  radiating  element  at  a 
standard  phase  and  amplitude  input  (i  e.,  for  1„ 
and  ('n  being  unity),  in  the  Gaussian  system  ol 
units 

X  =  the  vector  representing  th"  observat  jou  point  with 
respect  to  the  center  of  the  arra> 

X|,  =  the  vector  representing  the  center  of  the  nth  ra 
diating  element  with  respect  to  the  center  of  the 
array. 


4.  PHASE  AND  AMPLITUDE 

CONJUGATION^  AND  THREE- 
DIMENSIONAL  .'OCUSING 

The  field  patterns  presented  in  Eqs.  (22  2.'))  depend  on  the 
set  of  factors  lInCnI.  In  this  section  we  introduce  a  formula 
for  assigning  the  values  of  Cn  to  provide  a  maximum  relative 
field  at  a  desirable  focal  point.  The  idea  is  that,  if  one  wishes 
to  focus  the  field  of  the  arrav  at  a  point,  sav,  X|,  one  lan 
maximiz.e  the  field  at  that  (toint  by  aiiplymg  a  phase  and 
amplitude  taper  which  compensates  the  |)ropagation  loss 
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Fig.  2a-2c.  Relative  amplitude  of  the  vector  field  in  the  surface 
7.  =  7  cm  with  phase  and  amplitude  conjugations  focused  at  if  = 
(0. 0,  7)  cm.  The  element  spacings  are  (al  AX  =  0.381  cm  and  AY 
=  0.6,35  cm,  (b)  AX  =  0.762  cm  and  AY  =  1.270  cm.  (c)  AX  = 
0.9SO  cm  and  AY  =  1.6.50  cm.  The  3  dB  widths  of  the  main-beam 
in  both  X-  and  y-directions  are  as  indicated. 


Figs.  2d-2e.  Relative  amplitude  of  the  vector  field  in  the  surface 
7.  =  7  cm  with  phase  and  amplitude  conjugations  focused  at  Xf  = 
(4,  3. 7)  cm.  The  element  spacings  are  (dl  AX  =  0.762  cm  and  AY 
=  1.270  cm,  (e)  AX  =  0.990  cm  and  AY  =  1.650  cm.  The  3  dB 
widths  of  the  main-beam  in  both  x-  and  y-directions  are  as  indi¬ 
cated. 


and  phase  differential  from  each  array  element  to  the  focal 
point.  Thus,  to  focus  the  main  beam  at  the  point  Xf,  the  factor 
Cn  is  assigned  as 

Cn  =  |Xf-  Xn|exp(-ik|xf-  Xn|). 


Noting  that  k  is  a  complex  number,  the  exponential  factor 
in  the  above  equation  includes  a  phase  factor  and  an  am¬ 
plitude  factor  to  compensate  for  the  absorption  by  the  cou¬ 
pling  medium.  For  this  reason,  we  call  Cn  the  phase-ampli- 


(25) 
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BW  =  (0  70,  '  821  cm 

=  0  990 


Figs,  3a-3e.  Relative  amplitude  of  the  vector  field  in  the  yz-plane  those  in  Figs.  2a-2e.  respectively.  The  ,3  dB  widths  of  the  main 

at  X  =  0,  with  focusing  and  element  spacings  corresponding  to  beam  in  both  cm  y-  and  z-directions  are  as  indicated. 


tude  conjugation  factor.  The  factor  !„  in  Eqs.  (22-24),  which 
is  called  the  illumination  factor,  may  be  used  as  additional 
leverage  to  provide  optimal  resolution  and  minimize  the  side 
lobes. 

For  a  planar  array  of  N  X  M  radiating  dipole  elements, 
Eqs.  22  and  23  may  be  written  as 


e'^lx  -  Xnml 
|x  -  v,,ml 


and 

Cnm  =  [Xf-  X„n,|exp(-ik|Xf  -  Xnml).  (27) 

As  we  shall  see  below,  the  phase-amplitude  conjugation 
as  presented  in  Eq.  (2,5)  does  not  make  the  field  peak  at  ex¬ 
actly  the  point  Xf.  This  is  simply  due  to  the  fact  that  we  are 
in  the  local  field  region.  For  the  array  structure  described  in 


(26) 
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X  t  c  m )  X  ( V  • 
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X  (cm) 


JX  '  0  990 
-  1  66C 
Xf  -  4  0 
Yt  *  30 


X  (cm) 


Figs.  4a-4e.  Amplitude  contours  of  the  vector  field  in  isodecibels  corresponding  to  the  ca'-es  indicated  in  Figs  2a-2e  respectively. 


Section  2,  the  main  lobe  beamwidth  is  minimized  when 
is  taken  to  be  a  uniform  illumination. 

Now  we  present  some  result  of  the  field  patterns  using  the 
phase  and  amplitude  conjugations.  Only  the  absolute  value 
of  vector  potential  and  the  corresponding  field  character¬ 
istics  are  presented  here.  The  structure  of  the  array  lattice 
is  illustrated  in  Fig.  la,  and  the  equations  used  are  (26),  (27) 


and  (6),  with  illumination  factor  Inm  =  1,  and  polarization 
taken  to  be  in  y-direction.  In  all  the  figures,  the  plane  of  the 
array  is  taken  to  be  the  xy-plane  and  the  z-axis  is  perpen¬ 
dicular  to  the  array  plane  and  pointing  to  the  forward  di¬ 
rection.  In  the  following  discussion  and  in  all  figures,  the 
phrase  “mainbeam"  is  used  in  reference  to  field  character¬ 
istics  of  the  3-dimensional  focal  region.  Figures  2(a)  through 
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Y  <cm) 


Y  (cm) 


Y  (cm)  Y  (cm) 

Firs.  5a-5e.  Amplitude  contours  of  the  vector  field  In  isodecibelscorresponding  to  the  rases  indicated  in  Fit's  :ia  ie  ri-spect  ivelv 


2(e)  show  various  field  patterns  with  different  interelement 
spacings  in  the  transverse  plane  at  an  axial  distance  of  7  cm, 
which  is  also  the  axial  distance  of  the  desired  focal  point. 
Similarly,  Figs.  3(a)  through  3(e)  are  the  longitudinal  field 
patterns  at  different  interelement  spacings  and  focal  points. 
The  corresponding  3-dB  full  mainbeam  width  is  given  on 
each  figure.  Note  that  Figs.  2(b),2(c)  and  3(b),3(c)  are  for  the 
array  with  lattice  structure  shown  in  Fig.  la.  Figures  2(d)- 
,2(e)  are  for  arrays  that  are  similar  to  the  one  shown  in  Fig. 
la,  except  with  different  lattice  spacings.  Figures  2(a)  and 


3(a)  are  for  the  same  size  array  as  shown  in  Fig.  la  but  with 
419  instead  of  127  elements.  Figures  4  and  5  are  the  corre¬ 
sponding  pictures  of  Figs.  2  and  3,  and  plotted  in  isodecibel 
field  contours.  For  easy  comparison,  the  beam  characteristics 
of  these  figures  are  tabulated  in  Table  1.  These  data  show 
that  larger  interelement  spacings  result  in  narrower  3-dB 
beamwidths  and  more  accurate  focusing;  however,  larger 
interelement  spacings  also  bring  grating  lobes  closer  to  the 
mainbeam.  This  result  demonstrates  that,  with  respect  to 
beamwidth  reduction,  the  array  size  plays  a  more  important 
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Table  1 


ARRAY  BEAM  CHARACTERISTICS  OF  THE  FIELD 
PATTERNS  AT  3  GHz 


FOCAL  POINT 

(x  V  z)  cm 

element  spacing 

(dx  dy)  cm 

PEAK  AT  Z-7 

(x  y  z)  cm 

PEAK  AT  X=X, 

(x  y  z)  cm 

BEAMWIDTH  AT  Z=7 

(dx  dy)  cm 

BEAMWIDTH  AT  X=Xf 

(dy  dz)  cm 

0 

0 

7 

.381 

.635 

0 

0 

7 

0  0  6 

.60 

.60 

.81 

2.83 

0 

0 

7 

.762 

1.27 

0 

0 

7 

0  0  6 

.60 

.60 

.50 

1.60 

0 

0 

7 

.990 

1.65 

0 

0 

7 

0  0  6.5 

.53 

.53 

.20 

1.17 

4 

3 

7 

.762 

1.27 

4 

3 

7 

4  2.75  6.25 

.74 

.67 

.80 

1.87 

4 

3 

7 

.900 

1.65 

4 

3 

7 

4  3  6.5 

.60 

.60 

_ 1 

.70 

1.82 

role  than  the  number  of  elements  The  resolutions  of  '/-jK  in 
transverse  direction  and  IX  in  longitudinal  direction  can  be 
achieved  with  the  use  of  phase  and  amplitude  conjugation. 
Note  also  that  the  peaks  along  z  =  2.5  cm  shown  in  Figs. 
d(a)-3(e)  and  Figs.  5(a)-5(e)  are  due  to  the  single  element 
that  is  closest  to  each  of  these  peaks.  As  the  distance  becomes 
.so  close  to  an  individual  element,  the  coherent  addition  from 
other  elements  is  negligible  in  comparison.  As  long  as  the 
target  is  not  closer  than  3  cm  from  the  array,  these  peaks  will 
not  cause  any  problem  for  actual  applications  in  microwave 
imagery. 


.5.  THE  INVERSE  SCATTERING 

THEOREM  AND  ITS  APPLICATION 
TO  IMAGE  RECONSTRUCTION 


susceptibility  of  the  entire  system,  including  the  homoge¬ 
neous  medium  and  the  target,  so  that  x  ~  Xm  is  null  outside 
the  target.  It  is  assumed  that  there  is  no  free  charge  or  cur¬ 
rent  distribution  (including  ionic  charge  and  current)  in  the 
target  or  the  medium,  and  that  both  have  the  same  homo¬ 
geneous  magnetic  permeability  however,  we  remark  that 
the  theorem  may  be  generalized  to  include  free  charge-cur¬ 
rent  and  magnetization.  Let  the  target  be  illuminated  by  a 
plane  wave  Ejnc  of  frequency  ii)/2t,  which  induces  electric 
polarization  P  in  the  target  and  produces  a  scattered  field 
R»catt-  Let  Jw  be  some  weighing  function  and  Aw  be  the  cor¬ 
responding  field  derivable  from  Jw  in  the  same  way  as  the 
■ector  potential  is  derivable  from  a  current  density.  More 
precisely, 

(v^ -I- k„2)Aw  =  -  Jw,  (28) 

c 


The  objective  of  all  inverse  scattering  problems  is  to  recon¬ 
struct  the  target  from  the  scattered  field,  k’rom  Maxwell’s 
electromagnetic  theory,  if  one  knows  the  field  everywhere 
in  space,  the  polarization  charge-current  distribution  of  the 
siattering  source  can  be  derived  completely.  However,  in 
practice,  one  can  only  measure  the  scattered  field  at  a  limited 
number  of  points  in  space  which  are  often  confined  in  a  small 
region.  The  question  is  then  how  much  information  on  the 
scattering  target  one  can  infer  based  on  a  limited  knowledge 
of  the  .scattered  field.  Here  we  present  a  theorem  |8|.  which 
is  indeed  a  generalization  of  the  lyorentz  reciprocity  theorem 
|19j-|2l|,  and  show  that  it  can  provide  a  good  facility  in 
answering  the  above  question. 

In  this  section  we  shall  state  the  theorem  and  describe  its 
application  to  microwave  biological  imaging,  and  leave  the 
details  of  the  prmif  to  the  next  section.  Consider  a  dielectric 
target  immerged  in  a  homogeneous  and  dissipative  medium 
of  dielectric  constant  fm  and  dielectric  susceptibility  Xm 
(both  are  complex  numliers).  l,el  xlxl  descriln- the  dielectric 


or  its  reverse  equivalence, 

uw,  ^ 

Aw(x)  =  — J  - - rJw(x')dx'.  (29) 

c  |x-x'| 

Then  the  inverse  scattering  theorem  may  lie  stated  as 
below; 


- Vp, 

dk  =  ///  Jw  •  Ewralldx  (30) 

where  k^,  defined  as  vVni<in‘*t/c,  is  the  complex  wave  num¬ 
ber  in  the  homogeneous  medium,  and  p,  and  Jg  represent  the 
equivalent  charge  density  and  its  time  derivative,  the  current 
density,  due  to  the  polarization  of  the  target  in  excess  of  the 
homogeneous  background  polarization: 


JJ/A 


w(x)  • 


_  X  Xm  A  I  X  Xm  dP 

p.,  =  -  V  - - P,  ,L  = - — 

X  X  dt 


(31) 
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As  we  shall  see  in  the  next  section,  ik  and  .1,  are  indeed  the 
source  of  the  scattered  field. 

This  theorem  is  more  general  than  the  Lorentz  recipriK'ity 
theorem  (191  |‘J1 1  since  J,,,  as  well  as  the  associated  vector 
field  A»,  is  only  a  weighing  fum  tion.  They  need  not  be 
physical  quantities;  for  example,  .1^^  needs  not  satisfy  the 
equation  of  continuity  and  Aj^,  needs  not  satisfy  the  Lorentz 
gauge  condition  (or  any  alternative  gauge).  The  proof  of  this 
theorem  is  based  on  the  fact  that  the  medium  is  dissipative 
and  therefore  a  Hilltert  space  may  l>e  defined  in  which  all  the 
differential  operators  involved  may  have  their  Hermitian 
conjugations  defined.  The  great  facility  of  this  theorem  in 
applications  to  image  reconstuction  lies  on  the  flexibility  of 
chixcsing  the  weighing  function  -]„  If  the  scattered  field 
is  measured  at  a  set  of  spatial  [joints  lx„l,  which  are  the  lattice 
points  of  our  receiving  array,  and  if,  say,  only  the  y-polar 
ization  is  measured,  then  the  weighing  function  .l»(.x)  may 
be  chosen  to  be  only  in  the  y  direction  and  to  be  a  discrete 
distribution  over  the  points  Ix^I.  It  then  follows  from  Kq.  (29) 
that  the  resulting  weighing  vector  potential  Av,(x)  will  also 
be  in  the  y-direction.  One  may  further  adjust  the  phases  and 
amplitude.^  of  -iwIXnl  to  optimize  the  weighing  potential 
A„(x)  at  any  desired  point  of  the  target.  It  was  shown  in  the 
previous  section  that  applying  the  phase  and  amplitude 
conjugations  to  a  transmitting  array  could  provide  a  11- 
dimensionally  focused  radiation  in  the  neighborhood  of  the 
array.  For  imagery  application  with  an  incident  plane  wave, 
this  technique  may  he  applied  to  the  receiving  array  to 
"focus"  the  weighing  potential.  Thus,  by  setting  .i«(Xn)  to 
be  the  phase  and  amplitude  conjugation  factor,  one  may 
make  Aw(x)  negligibly  small  inside  the  target  organ,  except 
for  a  sharp  peak  near  any  desired  point,  say  X(,  which  we  shall 
call  the  focal  point.  Note  that,  since  vanishes  outside  the 
target  region  and  so  the  integration  in  the  left  hand  side  of 
Fq.  (30)  is  limited  to  only  the  target  region,  A„(x)  may  as¬ 
sume  any  value,  however  large,  outside  the  target  region. 
Then  the  right  hand  side  of  Kq.  (30)  is  obtainable  from  the 
measured  data,  whereas  the  integration  on  the  left  hand  side 
is  dominated  by  the  integrand  at  the  focal  point. 

It  is  emphasized  that  the  "fiK'Using"  described  above  does 
not  involve  active  focusing  of  the  transmitted  microwave 
energy  on  any  point  of  the  target.  Hather,  it  is  sinqilv  a 
mathematical  management  of  the  measured  .scattering  data 
so  that  the  retrieved  information  of  the  dielectric  properly 
may  be  ‘‘focused"  on  a  desired  point  inside  the  target  In  this 
sen.se.  it  may  be  considered  as  a  focusing  of  the  receiving 
array  Since  one  microwave  exposure  of  the  target  will  [iro 
vide  N  samples  of  the  complex  scattered  field,  where  N  is  the 
number  of  receiving  element-s,  it  is  then  theoretically  possible 
to  retrieve  N  estimates  of  the  spatial  distribution  of  ibe 
complex  dielectric  properties  of  the  target  by  o[)timai 
management  of  the  measured  data.  The  inveise  scattering 
theorem  described  above  is  to  provide  a  basis  for  such  pur 
pose.  VVe  shall  discuss  in  Section  7  the  maximum  limit  of 
retrievable  information  from  a  single  exposure  of  microwave 
radiation. 

To  express  this  imagery  application  more  clearly,  we  set 
J»(x)  =  2lni)(x  -  Xn)ydn  and  substitute  Fxj.  (31)  lor  /).  and  T 


in  Kq.  (30),  then  replace  P-'l’/Ot-  by  -u.'  !'  and  [lerlnnn 
partial  integration  on  the  first  term  on  the  left  hand  side. 
Noting  that  A^  vanishes  at  infinity,  one  then  gels 


<mMn 


JM'. 


r  ,  X  X  rn  r-. 
-(7  ■  A„)  V  ■ - 1 


,1  '.’A  k  ^ XI  rs 


dx  =  V  ,|„K,. 


With  the  discrete  distribution  ol  -Llxl.  Kq. 
to 


,.,(x,,l 

129 1  reduces 


.■\„(x)  =  V  — exp(ik,„|x  -  \„1  I 
n  |x  -  X„| 

w  hich  is  similar  to  the  Jcctor  [jolenlial  (iroduced  In  a  si  l  ol 
di|)oles  as  given  by  hxj.  ('2'2).  To  obtain  the  dielectrii  ininr 
mation  at  any  given  point  X|  inside  the  target  from  the 
measured  .scattered  field  at  the  (joints  |x„l.  one  sets  .)„  in  the 
right  hand  side  of  P'q.  (:V2l  to  be  the  phase-amplil iide  con 
jugation  factor: 

•  IntXf)  =  (c/pmll  Xn  ~  Xfl  ■  exp(  -ik|  X„  -  X,|  I  13)1 

Based  on  the  [jhase-amplitiide  conjugation  discussed  in  the 
[jrevious  .section,  the  resulting  .Aw(-x)  has  a  sharp  peak  at 
and  is  otherwise  negligible  in  the  target  region.  Owing  to  the 
factor  V  —  Xm  which  vanishes  outside  the  target  region,  all 
sidelobes  of  Aw(x),  however  large,  which  lie  outside  the  target 
region  will  not  contribute  to  the  integral  on  the  left  hand  side 
of  Kq.  (32).  The  main  contribution  of  the  integral  is  the  di¬ 
electric  characteristics  of  the  target  at  the  point  Xf.  which  is 
then  equal  to  the  weighted  version  of  the  measured  values 
-n'InKs,ait(x„).  For  this  reason  we  refer  to  the  point  X|  as  a 
focal  point  since,  by  Kq,  (32),  the  sum  of  the  products  of  the 
conjugation  factors  and  the  measured  scattered  fields. 
“n''nF.,„it(Xn).  "focuses"  the  result  to  the  dielectric  polar¬ 
ization  at  the  point  if  inside  the  target.  By  scanning  the 
vector  X|  in  P3q.  (34)  through  the  target  region,  one  then 
reconstructs  an  image  of  the  target.  The  response  [leak  in 
.Awix)  at  the  "focal  point"  determines  the  s[)atial  resolution 
and  (he  flat  response  of  A^(i)  within  the  target  elsewhere 
contributes  to  ana'ytical  interjiretation  ot  the  imagers  i7| 
The  size  of  the  target  is  limited  by  the  liK'ations  of  the  gral  me 
lobes,  since  these  grating  lobes  will  not  contribute  to  the 
integral  on  the  left  hand  side  of  Kq.  (.‘(21  as  long  as  Ihes  are 
outside  of  the  target  region.  Another  method  to  reduce  the 
grating  lobe  interference  is  to  limit  the  illumination  region 
so  (hat  the  liH'al  field  is  negligible  in  the  grating  lobe  region. 
For  a  larger  target,  this  will  require  dividing  the  target  into 
smaller  illumination  regions 


6.  PROOF  OF  THE  INVERSE 
SCATTERING  THEOREM 

Our  first  approach  in  proving  the  theorem  is  to  separate  ttu' 
target  from  the  surrounding  homogeneous  medium  so  th.it 
the  scattered  field  may  be  considered  as  due  to  a  localized 
charge  current  source.  We  shall  consider  the  entire  n.'leni 
to  be  a  su|>erposition  of  a  distribution  ol  dielectric  susce|> 
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tibility  x(i)  ~  Xm  ®nd  a  homogeneous  medium  of  dielectric 
susceptibility  Xml  thus  x  “  Xm  may  be  considered  as  the 
excess  susceptibility  of  the  system  over  the  homogeneous 
background.  Since  the  medium  is  homogeneous,  the  prop¬ 
agation  of  the  scattered  field  is  then  governed  by  a  set  of 
Maxwell  equations  similar  to  that  in  free  space  except  that 
the  free  space  permitivity  and  permeability  are  replaced  by 
that  of  the  homogeneous  medium.  Using  the  notations  and 
assumptions  that  were  described  at  the  beginning  of  the 
previous  section,  the  Maxwell’s  equations  for  the  total  field 
(6,  B)  may  be  written  c.s 

V-(l  -f4jrx)E  =  0 


Mm  C  dt 


(35) 


V  ■  B  =  0 


A  1 

vxE-f - =  0 

c  dt 

In  order  to  separate  the  source  of  the  scattered  field  from  the 
homogeneous  background,  we  transport  all  quantities  in¬ 
volving  excess  dielectric  susceptibility  to  the  right  hand  sides 
of  the  above  equations,  so  that  the  left  hand  sides  resemble 
that  of  the  Maxwell  equations  in  free  space.  The  results  are, 
with  the  definition  of  p,  and  J,  given  by  Eq.  (31). 

•  6  =  47rp, 

=  (36) 

Mm  C  dt  C 

v-S  =0 

vxetl2?.o 

c  dt 


where  (m  =  1  +  4)rXm  is  the  dielectric  permitivity  of  the 
homogeneous  medium,  and  P  =  xE  is  the  electric  polariza¬ 
tion  in  the  target.  Replacing  the  total  field  (fi.  8)  in  Eq.  (36) 
by  the  sum  of  the  incident  wave  and  the  scattered  field,  since 
the  incident  wave  satisfies  the  homogeneous  (viz.  sourceless) 
version  of  Eq.  (36),  one  finds  that  the  scattered  field  satisfies 
exactly  Eq.  (36).  Denote  the  scattered  field  by  a  subscript 
s,  one  then  has 

'  Ej  —  47rp,j 


—  V  X  B, 

Mm 


<n.  dE,  _  47r  - 
c  dt  c 


(.37) 


V  ■  8,  =  0 


V  X  8,  -t  =  0 

c  dt 


Thus  the  scattered  field  alone  may  be  considered  as  that 
produced  by  the  localized  charge-current  density  (p„  .1,). 

Due  to  the  dissipation  of  the  medium,  the  electromagnetic 
field  vanishes  exponentially  as  |x(  »  <».  Thus  it  is  possible 
to  define  a  Hilbert  space  in  which  the  electromagnetic  fields 
and  charge  densities  are  vectors.  We  define  the  scalar 
product  between  any  pair  of  vectors  (f.  g)  as  the  integral  over 
the  entire  space: 


fff  f*(x)g(x)dx, 

which  exists  owing  to  the  fact  that  f  and  g  diminish  expo¬ 
nentially  as  |x|  -»  ».  With  this  definition,  all  differential 
operators  involved  in  the  Maxwell  equations  as  well  as  ail 
derivative  equations  may  have  their  Hermitian  conjugations 
defined.  It  is  based  on  this  fact  that  the  inverse  scattering 
theorem  can  be  proved.  For  mathematical  simplicity,  we 
shall  prove  the  theorem  using  the  Hilbert  space  notations. 
Noting  that  the  Hilbert  space  product  combined  with  the 
hermiticity  of  linear  operators  is  equivalent  to  integration 
by  parts  involving  differential  operators,  the  proof  may  also 
be  made  equivalently  in  conventional  differential  equation 
form. 

We  shall  now  derive  the  inverse  scattering  theorem  from 
Eq.  (37).  First  note  that  the  third  and  fourth  equations  of 
Eq.  (37)  show  that  the  scattered  electric  field  and  magnetic 
induction  are  derivable  from  a  pair  of  scalar  and  vector  po¬ 
tentials  ((Ds,  As),  which  shall  be  called  the  scattered  poten¬ 
tials: 

B,  =  V  X  A, 

£,  =  -V(J,  -  -  dAs/dt  (38) 

c 

The  relationships  of  Eq.  (38)  still  leave  another  degree  of 
freedom  on  the  choice  of  the  potentials.  We  shall  use  the 
conventional  choice  of  Lorentz  gauge  condition 

v.A.  +  ^^  =  0.  (39) 

c  dt 

Then,  from  the  first  two  equations  of  Eq.  (37),  the  scattered 
potentials  satisfy  the  wave  equation  with  ps  and  .Is  respec¬ 
tively  as  the  source: 

(v'^  +  km2)0j  =  -47rp,/<m 

(v2-fk„^)A,=  MmJs  (40) 

c 

where  kn,^  =  (a)^/c^)Pnifm  and  we  have  replaced  d/dt  by 
-iu). 

To  facilitate  our  proof  of  the  theorem,  we  shall  denote  by 
I  As)  and  |Js)  respectively  the  4-dimensional  potential  and 
the  4-dimensional  current  density: 


Also  denote  by  Kn,  the  4-dimensional  Helmholtz  operator 
associated  with  kn,: 

K„,  =  (v2-tk,„2)[’°j  (41) 

Then  Eq.  (40)  may  be  expressed  as 

^ml  As)  =  I  Jg) 

c 


(42) 
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To  express  the  relationship  between  the  electric  field  and 
the  vector  potential,  we  define  the  4 -dimensional  E- 
vector. 


(43) 


and  the  4-dimensional  S-operator  and  its  hermitian  conju- 
pate. 


|J«,>  and  |Aw>  are  weighing  vectors  which  satisfy  only  the 
wave  equations,  but  not  the  equation  of  continuity  nor  the 
gauge  condition.  So  they  are  not  physical  quantities.  For 
instance,  the  weighing  current  |  J»,)  may  be  a  point  function 
with  only  the  vector  component,  viz.,  a  delta-function  cur 
rent  in  space  without  charge  density.  From  Eq.  (41),  the 
complex  conjugate  of  Km  is  equal  to  its  hermitian  conjugate, 
therefore  the  complex  conjugate  of  Eq.  (48)  gives  (Kn,)^|A„) 
=  -  4ir/c|Jw>,  where  t  denotes  the  hermitian  conjugation. 
The  hermitian  conjugate  of  this  equation  then  gives  ( Aw|  Km 
=  <'Twl  ( -47r/c).  Taking  the  scalar  product  of  both  sides  with 
|F’s>  and  utilizing  Eq.  (46),  it  then  gives 

-  —  <Jw|E.>  =  <A„|Km|E.>  =  -  —  <A.|F.) 

C 

'rherefore 


lE,>=SmlAs)  (4.^) 

Note  that  the  above  equation  actually  represents  two  rela¬ 
tionships:  the  scalar  component  is  equivalent  to  the  gauge 
condition  Eq.  (39).  and  the  vector  condition  gives  the  electric 
field  in  terms  of  the  potentials.  Since  the  operators  Km  and 
Sm  commute  with  each  other,  operating  both  sides  of  Fxj.  (42) 
by  Sm  yields  the  wave  equation  for  IE*): 

Km|E,>  =  -~Sm|J,)  =  -  — IF,)  (46) 

c  c 

which  shows  that  SmIJs)  may  be  considered  as  the  source  of 
( E),)  as  much  as  the  current  density  |  J,)  is  the  source  of  the 
potential  |  A,),  therefore  we  denote  it  by  (F,).  By  the  equa¬ 
tion  of  continuity  V-J,  -  ioip,  =  0,  we  then  have 

|F,)  Sm|3,>  =  ^ 

c  „  Pm  d,l,  I 

1  / 

'  (m  C  dt' 

So  far  we  have  rewritten  all  electrodynamic  equations  for 
a  dissipative  medium  in  their  Hilbert  space  representations. 
To  complete  our  description  of  the  Hilbert  space,  we  define, 
for  any  vector  |F),  a  complex  conjugate  vector  IF)  which 
corresponds  to  the  complex  conjugate  in  the  x-representa 
tion, 

(xlF)  =  ((xlF))* 

Also,  given  any  pair  of  4-dimensional  vectors  |F)  =  (f,  F)  and 
|G)  =  (g,  C),  their  Hilbert  space  scalar  product  is  defined 
as 

(FIG)  =  ///(fg-f  F*-G)dx 

With  these  Hilbert  space  representations,  the  proof  of  the 
inverse  scattering  theorem  becomes  rather  straightforward. 
I.et  |Jw)  be  any  4 -dimensional  vector,  and  |Aw>  be  the  4- 
dimensional  vector  derivable  from  |  Jw)  in  the  same  way  as 
I  A,)  is  derivable  from  (J,>.  That  is, 

KmlA,)  = -  — |J,)  (48) 

c 


(A„(FJ  =  <Jw|Ks>  (4;fi 

which  is  the  Hilbert  space  form  of  the  inverse  scattering 
theorem.  Take  |Jw>  to  be  (0.  J«,(x)),  so  that  the  scalar  com 
ponent  of  ( A* )  also  vanishes,  then,  with  the  aid  of  Eq .  ( 47 ) . 
the  above  equation  gives  the  inverse  .scattering  theorem 
stated  in  the  previous  section,  viz.,  Eq.  (,30). 

One  may  wish  to  express  Eq.  (49)  in  a  more  s>  mmelni 
form,  such  as  one  involves  |<))  and  |E)  in  both  sides  of  the 
equality.  To  do  this,  we  first  replace  |  F,)  in  the  left  hand  side 
of  t)q.  (49)  by  SmlJ,).  as  defined  in  Eq.  (47).  then  Eq.  (49) 
becomes 

<Sw|Sm|.J,)  =  (Jw|R,>  loOl 

If  we  also  define  a  weighing  electric  field  |F)„  )  in  the  same 
way  a.s  a  physical  electric  field  is  related  to  the  vector  po 
tential  through  Eq.  (4.'>),  then  |F)w)  =  Sm|.^w).  Taken  its 
complex  conjugation  followed  by  hermitian  conjugation,  this 
relationship  gives  <Aw((Sm)"*  =  (Ew|.  From  Eq.  (44)one,seps 
that  (Sm)*’  differs  from  S^,  by  having  opposite  signs  in  all 
off-diagonal  elements  while  being  equal  to  Sn,  in  all  diagonal 
elements,  therefore  the  left  hand  side  of  Eq.  (.50)  is  not  equal 
to  (E„|  Js).  as  what  would  have  been  expected  from  a  sym 
metric  expression.  Therefore  we  define  another  field  |E* 

=  (Sm)*’|A«,>,  so  that  its  hermitian  conjugation  followed  by 
complex  conjugation  is  (6*1  =  (Aw|Sm-  Then  Eq.  (oO)  may 
also  l)e  stated  in  the  following  form; 

(E*|J,>  =  (J*|E,>  (oU 

It  is  remarked  that,  while  |A„)  may  be  considered  as  the 
potential  due  to  the  source  current  density  |.I*),  neither 
|E*)  nor  (E*  )  defined  above  has  a  parallel  analogy.  From  Fap 
(43),  a  parallel  analogy  would  require  that  the  scalar  <  om 
ponents  of  (E*)  and  |E*>  vanish.  But,  since  |.4w)  may  not 
satisfy  the  Lorentz  gauge  condition,  the  scalar  components 
of  Sml  Aw)  and  (Sn,)'’|  A*)  do  not  vanish. 

Equation  (51)  appears  to  be  similar  to  the  Lorentz  reci 
procity  theorem  as  expressed  by  Carson  (19|-(21|.  The  dif¬ 
ference  here  is  that  jJ*)  and  jE*)  are  only  weighing  func¬ 
tions  and  they  need  not  be  physical  quantities.  Therefore 
this  theorem  may  be  considered  as  a  generaliz-ation  of  t  he 
reciprocity  theorem.  The  theorem  is  valid  only  if  a  Hilbert 
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Figure  6a 


Figs.  6a-6b.  Comparison  of  the  electric  fields  in  the  z  =  7  cm 
surface  generated  by  a  planar  array  and  a  volumetric  array,  with 
the  phase  and  amplitude  conjugation  fi>cusing  on  the  axis  at  7  cm 
from  the  center  of  the  array,  (a)  For  a  system  of.')  parallel  arrays 
spaced  at  I  cm  from  each  other,  the  3  dB  main-beam  width  and 


the  mainbeam  to  first  sidelobe  ratio  in  the  (x,  y)  directions  are 
respectively,  <0.56, 0.82)cm  and  (24,  )8)dB.  (b)  For  a  planar  array, 
the  ,3  dB  main-beam  width  and  the  mainbeam  to  first  sidelobe 
ratio  in  the  (x,  y)  directions  are  respectively,  (0.64,  0.82)cm  and 
(10,  16)dB. 


space  may  be  defined  in  which  the  fields,  the  sources,  and 
the  weighing  functions  are  Hilbert  space  vectors,  which  is 
the  ca.se  when  the  propagation  medium  is  dissipative. 


7.  REMARKS  AND  CONCLUSIONS 

Three  main  topics  have  been  presented  in  this  article;  ( 1 )  an 
instrumental  and  hardware  description  of  a  microwave  array 
system  being  developed  for  medical  application.  (2)  the 
pha.se  and  amplitude  conjugation  technique  which  may  be 
used  to  actively  ftKus  a  transmitting  microwave  array  or, 
with  the  help  of  the  inverse  scattering  theorem,  to  passively 
focus  a  receiving  array;  in  either  case,  it  may  achieve  a  .3- 
dimensional  focusing  in  the  neighborhood  of  the  array,  and 
(3)  an  inverse  scattering  theorem  on  retrieving  information 
of  a  scattering  target  from  limited  data  of  the  scattered  field, 
which,  if  applied  in  conjunction  with  the  phase  and  ampli¬ 
tude  conjugations,  may  provide  an  optimal  3-dimensional 
imaging  from  one  sampling  of  the  scattering  field  measured 
by  a  receiving  array. 

On  the  quality  of  3-dimensional  focusing  using  the  phase 
and  amplitude  conjugations,  we  have  presented  several  re¬ 
sults  graphically.  Summarizing  these  results,  the  following 
conclusions  may  be  drawn: 

1 .  Osing  the  method  of  phase  and  amplitude  conjugations, 
one  may  achieve  a  satisfactory  degree  of  3-dimensional  fo 


cusing  in  the  neighborhood  of  radiation  sources  in  a  lattice 
structure.  There  is  a  slight  shift  of  the  peak  point  of  the  field 
from  the  intended  focal  point  as  defined  in  the  phase-am¬ 
plitude  conjugation  factor.  The  shift,  as  outlined  in  Table 
1,  generally  points  toward  the  center  of  the  array.  This  should 
not  pose  any  problem  for  practical  applications  since  it  can 
be  calibrated. 

2.  Upon  applying  the  phase  and  amplitude  conjugations, 
the  field  patterns  and  the  beam  characteristics  in  the 
transverse  direction  appear  to  have  similar  dependency  upon 
the  lattice  structure  and  the  array  size  as  that  of  a  Fraun 
hofer  field.  'I'hus,  the  transverse  beamwidth  becomes  nar¬ 
rower  as  the  element  spacing  increases,  at  the  expense  of 
more  grating  lobes.  Interestingly,  this  behavior  also  applies 
to  the  longitudinal  beamwidth.  Therefore,  for  a  smaller 
target,  it  is  possible  to  improve  the  resolution  further  by 
increasing  the  element  spacing,  as  long  as  the  target  does  not 
extend  to  the  region  covered  by  grating  lobes.  Alternatively, 
the  regions  covered  by  the  grating  lobes  may  be  excluded 
from  illumination  by  active  control  on  the  transmitting 
array. 

3.  Along  the  longitudinal  direction,  the  field  patterns  and 
the  beam  characteristics  using  the  phase  and  amplitude 
conjugations  differ  considerably  from  that  of  a  Fraunhofer 
field.  A  Fraunhofer  field  is  invariant  in  the  longitudinal  di¬ 
rection  except  for  the  inverse-square  dependence,  whereas 
a  local  field  under  phase  and  amplitude  conjugations  has  a 
diffraction  structure  in  the  longitudinal  direction,  as  well  as 
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in  the  transverse  direction.  Therefore  some  degree  of  fo¬ 
cusing  of  the  local  field  along  the  longitudinal  direction  may 
be  achieved  using  the  phase  and  amplitude  conjugation 
technique.  It  must  be  remarked  that,  as  the  axial  distance 
increases,  the  sensitivity  of  longitudinal  focusing  to  the 
phase-amplitude  conjugation  decreases.  Our  analysis  of  the 
field  pattern  for  focusing  at  30  cm  axial  distance  indicates 
that  the  conjugation  factor  is  totally  overcome  by  the  ex¬ 
ponential  attenuation.  However,  the  transverse  focusing 
remains  good  even  at  this  axial  distance. 

4.  If  all  array  elements  may  be  represented  by  parallel 
dipoles,  then  clearly  the  vector  potential  everywhere  must 
be  polarized  in  the  dipole  direction.  However,  in  the  local 
region,  the  electric  field  will  still  have  strong  polarization 
dependency.  This  polarization  dependency  is  reduced  as  the 
axial  distance  increases,  and,  at  an  axial  distance  of  30  cm 
and  beyond,  the  electric  field  is  highly  polarized  in  the  di¬ 
rection  of  the  source  dipoles  (I3i. 

The  3-dimcrisional  resolution  and  the  allowable  target 
volume  may  further  be  optimized  by  varying  the  lattice 
structure  and  the  element  spacing  of  the  receiving  array. 
However,  based  on  information  theory,  there  is  a  theoretical 
limit  on  the  ratio  of  the  target  volume  and  the  resolution.  If 
the  total  number  of  the  array  elements  is  N.  then  each 
sampling  of  the  scattered  field  provides  an  information 
equivalent  to  2N  real  numbers,  where  the  factor  2  accounts 
for  the  measurements  of  both  phase  and  amplitude.  The 
value  2N  is  the  maximum  information  one  may  expect  from 
the  image  reconstruction.  The  theoretical  limit  of  the  3- 
dimensional  resolution  from  each  sampling  of  the  scattered 
field  may  then  be  described  by  the  equation: 


volume  of  the  target 

- TTTT — 

volume  of  the  focal  region 

Therefore,  if  the  target  is  smaller,  the  resolution  may  be 
further  improved  without  increasing  the  number  of  ele¬ 
ments.  For  example,  for  the  enlarged  interelement  spacings 
corresponding  to  Figs.  2c  and  3c,  the  resolution  is  improved 
to  about  half- wavelength  transversely  and  one-wavelength 
longitudinally. 

It  is  possible  to  acquire  more  target  information,  and 
thereby  improve  the  three-dimensional  resolution,  by 
making  multiple  views  from  different  angles  with  respect  to 
the  direction  of  the  incident  wave  as  suggested  by  the  models 
of  diffraction  tomography.  Multiple  views  at  different  ranges 
from  the  target  for  a  fixed  transmitter  position  are  also 
possible;  however,  analysis  of  such  a  volumetric  synthesis 
due  to  the  superposition  of  parallel  planar  arrays  discloses 
that  such  an  approach  may  be  of  limited  value,  because  much 
of  the  information  contained  in  parallel  samples  of  the 
scattered  field  are  redundant.  Figure  6(a)  shows  the  field 
pattern  of  five  parallel  arrays  separated  at  1  cm  from  one 
another.  Comparing  to  the  field  pattern  of  a  single  array  as 
shown  in  Fig.  6(bl,  the  five-array  system  provides  slightK’ 
narrower  main-beam  width  and  smaller  grating  lobes. 
However,  the  difference  may  not  repay  the  minimum  of  five 
times  the  data  acquisition  and  data  processing  time. 
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